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Abstract 

A  dual  fine  tracking  control  system  (FTCS)  is  developed  for  a  single  aperture 
optical  communication  receiver  to  compensate  for  high  frequency  disturbances  affect¬ 
ing  tracking  of  two  incident  laser  communication  beams.  The  optical  communication 
receiver  resides  within  a  payload  module  aboard  a  geosynchronous  satellite,  while 
each  laser  communication  transmitter  is  housed  within  a  module  aboard  a  high  alti¬ 
tude  unmanned  aerial  vehicle  (UAV).  In  addition  to  platform  specific  disturbances, 
the  impact  of  atmospheric  optical  disturbances  affecting  tracking  error  are  investi¬ 
gated.  The  system  dynamics  and  FTCS  are  modeled  and  evaluated  in  MATLAB 
and  SIMULINK.  An  optimal  controller  is  developed  to  mitigate  these  disturbances 
and  provide  tracking  errors  commensurate  with  a  bit  error  rate  (BER)  that  does  not 
exceed  10-6. 

Based  on  the  respective  optimal  state  estimates  of  each  beam,  the  dual  control 
technique  regulates  the  fine  tracking  error  for  each  beam  by  switching  in  time  between 
each  state  estimate  and  applying  linear  quadratic  regulator  (LQR)  control  to  a  single 
fast  steering  mirror  (FSM).  Optical  turbulence  induced  random  beam  wander  effects 
revealed  limited  impact  to  tracking  error  due  to  spatial  separation  between  transmit¬ 
ters,  due  to  the  defined  communication  architecture  geometry.  Moreover,  simulation 
results  indicated  that  dual  control  did  not  achieve  the  tracking  error  limit;  however, 
single  control  of  one  beam  at  the  cost  of  tracking  error  on  the  other  beam,  achieved 
a  sufficient  temporal  mean  tracking  error  to  meet  the  required  BER. 
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Dual  Fine  Tracking  Control 


of  A 

Satellite  Laser  Communication  Uplink 

I.  Introduction 

Satellite  communication  systems  predominately  operate  within  the  radio  fre¬ 
quency  (RF)  and  microwave  frequency  regions  of  the  electromagnetic  spectrum.  Laser 
communication  systems  operate  at  higher  frequency  regions,  infrared  (IR)  to  ultravio¬ 
let  (UV),  which  are  capable  of  providing  more  bandwidth.  On  the  other  hand,  due  to 
narrow  beamwidth,  laser  communication  systems  are  more  sensitive  to  propagation 
disturbances.  Propagation  disturbances  cause  degradation  of  incident  beam  energy  at 
the  optical  receiver,  which  adversely  affects  communication  link  bit  error  rate  (BER) 
performance. 

The  fundamental  objective  of  this  work  is  to  design  an  optical  communication 
receiver  fine  tracking  control  system  to  compensate  for  disturbances  affecting  line  of 
sight  (LOS)  uplink  beam  propagation  from  two  spatially  separated  aircraft  terminals, 
as  shown  in  Figure  1.1.  The  design  methodology  encompasses  four  steps:  define  the 
physical  system  and  control  requirements,  develop  a  system  model,  design  a  controller, 
and  simulate  the  design.  If  successful,  this  method  would  provide  a  possible  solution 
as  a  multiple  access  scheme  for  satellite  laser  communications. 

1.1  Dual  Fine  Tracking  Control  System 

The  receiver  dual  fine  tracking  control  system  (FTCS)  is  shown  in  Figure  1.2. 
Each  incident  beam  passes  through  a  single  telescope  aperture  and  is  transferred  by 
a  fast  steering  mirror  (FSM),  denoted  by  F(t)  subject  to  plant  disturbance  uq(t),  to 
signal  and  position  photodetector  sensors.  The  position  sensor,  known  as  a  quadrant 
cell  detector  (QD),  converts  the  received  optical  power  into  azimuth  and  elevation 
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Satellite 


Figure  1.1:  Laser  Communication  System  Architecture 
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Figure  1.2:  Dual  Fine  Tracking  Control  System 


tracking  error  signals,  and  is  subject  to  sensor  noise  v(t).  Within  the  optical  commu¬ 
nication  system  (OCS),  a  signal  sensor  converts  received  optical  power  into  electrical 
signals  for  signal  processing.  Based  on  the  measured  states  for  each  beam,  m^f)  and 
msit)  respectively,  the  dual  control  technique  regulates  the  tracking  error  for  each 
beam  by  switching  in  time  between  optimal  estimates  of  each  system  state  and  com¬ 
mands,  via  u(t),  a  single  FSM  to  center  the  incident  beam  on  the  position  detector 
for  maximum  available  received  signal  power.  Tracking  performance,  depending  on 
internal  and  external  disturbance  rejection,  is  quantified  by  maintaining  the  residual 
tracking  error,  Er(t),  within  the  fine  tracking  error  budget. 

1.1.1  Internal  Disturbances.  Control  system  internal  disturbances  are  caused 
by  sensor  noise,  sensor  bias,  actuator  nonlinearities,  and  satellite  plant  disturbances. 
Sensor  noise  is  caused  by  uncertainties  in  measured  states  due  to  quantum  shot  noise, 
Nq,  and  thermal  noise,  N0.  Sensor  bias,  unique  to  each  sensor,  is  a  constant  error 
between  the  actual  and  measured  state.  In  particular,  dark  current,  Nd,  is  a  bias 
inherent  to  a  photodetection  devices.  Plant  disturbances  are  vibrations  caused  by 
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Figure  1.3:  Uplink  Communication  Block  Diagram 


platform  structural  loads.  These  internal  disturbances,  categorized  as  receiver  effects 
as  shown  in  Figure  1.3,  induce  beam  tracking  loss,  denoted  as  Lr. 

1.1.2  External  Disturbances.  Uplink  propagation  of  each  beam  is  adversely 
affected  by  external  disturbances,  which  attenuate  the  transmitted  power,  Pt.  Subse¬ 
quent  to  propagation,  the  observed  received  power,  Pr,  determines  link  margin  and 
BER  performance,  which  are  fundamental  performance  measures  in  digital  commu¬ 
nication  systems.  These  external  disturbances  can  be  categorized  as  transmitter  and 
channel  effects,  as  shown  in  Figure  1.3,  and  described  below. 

1. 1.2.1  Traiismitter  Effects.  Transmission  of  a  signal  through  a  chan¬ 
nel,  the  propagation  medium,  requires  an  antenna  sized  to  the  signal  wavelength  [29, 
page  168].  Wavelength  is  inversely  proportional  to  frequency;  hence,  laser  commu¬ 
nication  yields  a  small  antenna  size  and  narrow  beamwidth.  Although  a  narrow 
beamwidth,  on  the  order  of  10  /irad,  may  reduce  interference  effects,  any  mechanical 
disturbance  adversely  affects  pointing  accuracy  [15,  page  137].  Typical  pointing  loss 
sources,  Ltl  are  aircraft  flight  structural  loads  arising  from  control  hardware,  tur¬ 
bulence,  and  unsteady  aerodynamics.  Unlike  RF  systems,  laser  communication  sys¬ 
tems  encounter  additional  signal  degradation  due  to  aerodynamic-optical  loss,  Lao. 
Aerodynamic-optical  loss  is  caused  by  beam  phase  distortion,  due  to  turbulent  bound¬ 
ary  layer  effects,  during  transmission  from  an  airborne  platform. 
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1. 1.2.2  Channel  Effects.  The  communication  channel  consists  of  the 
Earth’s  atmosphere  and  free  space.  Propagation  through  the  atmospheric  portion 
produces  three  main  disturbances:  absorption,  scattering,  and  optical  turbulence  [2, 
page  1].  The  atmosphere  is  a  nonhomogenous  medium  in  which  atoms  interact  with 
an  incident  electromagnetic  wave,  depending  upon  its  wavelength,  by  either  dissipa¬ 
tive  absorption  or  scattering  [13,  pages  86-87].  Optical  turbulence  produces  index 
of  refraction  fluctuations  which  induce  beam  spreading,  random  beam  wander,  and 
beam  scintillation  upon  a  propagating  optical  wave.  Based  on  relative  atmospheric 
conditions,  these  disturbances  yield  an  overall  atmospheric  loss,  La. 

As  a  function  of  propagation  distance,  zp,  and  wavelength,  free-space  diffrac¬ 
tion  affects  an  electromagnetic  wave  by  increasing  its  initial  beam  waist  radius,  W0] 
to  a  final  waist  radius,  W(zp).  In  turn,  the  initial  beam  peak  irradiance,  IQ  is  dis¬ 
persed  about  the  final  waist  radius  to  a  decreased  final  peak  irradiance,  I0(zp).  As  a 
consequence,  the  attenuation  in  transmitted  power  is  defined  as  a  path  loss,  Ls.  In 
addition,  optical  background  radiation,  Af,,  affects  photodetection  sensitivity  as  an 
independent  parameter  in  quantum  shot  noise  power. 

1.1.3  Fine  Tracking  Error  Budget.  The  fine  tracking  error  budget  is  com¬ 
prised  of  error  sources  arranged  in  a  hierarchical  relationship,  as  shown  in  Figure  1.4. 
Disturbance  error  contributors  are  external  and  internal  disturbance  sources,  as  de¬ 
fined  above.  Hardware  error  contributors  are  manufacture  and  assembly  of  system 
components,  which  are  assumed  to  be  corrected  or  minimized  during  analysis,  align¬ 
ment,  and  calibration.  Software  error  contributors  are  digital  communication  jitter, 
estimation  error,  and  control  system  sampling  jitter.  Sources  of  software  error  are 
assumed  to  be  corrected  or  minimized  during  real-time  or  off-line  calculations  applied 
to  the  closed-loop  control  algorithm.  Unforeseen  error  sources,  arising  during  mod¬ 
eling,  simulation,  or  experimentation,  will  be  analyzed  and  added  to  the  appropriate 
category. 
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Figure  1.4:  Fine  Tracking  Error  Budget  Tree 


1 . 2  Summary 

This  chapter  presented  the  dual  fine  tracking  control  problem,  which  investi¬ 
gates  the  feasibility  of  a  single  FSM  to  regulate  tracking  errors  for  two  uplink  beams. 
Chapter  II  provides  a  literature  review  of  common  methods  for  laser  beam  pointing 
and  tracking  control.  Chapter  III  presents  the  underlying  physics  of  laser  beam  up¬ 
link  propagation  and  development  of  the  dual  FTCS  model.  Chapter  IV  provides 
simulation  data  analysis  and  performance  results.  Lastly,  Chapter  V  summarizes  the 
research  and  contributions,  and  provides  recommendations  for  future  work. 
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II.  Literature  Review 


Maintaining  accurate  beam  pointing  in  the  presence  of  disturbances  has  been  a  fun¬ 
damental  goal  in  satellite  laser  communication  systems.  Typically  there  are  three 
types  of  pointing:  acquisition,  coarse  tracking,  and  fine  tracking.  Acquisition,  uti¬ 
lizing  satellite  ephemeris,  beacons,  or  navigation  data,  establishes  the  initial  line  of 
sight  (LOS)  communication  link  between  the  transmitter  and  receiver.  Coarse  track¬ 
ing  compensates  for  large  angular  deviations  and  low  frequency  disturbances,  while 
fine  tracking  compensates  for  high  frequency  disturbances.  These  functions  are  incor¬ 
porated  into  a  pointing,  acquisition,  and  tracking  (PAT)  control  system  to  maintain 
LOS  communications. 

Much  of  the  fine  tracking  control  research  focused  primarily  on  high  frequency 
disturbance  rejection,  also  known  as  jitter.  The  limiting  factor  is  closed-loop  con¬ 
trol  bandwidth  versus  disturbance  bandwidth  [30,  page  569].  Closed-loop  bandwidth 
is  determined  by  the  cutoff  frequency,  which  is  determined  by  the  -3  dB  point  of 
the  control  system  frequency  response  [21,  page  572].  Ideally,  2  to  3  KHz  control 
bandwidth  is  desired  [20,  page  72];  with  a  tracking  accuracy  of  1  p,rad,  due  to  the 
narrow  beamwidth  of  an  optical  link  [8,  page  252],  Several  methods  have  been  in¬ 
vestigated  for  high  frequency  disturbance  rejection:  platform  vibration  disturbance 
rejection,  acoustic  beam  steering,  nonlinear  optics,  and  antenna  gain  optimization.  A 
discussion  of  each  method  will  be  presented. 

2.1  Platform  Vibration  Disturbance  Rejection 

Tracking  errors  are  normally  caused  by  sensor  noise,  estimation  error,  noisy 
actuators,  and  platform  vibrations;  however,  platform  vibrations  are  the  dominant 
error  source  [8,  pages  252-253].  Platform  vibrations,  at  either  the  transmitter  or  re¬ 
ceiver,  produce  beam  position  fluctuations  on  the  receiver  position  detector.  On-orbit 
spacecraft  platform  vibrations,  characterized  by  measurement  data  from  the  European 
Space  Agency  (ESA)  OLYMPUS  satellite,  revealed  high  frequency  disturbances  above 
100  Hz  affecting  an  optical  communication  payload  [36,  pages  207-212],  In  particular, 
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power  in  the  frequency  responses  are  observed  approximately  up  to  500  Hz  from  solar 
array  actuators,  1  KHz  from  microwave  switches,  and  up  to  450  Hz  from  thruster 
firings. 

Through  various  studies,  comparison  of  linear  versus  nonlinear  tracking  con¬ 
trol  methods  yielded  varied  low  frequency  performance;  however,  nonlinear  methods 
achieved  better  high  frequency  performance.  Assuming  disturbances  are  zero  mean 
white  Gaussian  noise,  adaptive  direct  feedforward  control  achieved  approximately  9 
dB  greater  jitter  disturbance  rejection  than  linear  quadratic  Gaussian  (LQG)  control, 
at  the  cost  of  computation  resources  [25,  pages  191-194],  An  alternate  design,  incor¬ 
porating  a  self-tuning  feedforward  control  loop  with  linear  feedback  control,  achieved 
at  least  a  25  dB  disturbance  rejection,  and  without  self-tuning,  linear  feedback  con¬ 
trol  dropped  to  below  10  dB  of  disturbance  rejection  at  frequencies  above  300  Hz  [31, 
pages  572-574].  Investigation  of  frequency  dependent  noise,  commonly  referred  to  as 
colored  noise,  disturbance  effects  revealed  similar  high  frequency  disturbance  rejec¬ 
tion  results  for  LQG  control  compared  with  adaptive  least  mean  squares  (LMS)  and 
gradient  adaptive  lattice  (GAL)  algorithms,  however,  both  adaptive  control  methods 
resulted  in  added  noise  to  the  system  response  [33,  page  16]. 

2.2  Acoustic  Beam  Steering 

Acoustic  beam  steering  relies  on  the  acoustic  optic  phenomenon  using  sound 
waves  to  change  the  index  of  refraction  of  the  propagation  medium.  The  simplest 
form  of  this  acoustic  optical  effect  is  Bragg  diffraction:  an  incident  optical  beam  is 
diffracted  by  constructive  interference  based  on  [26,  pages  801-806] 


(2.1) 


8 


Bragg  Cell 


s 


t 

► 


k 


s 


Figure  2.1:  Bragg  Diffraction 


where  Ac  is  the  carrier  wavelength,  As  is  the  sound  wavelength,  and  9b  is  the  Bragg 
angle.  This  condition  is  equivalent  to  the  following  vector  relationship, 


(2.2) 


where  kt  is  the  transmitted  optical  wave,  ki  is  the  incident  optical  wave,  and  ks  is  the 
sound  wave.  The  vector  magnitudes  are  determined  by  wave  numbers  kc  =  2ir/Xc  and 
ks  =  2tt/\s,  for  optical  and  sound  wave  vectors  respectively.  An  optimum  amplitude 
reflectance  occurs  when  ks  =  2kc  sul(0b)  [26,  page  805]. 

A  common  device  used  to  implement  acoustic  beam  steering  is  a  Bragg  Cell, 
as  shown  in  Figure  2.1  [26,  page  801].  Utilizing  a  Bragg  Cell  and  a  position  sensing 
photodetector,  a  deflection  angle  can  be  calculated  as  [19,  page  2209] 


(2.3) 


where  n,  is  the  index  of  refraction.  Based  on  the  measured  deflection  angle,  used  as 
an  error  source,  the  sound  frequency  input  to  the  Bragg  Cell  is  adjusted,  which  steers 
the  optical  beam  in  the  desired  direction  for  zero  deflection. 
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2. 3  Nonlinear  Optics 


Linear  or  nonlinear  behavior  of  an  optical  wave  is  dependent  upon  the  prop¬ 
agation  medium  [26,  page  738].  Free-space  exhibits  properties  of  a  linear  medium, 
while  a  nonlinear  medium  affects  the  propagation  of  light  due  to  a  nonlinear  relation¬ 
ship  between  the  polarization  density  and  electric  field  [26,  page  739].  Based  on  the 
optical  phase  conjugation  principle,  wavefront  aberrations  can  be  removed  from  an 
incident  beam  [26,  page  761].  Using  nonlinear  optical  receiver  components,  an  adap¬ 
tive  laser  beam  tracking  scheme  reduced  incident  beam  angular  deflection,  resulting 
in  better  performance  than  classical  tracking  systems  against  milli-radian  amplitude 
disturbances  at  frequencies  >100  Hz  [9,  page  72], 

2-4  Antenna  Gain  Optimization 

Disturbances  on  an  incident  beam  degrade  the  received  signal  power,  which 
affects  the  communication  system’s  bit  error  rate  (BER)  performance.  Depending 
upon  the  tracking  error  statistics,  transmitter  or  receiver  antenna  gain  can  be  opti¬ 
mized  to  improve  system  performance  [8,  page  252],  Following  the  approach  of  Chen 
and  Gardner,  antenna  gain  optimization  method  relies  on  initial  determination  of  the 
theoretical  received  signal  power  and  corresponding  system  BER. 

2-4-1  Received  Signal  Power.  Received  signal  power,  Pri  is  generally  deter¬ 
mined  by  the  Friis  transmission  equation  [23,  page  110]. 


(2.4) 


where  Pt  is  the  transmitted  signal  power,  Gt  is  the  transmitter  gain,  Ls  is  the  path 
loss,  and  Gr  is  the  receiver  gain.  Path  loss,  depending  upon  carrier  wavelength,  Ac, 
results  from  divergence  of  an  electromagnetic  wave  as  it  propagates,  which  is  denoted 
by  [29,  page  254] 


(2.5) 
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where  zp  is  the  propagation  distance.  Accounting  for  additional  loss  factors,  similar 
to  RF  communications,  the  transmission  equation,  also  known  as  the  link  equation, 
can  be  modified  as  [15,  pages  94-105] 

Pr  —  ( VtPtGtLtLao)(LsLa)(LrGrr]r )  (2-6) 


where  r]t  is  transmitter  efficiency  loss,  Lt  is  transmitter  pointing  loss,  Lao  is  aerodynamic- 
optical  loss,  La  is  atmospheric  loss,  Lr  is  receiver  tracking  loss,  and  r/r  is  the  receiver 
efficiency  loss. 


2-4-2  Bit  Error  Rate  Performance.  Elevation  and  azimuth  angular  errors, 
denoted  by  x  and  y  respectively,  are  assumed  to  be  statistically  independent  Gaussian 
random  variables  [15,  page  100].  The  probability  density  functions  of  these  errors  can 
be  used  to  determine  signal  loss  and  corresponding  BER  performance  statistics  [27, 
page  296].  The  angular  error  joint  probability  density  function  is  expressed  as  [10, 
pages  101,  208] 


fxr(x,y ) 


O-TV  (Jx(7y 


exp 


lx 


v 

mx) 


2  al 


exp 


(y  ~  mvY 


(2.7) 


where  m*  and  of  are  respectively  the  corresponding  mean  and  variance  statistics 
of  the  errors.  Assuming  the  azimuth  and  elevation  errors  are  caused  by  combined 
disturbance  effects,  which  yield  similar  behavior  in  either  direction,  error  variances 
are  assumed  to  be  equivalent  [27,  page  297]. 

In  addition,  the  control  system  regulates  error  about  a  zero  set  point,  thus  mean 
error  is  assumed  to  be  zero.  Based  on  both  assumptions,  the  joint  probability  density 
function  can  be  modified  as, 


fxv(x,y ) 


~-( x2  +  y 2y 

2cr2 


(2.8) 
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Figure  2.2:  Coordinate  Transformation 


Performing  a  transformation  into  polar  coordinates,  as  shown  in  Figure  2.2,  a  radial 
error  probability  density  function  can  be  determined,  since  the  transformed  prob¬ 
ability  is  approximately  equivalent  to  the  original  [10,  pages  227-230].  Hence  the 
probability  density  functions  are  related  by, 

fxy(x,  y)dxdy  =  fR&(p ,  0)dA  (2.9) 


where 


x  =  pcos[6) 

y  =  psin($) 

2  2,2 
p  =  x  +  y 

A  =  area 

Matching  the  derivatives  to  the  transformation  variables,  a  relation  is  given  by  the 
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transformation  Jacobian  [10,  page  229] 


dxdy 

=  J(p,  9)  =  det  | 

f  dx 
dp 

dx  \ 

=  det 

cos^) 

— psin^) 

dA 

dy 

V  dP 

dy  J 
do  / 

sin(0) 

pcos(9) 

=  p  [cos2(0)  +  sin2(#)] 

=  p  (2.10) 

Substituting  Equations  2.8  and  2.10  into  Equation  2.9,  the  polar  coordinate  joint 
probability  density  function  is  expressed  as, 

fm(p,  9)  =  fxY(x,y)J(p,0 ) 

=  (2'n) 

Assuming  6  is  equally  likely  to  be  any  value  within  [0,27t],  the  polar  coordinate 
joint  probability  density  function  can  be  converted  to  a  marginal  probability  density 
function  over  p,  which  has  Rayleigh  statistics  defined  as  [10,  page  230] 

/2tt  /  2  A 

fm(p,  0)d8  =  exp  \~“^J  (2.12) 

A  general  expression  for  the  average  BER  for  a  communication  system  is  a  function 
of  the  signal  modulation  probability  of  bit  error,  Pmod  defined  in  Appendix  A,  and  an 
arbitrary  radial  error  probability  density  function  defined  as  [27,  page  298] 

poo 

BER  =  /  [PmodfR(p)}  dp  (2.13) 

Jo 

The  domain  of  integration  assumes  the  radial  tracking  error  resides  within  [0,  oo], 
based  on  the  Rayleigh  statistics.  However,  receiver  physical  dimensions  can  limit  the 
domain  to  reside  within  the  receiver  aperture  radius. 

2.4.3  Transmitter  and  Receiver  Gain  Optimization.  Minimizing  BER  as  the 
cost  function,  optimum  gain  can  be  determined  by  evaluating  Pmod  as  a  function  of 


13 


received  signal  power,  Pr.  Assuming  the  incident  beam  behaves  as  a  Gaussian  beam, 
transmitter  gain,  receiver  gain,  and  corresponding  transmitter  pointing  and  receiver 
tracking  losses  are  defined  as  [8,  page  253] 


(2.14) 


=  l^fj 


(  nPgpY 


(2.15) 


Lt  -  exp  (-Gtpf) 

Lr  —  exp  (  —  GrPr) 


(2.16) 


(2.17) 


where  Wa  is  the  initial  beam  waist  radius,  Dap  is  the  receiver  aperture  diameter,  pt 
is  the  transmitter  radial  pointing  error,  and  pr  is  the  receiver  radial  tracking  error. 
The  radial  pointing  and  tracking  errors,  subject  to  the  previously  derived  Rayleigh 
statistics,  are  defined  respectively  as 


(2.18) 


(2.19) 


Incorporating  Equations  2.6  and  2.18  into  Equation  2.13,  an  optimal  transmitter  gain 
is  determined  by  [8,  page  255] 


(2.20) 


Solving  Equations  2.6,  2.18,  and  2.20  simultaneously  for  Gt  and  Pr  results  in  pointing 
error  variance  being  inversely  proportional  to  maximum  transmitter  gain  but  pro¬ 
portional  to  minimum  received  power  [8,  page  255].  Utilizing  a  similar  method,  an 
optimum  receiver  gain  can  be  determined  by  solving  Equations  2.6,  2.19,  and  2.20 
simultaneously  for  Gr  and  Pr. 
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2. 5  Summary 

Several  fine  tracking  control  methods  have  been  presented  for  an  incident  laser 
beam.  For  the  purposes  of  this  study,  the  fine  tracking  control  methodology  will  incor¬ 
porate  techniques  from  platform  vibration  disturbance  rejection  and  the  underlying 
theory  in  antenna  gain  optimization.  These  methods  were  chosen  to  limit  implemen¬ 
tation  complexity  and  take  advantage  of  using  a  single  fast  steering  mirror  (FSM) 
within  the  dual  fine  tracking  control  system  (FTCS)  described  in  Chapter  I.  Chap¬ 
ter  111  will  develop  the  dual  FTCS  model,  based  on  LQG  control,  and  apply  received 
signal  power,  radial  tracking  error,  receiver  tracking  loss,  and  BER  performance  as 
constraints  on  performance. 
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III.  Dual  Fine  Tracking  Control  System 

The  laser  communication  system  architecture,  illustrated  in  Figure  1.1,  consists  of 
two  communication  links,  one  from  each  aircraft  transmitter  to  the  satellite  receiver, 
commonly  referred  to  as  a  terminal-to-satellite  communication  uplink.  The  proposed 
method  developed  for  this  work  investigates  the  performance  of  the  receiver  terminal, 
using  a  single  fast  steering  mirror  (FSM)  providing  fine  tracking  control,  to  main¬ 
tain  line  of  sight  (LOS)  with  two  spatially  separated  transmitting  terminals.  The 
control  concept  relics  on  quadrant  detectors  measuring  relative  tracking  errors  for 
each  incident  beam,  while  the  FSM  corrects  for  these  errors  depending  on  controller 
commands.  Control  requirements  are  driven  by  the  maximum  allowed  radial  tracking 
error  to  maintain  a  communication  bit  error  rate  (BER)  of  10~6,  which  is  a  commonly 
used  performance  limit  in  satellite  communications  [1,  page  25]. 

The  satellite  is  assumed  to  be  a  three-axis  stabilized  rigid-body  platform,  derived 
in  Appendix  B,  located  in  a  geosynchronous  orbit  at  35,786  km  altitude  [16,  page  117]. 
The  receiver  terminal  module,  mounted  on  the  nadir  panel  and  aligned  with  the  yaw 
axis,  utilizes  a  Cassegrain  telescope  antenna  with  an  assumed  aperture  efficiency,  r]r, 
of  0.8  [23,  page  110].  The  telescope,  with  arbitrarily  selected  specifications  of  300  mm 
primary  mirror,  15  focal  ratio,  and  10  magnitude  ratio,  collects  all  incident  optical 
sources  within  its  field  of  view  (FOV)  [13,  pages  174,  221],  A  1.0  nm  optical  bandpass 
filter  is  used  to  limit  the  amount  of  incident  optical  background  radiation  [5,  page  680]. 
All  receiver  electronic  components  are  maintained  at  a  noise  temperature  of  400  K  [23, 
page  118]. 

Each  aircraft  platform  is  a  high  altitude  unmanned  aerial  vehicle  (UAV)  loitering 
within  the  satellite  receiver  field  of  view  (FOVr)  at  60,000  ft  (18.3  km)  altitude  above 
mean  sea  level  under  clear-sky  conditions.  The  transmitting  terminal  is  assumed  to 
be  mounted  above  the  aircraft  center  of  mass  (cm)  within  a  dorsal  module  aligned 
with  the  yaw  axis,  with  a  minimum  20  degree  elevation  angle.  For  the  purposes  of  this 
study,  each  UAV  is  assumed  to  be  a  three-axis  stabilized  rigid-body  platform  subject 
to  platform  disturbances  which  induce  pointing  errors  and  aerodynamic-optical  effects 
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on  the  uplink  beam.  The  transmitter  antenna  aperture  efficiency,  rjt,  is  assumed  to 
be  0.8  [23,  page  110]. 


3. 1  Transmitter 

3.1.1  Aircraft  Dynamics.  Based  on  the  geocentric  latitude  and  longitude 
of  each  aircraft,  (</>*,  Aj),  relative  to  the  GCI  frame  as  derived  in  Appendix  C,  the 
respective  position  and  velocity  vectors,  are  determined  by 


r,  =  rn 


cos(Aj)  cos(0j)i  +  sin(Aj)  cos +  sin (<f>i)k 


Vi  = 


dfj 

dt 


(3.1) 

(3.2) 


where  ra  =  18.3+.Re  =  6,  396  km  is  the  aircraft  position  magnitude,  Re  =  6, 378  km  is 
the  Earth  equatorial  radius,  and  c u;  is  the  angular  velocity  vector.  Attitude  dynamics 
for  each  aircraft,  as  derived  in  Appendix  C,  are  denoted  by 


—  -I. i  iUii  +  Ij  1Tj 


(3.3) 


where  Ij  is  the  aircraft  moment  of  inertia  matrix,  is  the  angular  velocity  vector 
skew  symmetric  matrix,  and  Tj  is  the  total  torque  vector. 

3.1.2  Laser  Source.  Each  aircraft  transmits  an  On-Off  Keying  (OOK)  band¬ 
pass  modulated  100  mW  near  infrared  (Near  IR)  laser  beam,  as  derived  in  Appen¬ 
dices  A  and  D,  with  an  initial  beam  radius  of  40  mm,  at  wavelengths  of  1510  nm 
and  1550  nm  respectively.  The  data  rate  from  each  transmitter  is  100  Kbps,  which 
equates  to  a  signal  bandwidth  of  200  KHz.  Based  upon  theoretical  uplink  beam  prop¬ 
agation  characteristics,  derived  in  Appendices  D  and  E,  the  initial  beam  radius  and 
transmission  wavelengths  were  chosen  to  maximize  transmission  gain  and  minimize 
the  adverse  effects  of  atmospheric  loss.  Assuming  a  lowest  order  transverse  elec¬ 
tromagnetic  (TEMqo)  collimated  Gaussian  beam,  as  derived  in  Appendices  A  and  D, 
each  transmitted  beam  can  be  described  by  its  optical  field  equation,  diffractive  beam 


17 


radius,  and  irradiance  profile  respectively  as 


(3.4) 

(3.5) 

(3.6) 

E0  =  signal  amplitude 
g(t)  =  modulation  baseband  waveform 
V  =  modulation  voltage  amplitude 
a ;c  =  carrier  frequency 
0C  =  carrier  phase 
W0  =  initial  beam  radius 
zp  =  propagation  distance 
kc  =  wave  number 

pz  =  radial  distance  from  propagation  axis 
Pt  =  transmitted  power 
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As  the  beam  propagates,  its  irradiance  profile  expands  based  on  the  diffractive 
beam  radius,  W(zp).  Hence,  the  received  power,  Pr,  is  calculated  from  the  received 
peak  irradiance,  Io(zp),  as  illustrated  in  Figure  3.1.  Based  on  the  FOVr  geometry 
defined  in  Appendix  D,  the  transmitter  elevation  angle  is  limited  from  77.3°  to  90°. 
Thus,  as  shown  in  Figure  3.2,  the  propagation  distance  decreases  as  elevation  angle 
increases  [16,  page  111].  Given  the  predefined  laser  source  parameters  and  elevation 
angle  range,  the  theoretical  uplink  beam  propagation  characteristics  are  shown  in  Fig- 
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Figure  3.1:  Beam  Propagation 


ures  3.3  and  3.4.  In  particular,  an  increasing  elevation  angle  results  in  less  diffraction 
and  higher  irradiance,  due  to  less  propagation  through  the  atmosphere. 


Figure  3.2:  Propagation  Distance  vs.  Elevation  Angle 


3.1.3  Transmission  Effects.  For  the  purposes  of  this  study,  emphasis  is 
placed  on  the  receiver  tracking  performance  and  it  is  assumed  each  UAV  is  a  rigid 
body  in  the  steady  flight  regime  and  loitering  within  FOVr.  Hence,  attitude  stabiliza- 
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Figure  3.3:  Uplink  Beam  Diffractive  Beam  Radius 


x  10'7 


Figure  3.4:  Uplink  Beam  Free-Space  Peak  Irradiance 
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tion  and  constant  flight  parameters  are  controlled  by  an  autopilot  [39,  page  307].  In 
addition,  since  the  transmitting  terminal  is  aboard  an  aircraft,  aerodynamic-optical  ef¬ 
fect  is  an  additional  phenomenon  which  degrades  the  initial  beam  propagation.  These 
adverse  transmission  effects,  defined  in  Appendix  C,  are  quantified  by  the  transmit¬ 
ter  pointing  loss,  Lt  =  —0.969  dB,  and  aerodynamic-optical  loss,  Lao  =  —1.549  dB. 
Based  on  the  associated  transmitter  antenna  gain,  this  pointing  loss  creates  indepen¬ 
dent  random  azimuth  and  elevation  pointing  error,  each  assumed  to  have  zero-mean 
white  Gaussian  statistics,  with  a  variance  calculated  in  Appendix  C  respectively  for 
each  beam  as 


(a2)  4  =  4.0274  x  10~12 

[rad2] 

(3.7) 

(of)  R  =  4.2436  x  10“12 

[rad2] 

(3.8) 

3.2  Channel 

As  illustrated  in  Figure  1.3,  the  propagation  channel,  comprised  of  the  Earth’s 
atmosphere  and  free-space,  produces  three  adverse  effects:  atmospheric  loss,  path  loss, 
and  optical  background  radiation.  As  defined  and  modeled  in  Appendix  E,  theoretical 
atmospheric  loss  and  path  loss  for  the  given  laser  source  parameters  are  illustrated 
in  Figures  3.5  and  3.6.  A  longer  wavelength  is  inherently  less  susceptible  to  optical 
turbulence;  however,  as  shown  in  Figure  3.7,  optical  turbulence  induces  random  beam 
wander  about  the  propagation  axis,  which  is  independent  of  wavelength  but  dependent 
on  elevation  angle.  Optical  background  radiation,  as  derived  in  Appendices  D  and  E, 
is  modeled  as  additive  white  gaussian  noise  (AWGN)  which  adversely  impacts  optical 
signal  detection. 
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Figure  3.5:  Atmospheric  Loss 


Figure  3.6:  Path  Loss 
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Figure  3.7:  Random  Beam  Wander 


3.3  Receiver 

3.3.1  Satellite  Dynamics.  As  shown  in  Appendix  B,  the  geocentric  inertial 
frame  (GCI)  satellite  position  and  velocity  vectors  are  determined  as 


=  RGErxyz 

(3.9) 

drs 

— —  =  cm  X  rs 
dt 

(3.10) 

where  rxyz  is  the  Earth  centered  fixed  (ECF)  frame  satellite  position,  RGE  is  the  ECF 
to  GCI  rotational  transformation,  and  ujs  is  the  satellite’s  angular  velocity  vector. 
Based  on  the  predefined  operational  parameters,  the  position  vector  magnitude  is 
rs  =  Re  +  Rh  =  42, 164  km,  where  Rh  is  the  satellite  altitude. 

3. 3. 1.1  Field  of  View.  The  satellite  FOV,  assuming  a  spherical  Earth, 
is  the  visible  circular  area  about  the  subsatellite  point  (SSP)  [16,  page  161].  Based 
on  the  operational  parameters  for  satellite  location  and  aircraft  transmitter  mini- 
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mum  elevation  angle,  as  calculated  in  Appendix  B,  the  satellite  FOV  is  determined 
to  be  FOVs  =  16.34°  [0.285  rad],  which  covers  a  circular  area  ±61.82°  [±1.078  rad] 
along  the  Earth’s  surface  about  the  SSP.  The  receiver  FOV,  as  calculated  in  Ap¬ 
pendix  D,  is  determined  to  be  FOVr  =  3.82°  [0.066  rad],  which  covers  a  circular  area 
±10.81°  [±0.189  rad].  The  relationship  between  satellite  FOV  and  receiver  FOV  is 
depicted  in  Figure  3.8. 


3.3. 1.2  Satellite  Attitude  Dynamics.  Assuming  small  attitude  motion, 
as  shown  in  Appendix  B,  satellite  orientation  is  defined  by  Enler  angles:  roll  angle, 
( t>s ,  pitch  angle,  6S,  and  yaw  angle,  ips.  Applying  attitude  kinematics,  the  satellite 
angular  velocity  vector  is  defined  by 


r  <4  -  O ]0lps 

LUS  = 

IVy 

= 

0S  id  0 

V  UJz  ) 

[  ips  ±  u0(ps  ) 

(3.11) 


where  is  the  orbital  angular  velocity.  However,  due  to  the  geosynchronous  orbit,  the 
orbital  period  is  equivalent  to  the  Earth’s  rotational  period,  which  is  u>e  =  7.289  x  10-5 
rad/s.  Hence,  satellite  attitude  dynamics  is  defined  by  linearized  Enler’s  equations  of 
rotational  motion  as, 


0s 

Pps 


LO  n 


T 

-L  II 


1± 


(4 


T~ 

-ido4>s  ±  y~ 

J-  Z 


'•Ps  +  ^O 


(4  ~  4) 
4 


(3.12) 


where  4  and  4  are  the  mass  moments  of  inertia  of  the  axisymmetric  satellite  body 
and  Tx,  Ty,  and  Tz  are  the  total  torques  about  each  principal  axis. 


3.3.2  Optical  Communication  System.  The  optical  communication  system 
(OCS),  as  derived  in  Appendix  A,  performs  OOK  demodulation  by  direct  detection 
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Figure  3.8:  Satellite  and  Receiver  Field  of  View 
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of  the  received  optical  signal  power,  Pr.  Received  power  is  used  to  determine  com¬ 
munication  link  BER,  a  fundamental  performance  measure  for  digital  communication 
systems.  As  a  function  of  radial  tracking  error,  Pr  is  defined  as 


Pr(Pr)  =  (arGr)Lr  =  (arGr)  exp  (—Grpl)  (3.13) 

where  pr  is  the  radial  tracking  error,  ar  incorporates  a  priori  known  parameters  of  the 
Friis  transmission  equation,  Gr  is  the  receiver  gain,  and  Lr  is  the  tracking  loss.  The 
received  optical  signal  power  is  converted  by  an  avalanche  photodiode  (APD),  with  a 
multiplicative  gain  of  20,  into  a  received  current  signal,  defined  as 

P  idc  T  in 

—  RphPr+in  [A]  (3-14) 

where  idc  is  the  detector  DC  current,  Rph  is  the  detector  responsivity,  and  in  is  the 
detector  noise.  After  filtering,  the  received  bit  energy  to  average  noise  power  spectral 
density  ratio  Eb/N0,  as  derived  in  Appendix  A,  is  determined  as 

^  =  [KrarGrexp  (-Grp2r)]2  (3.15) 

where  Kr  is  a  function  of  Rph  and  noise  power.  Based  on  the  previously  defined 
transmitter,  channel,  and  receiver  parameters,  the  relationship  between  Eb/N0  and 
pr  is  shown  in  Figure  3.9.  In  particular,  the  longer  wavelength,  less  affected  by 
atmospheric  loss,  has  a  slower  attenuation  rate  as  pr  increases. 

The  associated  probability  of  bit  error,  as  defined  in  Chapter  II,  is  determined 
by  [27,  page  298] 

/»oo 

PB  =  BER  =  /  {Prnod.fn{pr)\  dpr  (3.16) 

Jo 
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Figure  3.9:  Eb/N0  Performance  due  to  Radial  Tracking  Error 


where  Pmod  is  the  modulation  BER  and  fii(pr)  is  the  radial  tracking  error  probability 
distribution  function,  defined  in  Chapter  If  by  Equation  2.19.  After  incorporation  of 
the  OOK  modulation  BER,  defined  in  Appendix  A  by  Equation  A. 38,  the  probability 
of  bit  error  becomes 

/»oo 

PB=  Q[KrarGrexp(-Grp2r)]  fR(pr)dpr  (3.17) 

Jo 

where  Q[x]  is  the  complementary  error  function.  Based  on  the  previously  defined 
transmitter,  channel,  and  receiver  parameters,  the  theoretical  probability  of  bit  error 
as  a  function  of  radial  tracking  error  is  shown  in  Figure  3.10.  In  particular,  in  order 
to  maintain  a  BER  =  10-6  the  maximum  radial  tracking  error  for  Beam  A  (1510 
nm)  is  approximately  1.68  /rrad  and  Beam  B  (1550  nm)  is  approximately  1.74  yurad. 
Thus,  for  the  purposes  of  this  study,  1.68  yurad  will  be  the  tracking  error  performance 
limit. 
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Radial  Tracking  Error  [grad] 

Figure  3.10:  Probability  of  Bit  Error  Performance  Curve  due  to  Radial  Tracking 
Error 

3.3.3  Receiver  Effects. 

3.3.3. 1  Platform  Disturbances.  Since  attitude  disturbances  are  the 
primary  concern,  it  is  assumed  the  satellite  structure  and  associated  ADCS  adequately 
compensates  for  expected  environmental  disturbances.  To  account  for  common  on- 
orbit  internal  disturbances,  as  described  in  Appendix  B,  a  zero  mean  Gaussian  noise 
power  spectral  density  (PSD)  model  was  developed  by  the  European  Space  Agency 
(ESA)  OLYMPUS  satellite  experiment  as  [36,  page  206] 

Ss(f)  =  [/^ad 2/Hz\  (3.18) 

3. 3. 3. 2  Detector  Noise.  As  derived  in  Appendix  D,  photodetector 
noise  is  comprised  of  three  independent  noise  sources:  quantum  shot  noise,  thermal 
noise,  and  dark  current.  As  such,  the  equivalent  current  noise  power  within  a  pho- 


todetector  is  defined  as  [1,  page  154] 


il  =  i2qsh  +  i2th  +  2<leiduBWd  [A2]  (3.19) 

where  i2qsh  is  quantum  shot  noise  current  power,  i2h  is  thermal  noise  current  power,  qe 
is  the  electronic  charge,  idu  is  unmultiplied  dark  current.  BWd  is  the  detection  band¬ 
width,  which  is  a  function  of  the  null-to-null  signal  bandwidth,  BWnn  =  200  KHz, 
determined  as  [17,  pages  943-944] 


BWd  ~  0.886BWnn  =  (0.886)(200) 

«  177.2  [KHz]  (3.20) 

Thus,  as  shown  in  Appendix  D,  the  the  maximum  theoretical  detector  equivalent 
noise  power  occurs  within  a  range  of  3.9162  x  10“18  A2  to  3.9164  x  10“18  A2,  based 
on  wavelength  and  elevation  angle. 


3-4  Dual  Fine  Tracking  Control  System  Model 

Derived  from  Figure  1.2,  the  dual  fine  tracking  control  system  (FTCS)  model  is 
presented  in  the  closed  loop  diagram  shown  in  Figure  3.11.  Plant  disturbance  signals, 
wdA{t )  and  wdB(t),  are  caused  by  transmitter,  channel,  and  receiver  effects.  F(t)  is 
the  FSM  output  x-axis  and  y-axis  angle  used  to  steer  the  incident  beams.  Residual 
tracking  error,  for  each  beam,  is  determined  respectively  by 


ErA(t)  =  wdA(t)  -  F{t) 


ErB(t )  =  wdB(t )  -  F{t) 


wdAx(t )  -  Fx(t ) 

WdAy{t)  -  Fy(t) 

wdBx(t )  -  Fx(t) 

wdBy{t)  -  Fy(t) 


(3.21) 


(3.22) 


Each  quadrant  detector  (QD),  assumed  to  have  an  identity  gain  matrix  transfer  func¬ 
tion,  measures  the  respective  beam  centroid  azimuth  (x-axis)  and  elevation  (y-axis) 
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Figure  3.11:  Dual  Fine  Tracking  Control  System  Closed  Loop  Diagram 


tracking  error.  After  corruption  by  zero-mean  white  Gaussian  noise,  denoted  by  vA(t) 
and  vB(t)  respectively,  the  measured  states  are  determined  by 

ErAx(t )  +  VAx{t) 

ErAy(t )  +  vAy(t) 

ErBx(t )  +  uBa;(f) 

Er  By(t)  +  vBy(t ) 

A  Kalman  filter  is  used  to  estimate  the  system  states,  ahi(f)  and  xb(£)  respectively 
for  each  beam,  based  on  the  noise  corrupted  measured  states  and  control  signal,  u(t). 
Based  on  an  as  yet  to  be  determined  switching  time,  Tsw,  the  controller  switches 
between  the  estimated  system  states  and  generates  the  control  signal  based  on  the 
linear  quadratic  regulator  (LQR)  algorithm,  which  commands  the  FSM. 


(3.24) 


mB(t)  =  ErB(t)  +vB(t)  = 


mA{t)  =  Er  A(t)  +  vA(t) 


(3.23) 
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3-4-1  Plant  Disturbance.  Plant  disturbance  is  comprised  of  the  transmitter 
pointing  error,  satellite  internal  disturbances,  and  optical  turbulence  induced  random 
beam  wander.  Transmitter  point  error  and  beam  wander  are  assumed  to  be  zero-mean 
white  Gaussian  noise  inputs.  However,  satellite  disturbance  is  a  frequency  dependent 
zero-mean  Gaussian  distribution,  hence  nonwhite  noise,  also  known  as  colored  noise. 


3-4-1 -1  Colored  Noise  Source.  A  shaping  filter,  combined  with  a  unit 
spectral  density  white  noise  source,  is  used  to  model  colored  noise,  defined  as  [7, 
pages  78-79] 


Gsf{s) 


bF 

s  +  dp 


(3.25) 


where  the  filter  parameters,  Of  and  bp,  are  determined  from  the  filter  output  power 
spectral  density  as  [25,  page  190] 


Ssf(u) 


b-P  =  (  1  ^ 

ui2  +  a2F  \l  +  f2J 

(2vr)2 

UJ2  +  (27t)2 


(3.26) 


Thus,  the  shaping  filter  transfer  function  is  defined  as 


GSf{s) 


2tt 

s  +  2tt 


(3.27) 


The  associated  satellite  disturbance  white  noise  source,  defined  by  zero  mean  Gaussian 
statistics,  has  a  power  spectral  density  defined  by  [25,  page  190] 


Ss{f)  =  —  =  80  x  10-12  [rad2 /Hz] 


(3.28) 


3-4-1-2  White  Noise  Sources.  Transmitter  pointing  errors  and  optical 
turbulence  induced  random  beam  wander,  are  assumed  to  be  frequency  independent. 
Thus,  an  equivalent  white  noise  power  spectral  density  for  each  noise  source  is  deter- 
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mined  respectively  as 


[St(f)]A  =  {at)A  —  4-0274  x  iO-12  [rad2/ Hz]  (3.29) 

[St(f)]B  =  (of )B  =  4.2436  x  IO’12  [rad2/ Hz]  (3.30) 


Sbw(f)  =  [rms(bw)]2  [rad 2 /Hz\ 


(3.31) 


3-4- 1-3  Total  Plant  Disturbance.  The  total  plant  disturbance  white 
power  spectral  density,  for  azimuth  and  elevation  tracking  error,  is  defined  respectively 
for  each  incident  beam  as 


[SM)]A  =  ss(/)  +  [s,(/)]x  +  [sta(/)L, 

=  (80  +  4.0274)  x  10~12  +  [rmsfbw)^]2  (rad 2 /Hz]  (3.32) 

[S„(/)]B  =  Ss(f)(f)  +  {St(f)}B  +  {Saf)}B 

—  (80  +  4.2436)  x  10-12  +  [rms(bw)B]2  [rad2//fa]  (3.33) 

where  rms(bw)  for  each  beam  is  determined  by  its  spatial  location  within  the  receiver 
FOV.  Thus,  based  on  the  previously  determined  random  beam  wander  profile,  the 
satellite  disturbance  is  the  dominant  plant  disturbance  source. 

3-4-2  Measurement  Noise.  Measurement  noise  is  caused  by  the  QD  noise, 
assumed  to  be  zero-mean  white  Gaussian  noise.  However,  QD  SNR  is  inversely  pro¬ 
portional  to  the  azimuth  and  elevation  error  variance,  as  derived  in  Appendix  D  [15, 
page  182],  Thus,  based  on  the  predetermined  maximum  theoretical  detector  noise 
power  and  assuming  a  unit  signal  power,  the  measurement  noise  power  spectral  den- 
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sity,  is  defined  respectively  for  each  incident  beam  as 


[Sv{f)]A  «  3.9163  xlO"18 
[Sv(f)]B  «  3.9164  xlO”18 


[rad  2 /Hz] 

(3.34) 

[rad  2 /Hz] 

(3.35) 

3-4-3  Plant.  The  FSM  is  a  linear  mechanical  device  subject  to  rotational 
motion  about  two  orthogonal  axes,  commonly  referred  to  as  tip  (x-axis)  and  tilt  (y- 
axis)  axes.  Assuming  uncoupled  dynamics,  motion  about  either  axis  is  modeled  by  a 
second  order  transfer  function  as  [21,  page  226] 


G, 


r< 

VJyy 


s2  ±  2(u JnS  +  LU2 


(3.36) 


where  u)n  is  the  natural  frequency  and  (  is  the  damping  ratio.  For  the  purposes  of 
this  study,  each  axis  is  assumed  to  have  the  same  natural  frequency  and  damping 
ratio,  defined  respectively  as  ivn  =  5000  rad/s  and  (  =  0.5,  based  on  experiment  and 
analysis  of  a  Newport  Corporation  FSM  [33,  page  3].  Thus,  the  FSM  transfer  function 
is  defined  as 


Gfsm  — 


G 


XX 

0  G 


0 


yy 


(3.37) 


Since  (  <  1,  each  axis  behaves  like  a  single-input  to  single-output  (SISO)  under¬ 
damped  system.  Hence,  the  FSM  poles  are  complex  conjugates  and  lie  in  the  left-half 
s-plane,  calculated  as  [21,  pages  225-226] 


P!,P2  =  -2.5  ±4.3301?  (3.38) 

p3,p4  =  -2.5  ±4.3301?  (3.39) 
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Thus,  each  FSM  axis  is  a  stable  SISO  system  with  an  initial  oscillatory  transient 
response.  The  associated  damped  natural  frequency  is  determined  as  [21,  page  227] 


^ d  n  C~ 

=  (5000)  a/ 1  -  (0.5)2  =  4330.1  [rad/s] 


(3.40) 

(3.41) 


3.4.3. 1  Time  Response.  Applying  the  MATLAB  Control  System  tool¬ 
box  step  response  algorithm,  the  FSM  step  response  is  shown  in  Figure  3.12.  Transient 
response  characteristics  for  each  FSM  axis,  rise  time,  maximum  percent  overshoot, 
and  2%  settling  time,  were  calculated  respectively  as  [21,  pages  231-233] 


—  )  tan  1  (  4^-  ]  =  0.0003  [sec] 
Ud) 


(3.42) 


MP%  =  exp  [-ir(£un)/ud]  x  100%  =  16.3  %  (3.43) 


ts 


4 


0.0016  [sec] 


(3.44) 


3. 4-3. 2  Frequency  Response.  Since  the  tip  and  tilt  axes  are  uncoupled 
and  have  the  same  response  behavior,  applying  the  MATLAB  Control  System  toolbox 
Bode  algorithm,  the  FSM  frequency  response  for  one  axis  is  shown  in  Figure  3.13.  The 
FSM  bandwidth,  based  on  the  —3  dB  cutoff  frequency,  is  determined  approximately 
as 


FSMC  =  6.36  x  103  [rad/sec] 


(3.45) 


=  1.012  [KHz] 


(3.46) 
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X-axis  XY-axis  Cross-Coupling 


Time  [msec]  Time  [msec] 


Figure  3.12:  Fast  Steering  Mirror  Step  Response 


Figure  3.13:  Tip/Tilt  Axis  Frequency  Response 
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3. 4  -3. 3  State  Space  Model. 


The  plant  state  space  representation  is 


defined  as 


X(t)  =  AX(t)  +  Buu(t)  +  B  Wwd(t ) 
m{t)  =  CmX(t)+v(t) 


where  the  state  vector  is  defined  as 


m 


Fx(t )  Fx{t)  Fy(t ) 


(3.47) 

(3.48) 


(3.49) 


Thus,  the  control,  plant  disturbance,  and  measurement  noise  input  vectors  are  defined 
respectively  as 


r  i T 

=  Ux{t)  Uy(t ) 

(3.50) 

r  i  T 

wd(t )  =  wdx(t)  wdy(t ) 

(3.51) 

r  i  T 

V(t)  ~  Vx(t)  Vy(t) 

(3.52) 

Based  on  the  system  state  vector,  input  vectors,  and  FSM  transfer  function,  the  linear 
time- invariant  state  matrices  are  defined  as 


A  = 


f  0 

-2.5  x  107 
0 

V  0 


1 

-5000 

0 
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0 

0 

0 

-2.5  x  107 


B„  = 


l  0 

2.5  x  107 
0 

V  0 


0 
0 
0 
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0  ^ 
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3-4-4  Kalman  Filter.  Since  the  LQR  algorithm  requires  full  state  feed¬ 
back,  an  estimator  is  required  to  estimate  the  tip  and  tilt  slew  rate,  Fx(t)  and  Fy(t) 


respectively.  Based  on  the  measured  states,  subject  to  independent  plant  noise  and 
measurement  noise  being  white,  a  Kalman  filter  is  the  optimal  estimator  [7,  page  231]. 
However,  the  plant  disturbance  is  nonwhite  noise,  thus  the  plant  model  must  be  aug¬ 
mented  with  a  shaping  filter  state  model  [7,  page  266].  The  Kalman  filter  is  designed 
based  on  the  augmented  plant  model,  with  the  filter  gain  being  determined  using 
algorithms  available  in  MATLAB. 

3-4-4- 1  Augmented  Plant  State  Model.  The  shaping  filter  state  space 
model  is  defined  as 


Xp(t)  —  Ap’A _p(t)  T  B/.'ici (f) 
wd(t)  =  CpXF  (t) 


(3.53) 

(3.54) 


where  uq  (t)  is  a  unit  spectral  density  two  input  matrix,  and  the  shaping  filter  state 
vector  is  defined  as 


T 


Xp(t)  =  wdx(t )  wdy(t ) 


(3.55) 


The  resulting  shaping  filter  matrices  are  defined  as 
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Thus,  the  augmented  plant  model  is  defined  as 


XG(t)  =  AGXG(t)  +  BuGu(t)  +  BwGw1(t)  (3.56) 

m(t)  =  CmGXG(t)  +  v(t)  (3.57) 


where  the  augmented  state  vector  is  defined  as 


Xg  (t) 


x(t)  :  xF(t) 


T 


(3.58) 


Based  on  the  augmented  state  vector  and  input  vectors,  the  time-invariant  state 
matrices  are  defined  as 


A 

B.WC  f 

B  u 

04x2 

;  Bug  = 

;  Bwg  = 

02x4 

1 

< 

02x2 

- 1 

PQ 

CmG  =  (  -Cm  ;  CF  ) 

The  augmented  output  matrix  is  chosen  to  yield  a  tracking  error,  defined  as 


Er(t)  = 

{  -Cm 

CF  ) 

x(t) 

=  -CmX(t)  +  C  FXF(t) 

= 

wdx(t ) 

wdy(t ) 

I  I 

I~I~|  hq 

c-+» 

_  XF(t)  _ 

(3.59) 

S.4-4-%  Steady-State  Kalman  Filter.  Although  the  plant  2%  settling 
time  is  0.0016  seconds,  the  fine  tracking  control  system  (FTCS)  is  envisioned  to 
operate  over  a  longer  time  period.  Thus,  the  Kalman  filter  is  designed  to  operate 
over  a  time  period  longer  than  the  plant  settling  time.  Hence,  a  steady-state  Kalman 
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filter  will  be  used  to  estimate  the  states  for  each  beam,  which  is  nearly  optimal  when 
estimating  states  over  a  long  time  period  [7,  page  256] .  The  steady-state  Kalman  filter, 
defined  by  the  predictor/corrector  form,  is  implemented  for  each  beam  respectively 
as  [7,  page  257] 


%A{t) 

h{t) 


A GXA(t)  +  BuGtt(f)  +  La  mA(t)  -  Cm GXA(t) 
AGXB(t)  +  BuG«(f)  +  Lb  mB(t )  -  CmGiB(f) 


(3.60) 

(3.61) 


where  the  state  estimate  vectors  are  defined  respectively  as 


XA(t) 

XB(t) 


X(t) 

m 


wdAx(t) 

wdBx(t) 


WdAy(t) 

WdBy(t) 


T 

T 


(3.62) 

(3.63) 


La  and  LB  are  the  steady-state  Kalman  gains,  defined  respectively  by  [7,  page  242] 


La  =  [EeMUCm^S^1  (3.64) 

Lb  =  [Ee(oo)]BCmg(S„)g1  (3.65) 


where  [S6(cx))]a  and  [£e(oo)]s  are  the  steady-state  estimation  error  covariance  matri¬ 
ces  respectively  for  each  beam  and  (S„)a  and  (S.„) B  are  the  respective  positive  definite 
diagonal  measurement  noise  spectral  density  matrices.  Each  steady-state  estimation 
error  covariance  matrix  is  computed  by  solving  for  the  steady-state  solution  of  the 
algebraic  Riccati  equation  defined  as  [7,  pages  243] 

0  =  Ee(oo)Ag  +  AGEe(oo)  +  [Bwg]  [Sto]  [Bwg] 

-  Ee(oo)  [CmG]  [S"1]  [CmG]  Ee(oo)  (3.66) 

where  S„,  is  the  associated  positive  definite  diagonal  plant  disturbance  noise  spectral 
density  matrix. 
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3-4-5  Controller.  As  determined  for  the  Kalman  filter,  the  dual  FTCS  is 
designed  to  operate  over  a  longer  time  period  than  the  transient  response;  hence,  the 
LQR  controller  is  designed  as  a  steady-state  controller  [7,  page  196].  However,  the 
control  system  must  regulate  the  tracking  error  occurring  from  both  beams.  Thus, 
with  only  a  single  FSM,  the  controller  is  designed  to  switch  between  the  estimated 
states  from  each  filter. 


3. 4-5.1  Linear  Quadratic  Regulator  Control.  Since  the  linear  plant  is 
driven  by  random  conditions  and  random  disturbance  inputs,  a  stochastic  regulator 
is  the  optimal  controller,  which  is  designed  to  minimize  a  cost  function  defined  as  [7, 
pages  218,  221] 


Jsr  —  vA 


Er  (t)ZEr(t)  +  nT(f)Ru(f) 


L/o 


(3.67) 


where  Z  is  the  tracking  error  weighting  matrix  and  R  is  the  control  weighting  matrix, 
in  which  both  are  initially  assumed  to  be  identity  matrices.  However,  due  to  being 
in  steady-state,  the  stochastic  regulator  gain  equals  the  steady-state  LQR  gain  [7, 
page  221],  Thus,  in  terms  of  the  state  estimates,  the  controller  cost  function  is  defined 
as 


J  = 


-E 
2 

=  -E 
2 


L/o 


Er  (t)ZEr(t)  +  uT(t)Ku(t) 


CmGXG(f)  Z  CmGXG(f)  +  uT(t)Ru(t) 


=  \e{J  XG(f)  [CmgZCmG]  XG(t)  +  uT (t)Hu(t) 


=  -E 
2 


X£(f)QXG(f)+MT(f)Ru(f) 


L/o 


(3.68) 


where  Q  is  the  state  weighting  matrix.  Based  on  the  chosen  state  and  control  weight¬ 
ing  matrices,  the  optimal  steady-state  feedback  gain  is  determined  by  [7,  page  206] 


K  =  (R"1)  (Bug)  (P) 


(3.69) 
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where  P  is  the  steady-state  solution  of  the  algebraic  Riccati  equation,  which  is  defined 
as  [7,  pages  206] 


PAg  +  A^P  -  P  [Bug]  [RT1]  [Bug  P  +  Q  =  0  (3.70) 

Thus,  the  control  input  to  the  FSM  is  defined  as  [7,  page  187] 

u(t)  =  -K  ■  X(t)  (3.71) 

3. 4-5. 2  Switching  Time.  Since  the  dual  FTCS  regulates  the  tracking 
error  for  each  beam  by  switching  in  time  between  each  state  estimate,  the  system  can 
be  considered  a  switched  system.  Assuming  a  time  period  defined  by  an  arbitrary 
number,  Tf,  the  full-state  feedback  estimate  is  varied  based  on  a  time  interval  defined 
by 


t  e  [■ nTsw ,  (n  +  1  )TSW\  C  [0,  Tf\ 


(3.72) 


where  n  is  an  arbitrary  nonnegative  integer  and  Tsw  is  the  switching  time.  Thus,  the 
control  input  is  defined  as 


u{t) 


— K  •  X^(t)  :  n  —  1,3,5,... 
— K  ■  Xfi(t)  :  n  =  0,2,4,... 


(3.73) 


Although,  both  subsystems  are  linear  with  quadratic  cost  functions,  a  closed-form 
solution  for  an  optimal  switching  time  is  still  an  open  problem  [38,  page  2].  Ideally,  a 
switched  system  must  not  exhibit  the  Zeno  phenomenon,  defined  by  infinitely  many 
switchings  in  a  finite  time  period  [37,  page  2684] .  It  is  envisioned  during  the  analysis 
process  to  determine  if  the  dual  FTCS  satisfies  this  criteria. 
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3. 5  Summary 

This  chapter  developed  the  dual  fine  tracking  control  problem.  Based  on  the 
presented  theoretical  uplink  beam  propagation  characteristics,  a  maximum  residual 
radial  tracking  error  of  1.68  /irad  is  the  performance  requirement  to  achieve  a  commu¬ 
nication  BER  of  ICE6.  Since  the  plant  is  envisioned  to  operate  longer  than  its  settling 
time  and  is  driven  by  non- white  disturbance  noise,  a  steady-state  Kalman  filter  is  used 
to  estimate  the  plant  states  augmented  with  the  disturbance  noise  for  each  incident 
beam.  Based  on  these  estimated  augmented  states,  linear  plant,  and  quadratic  cost 
function,  a  steady-state  LQR  based  controller  is  desired.  However,  due  to  two  sets  of 
state  estimates,  the  controller  must  incorporate  some  type  of  switching  time  sequence 
to  regulate  tracking  error  on  both  beams.  Appendix  G  shows  the  implemented  dual 
FTCS  model  designed  in  SIMULINK. 
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IV.  Analysis  and  Results 

Analysis  of  FTCS  performance  is  based  on  applying  Monte  Carlo  analysis,  taking 
samples  generated  from  multiple  simulation  runs,  and  computing  an  ensemble  and 
temporal  average  of  the  outputs.  In  order  perform  ensemble  averaging,  the  system 
model  is  converted  to  a  fixed-step  system  using  the  SIMULINK  Runge-Kutta  solver 
solution.  After  trial  and  error,  a  fixed  time  step  of  Ts  —  5  x  10-5  sec  was  used  to  ensure 
the  solver  did  not  encounter  a  singularity.  Only  ten  Monte  Carlo  simulation  runs  were 
performed  to  ensure  against  creating  memory  overloads,  due  to  the  small  fixed  time 
step  and  final  time  selection,  in  processing  the  data  in  MATLAB.  In  addition,  a  rate 
transition  block,  defaulted  to  act  as  a  unit  delay  for  continuous  sample  time,  is  used 
on  all  plant  disturbance  and  measurement  noise  inputs  to  ensure  proper  timing  of 
these  inputs  into  the  system.  Correlation  time  for  for  the  bandlimited  white  noise 
blocks,  due  to  SIMULINK  requirements,  is  changed  to  a  integer  factor  of  the  fixed 
time  step  as  tc  =  20 Ts  =  0.001  sec. 

As  determined  in  Appendix  F,  simulation  results  reveal  optimum  regulator  per¬ 
formance  is  achieved  with  a  state  weighting  scale  factor  of  10  and  control  weighting 
scale  factor  of  0.1.  In  general,  initial  transients  occurred  at  intervals  of  approximately 
0.001  sec,  due  to  the  correlation  time  for  the  bandlimited  white  noise  blocks.  On 
the  other  hand,  neither  candidate  switching  time  case,  defined  in  Appendix  F,  posi¬ 
tively  or  negatively  affected  controller  performance.  Hence,  for  the  purposes  of  this 
study,  the  defined  optimum  regulator  weighting  parameters  will  be  implemented  with 
a  switching  time  of  0.02s,  which  allows  the  noise  correlation  time  to  remain  an  order 
of  magnitude  faster. 

In  order  to  determine  relative  capability,  dual  FTCS  performance  was  compared 
to  observed  tracking  error  without  control  compensation  and  single  fine  tracking  con¬ 
trol  of  either  beam.  Noise  effects  are  characterized  for  measurement  and  plant  distur¬ 
bance  inputs.  In  particular,  UAV  spatial  separation  affects  upon  plant  disturbance 
were  analyzed  against  the  observed  tracking  error.  Disturbance  rejection  analysis  was 
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performed,  for  either  control  case,  by  comparison  of  each  frequency  response  function 
(FRF). 

4-1  Noise  Effects 

Noise  affects  on  tracking  error  were  observed  for  each  axis,  based  on  a  5  second 
simulation  time.  Each  noise  source  was  modeled  by  a  SIMULINK  bandlimited  white 
noise  block  with  a  correlation  time  of  0.001  sec.  Since  the  satellite  plant  disturbance  is 
a  colored  noise  source,  a  shaping  filter  is  applied  after  the  associated  equivalent  white 
noise  source  block.  The  output  variance  of  each  block  was  scaled  to  a  unit  variance 
and  multiplied  by  a  corresponding  gain  to  yield  the  desired  noise  variance,  as  defined 
in  Chapter  III.  Each  block  is  assigned  a  unique  random  seed  value  which  is  reset 
every  simulation  case.  Frequency  responses  were  calculated  for  the  measurement  and 
plant  disturbance  noise. 

4-1.1  Measurement  Noise.  The  measurement  noise  time  response,  within 
a  1  second  time  domain,  is  shown  in  Figure  4.1.  In  particular,  the  modeled  sensor 
noise,  exhibits  similar  behavior  about  each  axis.  Based  on  the  time  response,  the  cor¬ 
responding  power  spectral  density  (PSD)  was  calculated  and  illustrated  in  Figure  4.2. 
Up  until  approximately  1  KHz  the  PSD  is  relatively  flat,  which  is  consistent  with  the 
chosen  noise  block  correlation  time. 

4-1.2  Plant  Disturbance.  As  defined  in  Chapter  III,  plant  disturbance  is 
a  function  of  the  transmitter  pointing  error,  satellite  disturbance,  and  optical  tur¬ 
bulence  induced  random  beam  wander.  Satellite  disturbance  is  modeled  as  colored 
noise,  while  pointing  error  and  beam  wander  effects  were  modeled  by  zero-mean  white 
Gaussian  statistics.  However,  as  determined  in  Appendix  E,  beam  wander  increases 
as  propagation  distance  increases,  which  is  a  function  of  decreasing  transmitter  eleva¬ 
tion  angle.  Elevation  angle  changes  based  upon  either  changing  transmitter  altitude 
or  Earth  central  angle  (ECA),  as  shown  in  Appendices  B  and  D.  Since,  transmission 
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Beam  B 


Figure  4.1:  Measurement  Noise  Time  Response 


Figure  4.2:  Measurement  Noise  Power  Spectral  Density 
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altitude  is  assumed  to  be  fixed,  plant  disturbance  will  vary  due  to  changing  EGA, 
which  manifests  as  a  changing  spatial  separation  between  aircraft  platforms. 

4- 1.2.1  Spatial  Separation  Effects.  Assuming  UAVa  is  located  at  the 
subsatellite  point  (SSP),  corresponding  to  an  elevation  angle  of  90°,  as  UAVb  closes 
spatial  separation  between  aircraft,  its  elevation  angle  increases,  as  determined  for  13 
arbitrarily  chosen  elevation  angles  within  the  receiver  field  of  view  (FOVr)  illustrated 
in  Figure  4.3.  In  particular,  the  corresponding  rms  beam  wander  decreases  as  the 
aircraft  elevation  angle  increases,  due  to  less  atmospheric  propagation. 
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Figure  4.3:  Elevation  Angle  Affects  on  Spatial  Separation  and  Beam  Wander 


Based  on  simulation  final  time,  Tf  —  1  second,  the  ensemble  average  mean  and 
standard  deviation  data  were  calculated  over  the  10  Monte  Carlo  runs  for  tracking 
error  and  control  input  subject  to  changing  rmsbw  resulting  from  each  arbitrarily 
chosen  elevation  angle.  Figures  4.4  thru  4.6  show  a  magnified  view  of  the  statistics 
on  the  time  response  data  for  each  of  the  13  elevation  angles.  Specifically,  Beam 


46 


10 


Beam  A 


Beam  B 


Figure  4.4:  Beam  Wander  Ensemble  Average  Mean  Tracking  Error  Response  for 

13  Elevation  Angles 

B  is  regulated  from  0  to  0.02  sec  resulting  in  reduced  mean  and  standard  deviation 
values.  Similarly,  Beam  A  is  regulated  from  0.02  to  0.04  sec  resulting  in  reduced 
mean  and  standard  deviation  values.  Control  input  response  shows  similar  behavior 
for  each  FSM  axis,  which  is  due  to  uncoupled  dynamics  and  each  axis  modeled  by  the 
same  SISO  system.  In  particular,  both  tracking  error  and  control  input  response  data 
exhibit  transients  0.001s,  due  to  noise  source  correlation  time,  which  behave  as  random 
amplitude  step  inputs.  The  data  indicates  similar  controlled  system  performance  at 
all  13  elevation  angles.  This  is  clearly  illustrated  by  taking  a  temporal  average  of 
the  ensemble  averaged  mean  and  standard  deviation  data  for  each  elevation  angle. 
Figures  4.7  thru  4.9  reveal  that  each  rmsbw  associated  with  each  elevation  angle  yields 
relatively  similar  temporal  mean  and  standard  deviation  results  for  each  beam.  In 
addition,  temporal  mean  and  standard  deviation  statistics  of  control  input  exhibits 
a  similar  trend.  Based  upon  these  results  for  the  range  in  elevation  angles  within 
FOVr,  spatial  separation  does  not  adversely  affect  FTCS  performance. 
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Beam  A 


Beam  B 


Figure  4.5:  Beam  Wander  Ensemble  Average  Tracking  Error  Standard  Deviation 

Response  for  13  Elevation  Angles 

X-axis  Y-axis 


Figure  4.6:  Beam  Wander  Ensemble  Average  Control  Input  Response  for  13  Ele¬ 
vation  Angles 
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Beam  A 


Beam  B 


Figure  4.7:  Beam  Wander  Temporal  Average  Mean  Tracking  Error  Response 
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X-axis  Y-axis 


Figure  4.9:  Beam  Wander  Temporal  Average  Control  Input  Response 

4- 1.2.2  Plant  Disturbance  Power  Spectral  Density.  Based  on  the  max¬ 
imum  theoretical  beam  wander,  plant  disturbance  noise  time  response  over  a  1  second 
time  frame  is  shown  in  Figure  4.10.  I11  particular,  the  modeled  plant  disturbance,  ex¬ 
hibits  similar  behavior  about  each  axis.  Based  on  the  time  response,  the  corresponding 
power  spectral  density  (PSD)  was  calculated  and  is  illustrated  in  Figure  4.11.  As  seen 
with  the  measurement  noise  response,  the  plant  disturbance  noise  PSD  decays  above 
approximately  1  KHz,  due  to  the  noise  block  correlation  time. 


4-2  No  Control 

Using  a  5  second  simulation  time,  tracking  error  data  was  observed  for  each  axis 
and  the  associated  radial  tracking  error,  as  derived  in  Chapter  II  as  p  =  Vx'2  +  V 2  ■ 
Figures  4.12  thru  4.13  illustrate  a  magnified  view  of  the  tracking  error  time  response 
ensemble  average  (10  runs)  mean  and  standard  deviation  data.  Clearly,  without 
regulation,  the  tracking  error  has  a  similar  time  response  as  the  plant  disturbance. 


50 


0  0.5  1  0  0.5  1 


Time  [sec]  Time  [sec] 

Figure  4.10:  Plant  Disturbance  Time  Response 


Beam  A  Beam  B 


Figure  4.11:  Plant  Disturbance  Power  Spectral  Density 
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Beam  A  Beam  B 
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Figure  4.12:  No  Control  Ensemble  Average  Mean  Tracking  Error  Time  Response 


In  terms  of  radial  tracking  error,  the  associated  ensemble  average  mean  and  standard 
deviation  is  shown  in  Figures  4.14. 


4-3  Single  Fine  Tracking  Control 

As  determined  in  Appendix  F,  optimum  regulator  performance  was  achieved 
with  a  state  weighting  scale  factor  of  10  and  control  weighting  scale  factor  of  0.1. 
Using  the  same  simulation  set  up  as  described  for  the  no  control  case,  radial  tracking 
error  time  response  data  was  observed  for  single  tracking  control  applied  for  each 
beam.  In  addition,  based  on  the  time  response,  the  FSM  frequency  response  was 
calculated  and  analyzed  for  disturbance  rejection  capability. 

4-3.1  Time  Response.  As  illustrated  in  Figures  4.15  and  4.16,  the  controller 
regulates  beam  tracking  error,  subject  to  random  amplitude  noise  inputs  at  intervals 
of  the  noise  correlation  time.  In  either  case,  there  is  no  apparent  benefit  in  regulating 
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Figure  4.13:  No  Control  Ensemble  Average  Tracking  Error  Standard  Deviation 

Time  Response 

Beam  A  Beam  B 


Figure  4.14:  No  Control  Ensemble  Average  Radial  Tracking  Error  Time  Response 
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Beam  A 


Beam  B 


Figure  4.15:  Beam  A  Control  Only:  Ensemble  Average  Radial  Tracking  Error  Time 
Response 


tracking  error  for  a  single  beam  over  the  other.  Any  variation  is  simply  due  to  the 
random  noise  samples. 


4-3.2  Frequency  Response.  Based  on  the  observed  simulation  time  response 
data,  each  signal  is  filtered  by  a  Hann  window,  processed  into  a  frequency  spectrum 
by  Fast  Fourier  transform  (FFT),  and  averaged  20  times  to  achieve  an  adequate  linear 
estimate  of  the  frequency  response  data.  As  illustrated  in  Figures  4.17  and  4.18,  the 
frequency  response  function  (FRF)  was  calculated  for  each  single  fine  tracking  control 
scheme  based  on  the  single  fine  tracking  control  input  and  corresponding  FSM  output 
response.  I11  either  case,  the  FSM  achieved  similar  magnitude  and  phase  distributions. 
Thus,  there  is  no  apparent  benefit  in  regulating  tracking  error  for  a  single  beam  over 
the  other. 


4-3.3  Disturbance  Rejection.  The  corresponding  disturbance  rejection  is 
calculated  and  illustrated  in  Figures  4.19  and  4.20.  Although  the  accuracy  of  this 
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Figure  4.16: 
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Beam  B  Control  Only:  Ensemble  Average  Radial  Tracking  Error  Time 
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Figure  4.17:  Beam  A  Control  Only: 

Function 


Fast  Steering  Mirror  Frequency  Response 
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Figure  4.18: 
Function 


Beam  B  Control  Only:  Fast  Steering  Mirror  Frequency  Response 


data  is  suspect  do  the  chosen  noise  correlation  time,  approximately  11.15  dB  power 
rejection  is  observed  for  either  regulated  beam  tracking  error.  Thus,  fine  tracking 
control  does  somewhat  reduce  line  of  sight  beam  propagation  effects. 


4-4  Dual  Fine  Tracking  Control 

Using  the  same  simulation  set  up  as  described  for  the  single  fine  tracking  control 
case,  a  switching  time  of  0.02s  was  implemented  for  dual  fine  tracking  control  of  both 
incident  beams.  The  ensemble  average  tracking  error  time  response  data  was  observed 
and  used  to  calculate  the  associated  FSM  frequency  response.  Based  on  the  observed 
plant  disturbance  frequency  response  spectrum,  disturbance  rejection  was  calculated 
for  each  beam. 

4-4-1  Time  Response.  As  illustrated  in  Figure  4.21,  the  controller  regulates 
radial  tracking  error  for  both  beams,  subject  to  random  amplitude  noise  inputs  at 
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Beam  A 


Beam  B 


Figure  4.19:  Beam  A  Control  Only:  Disturbance  Rejection 

Beam  A  Beam  B 


Figure  4.20:  Beam  B  Control  Only:  Disturbance  Rejection 
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Beam  A 


Beam  B 


Figure  4.21:  Dual  Fine  Tracking  Control  Ensemble  Average  Radial  Tracking  Error 
Time  Response 

intervals  of  the  noise  correlation  time.  Regulation  is  alternated  between  each  beam, 
based  on  the  chosen  switching  time  of  0.02s.  Comparison  with  either  single  track¬ 
ing  control  scheme,  Figures  4.15  and  4.16  respectively,  reveals  dual  tracking  control 
achieves  a  similar  regulation  performance,  sharing  the  negative  disturbance  effects. 

4-4-2  Frequency  Response.  The  frequency  response  function,  illustrated 
in  Figure  4.22,  was  calculated  for  the  dual  fine  tracking  control  scheme  using  the 
same  approach  discussed  in  section  4.3.2.  In  either  case,  the  FSM  achieved  similar 
magnitude  and  phase  distribution.  Although  highly  corrupted  by  noise  at  higher 
frequencies,  the  dual  control  FRF  exhibits  similar  behavior  as  single  fine  tracking 
control. 

4-4-3  Disturbance  Rejection.  The  corresponding  disturbance  rejection  is 
calculated  and  illustrated  in  Figure  4.23.  As  revealed  for  single  fine  tracking  control 
scheme,  the  accuracy  of  this  data  is  suspect  due  to  the  chosen  noise  correlation  time. 
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X-axis 


Y-axis 


Frequency  [Hz] 


Figure  4.22:  Dual  Fine  Tracking  Control  Control  Fast  Steering  Mirror  Frequency 
Response  Function 


Assuming  the  observed  frequency  response  is  accurate  up  to  a  1000  Hz,  slight  power 
rejection  is  observed  for  each  axis. 


4-5  Bit  Error  Rate  Performance 

Taking  a  temporal  average  of  the  ensemble  averaged  radial  tracking  error  data 
for  each  control  scheme,  illustrated  in  Figure  4.24,  yields  remarkably  different  results 
for  each  scheme.  Without  regulation  compensation,  the  observed  temporal  aver¬ 
age  mean  radial  tracking  error  for  both  beams  is  approximately  7.268  /irad.  While 
single  control  of  either  beam  yields  mean  error  of  approximately  1.6  /irad,  the  non- 
controllcd  beam  suffers  a  mean  error  of  approximately  10.11  /irad.  Dual  control 
yields  approximately  5.85  /irad  mean  tracking  error.  I11  addition,  standard  deviation 
statistics  exhibit  the  same  trend  between  each  control  scheme.  Both  statistics  reveal 
increased  dispersion  about  the  respective  means  for  the  regulated  beams  versus  the 
non-regulated  beams.  As  shown  in  Figure  3.10,  a  radial  tracking  error  greater  than 
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Figure  4.23:  Dual  Fine  Tracking  Control  Disturbance  Rejection 


approximately  4.2  /irad  yields  a  BER  of  1.  Thus,  only  regulation  of  a  single  beam, 
at  the  cost  of  the  other  beam,  satisfies  the  tracking  error  requirement  of  1.68  jurad. 
Chapter  V  will  discuss  a  possible  implementation  scheme  that  could  still  exploit  the 
dual  tracking  approach. 
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Beam  A  Beam  B 


Figure  4.24:  Temporal  Average  of  Ensemble  Radial  Tracking  Error  Performance 
Off:  No  Control,  A:  Beam  A  Control  Only,  B:  Beam  B  Control  Only,  A+B:  Dual 
Control 
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4-6  Summary 

Dual  fine  tracking  control  system  simulation  data  was  observed  and  compared 
to  non-control  and  single  fine  tracking  control  simulation  data.  Based  on  Monte  Carlo 
analysis,  ensemble  average  mean  and  standard  deviation  statistics  were  computed  for 
the  time  response  and  revealed  similar  tracking  error  regulation  performance  between 
single  and  dual  fine  tracking  control.  In  either  case,  regulation  occurred  at  regular 
time  intervals  due  to  the  noise  correlation  time.  However,  the  chosen  time  of  0.001s 
appeared  to  affect  the  frequency  response  of  the  disturbance  inputs.  Hence,  the 
control  system  frequency  response  accuracy  was  limited. 
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V.  Conclusions 


5. 1  Thesis  Summary 

A  conceptual  dual  fine  tracking  control  system  (FTCS)  was  designed  and  simu¬ 
lated  to  evaluate  the  effectiveness  of  a  single  fast  steering  mirror  (FSM)  to  compensate 
for  disturbances  affecting  line  of  sight  propagation  of  two  uplink  laser  communica¬ 
tion  beams.  The  optical  receiver  resided  in  a  satellite  in  geosynchronous  orbit  while 
each  incident  beam  originated  from  two  spatially  separated  unmanned  aerial  vehicles 
(UAVs).  Optical  frequencies  are  desired  for  their  low  power,  high  bandwidth,  and 
narrow  beamwidth  characteristics.  On  the  other  hand,  compared  to  radio  frequen¬ 
cies,  optical  frequencies  are  highly  susceptible  to  propagation  disturbances,  which 
attenuates  received  signal  power.  Degradation  in  received  signal  power  adversely  af¬ 
fects  communication  link  bit  error  rate  (BER)  performance,  which  is  a  fundamental 
performance  measure  in  digital  communications.  As  defined  in  Chapter  I,  for  the 
purposes  of  this  study,  the  primary  parameter  of  interest  for  disturbances  affecting 
uplink  beam  propagation  is  residual  tracking  error. 

Chapter  II,  provided  a  review  of  previous  works  covering  various  methods  for 
laser  beam  pointing  and  tracking  control.  In  particular,  optical  power  loss  can  be 
modeled  as  a  decreasing  exponential  function  of  radial  tracking  error  and  BER  per¬ 
formance  is  subject  to  tracking  error  probability  distribution.  Given  this  observation, 
a  dual  fine  tracking  control  system  was  designed  in  Chapter  III  based  on  a  residual 
radial  tracking  error  limit  of  1.68  /urad,  which  yields  a  BER  of  approximately  10~6. 
This  error  limit  was  determined  from  the  theoretical  propagation  characteristics  of 
each  uplink  beam. 

Propagation  of  each  beam  was  affected  by  transmitter  effects,  channel  effects, 
and  receiver  effects.  Other  than  the  power  loss  subject  to  propagation  distance,  the 
greatest  power  attenuation  occurred  from  atmospheric  loss.  To  minimize  atmospheric 
loss,  each  transmitted  beam  was  modeled  with  a  wavelength  of  1510  and  1550  nm 
respectively.  In  particular,  the  largest  source  of  atmospheric  loss,  approximately 
between  -4.5  to  -5  dB,  was  attributed  to  random  beam  wander,  which  is  caused 
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by  optical  turbulence.  Random  beam  wander  varied  as  a  function  of  propagation 
distance,  which  was  determined  from  the  elevation  angle  of  the  UAV  transmitter. 
However,  the  defined  satellite  receiver  field  of  view  (FOVr)  and  UAV  altitudes  limited 
the  range  of  elevation  angles  between  76  and  90  degrees.  This  range  in  elevation  angle, 
corresponded  to  UAV  spatial  separation  limited  to  ±10.81°  latitude/longitude. 

Chapter  IV  characterized  and  analyzed  dual  FTCS  performance.  Based  on  the 
respective  Kalman  filter  state  estimates  of  each  beam,  a  linear  quadratic  regulator 
(LQR)  controller  applied  steady  state  regulator  gain  and  switching  time  to  compen¬ 
sate  for  tracking  errors  upon  each  beam.  Initial  simulation  results  of  various  weight¬ 
ing  matrix  combinations,  revealed  an  optimum  regulator  performance  achieved  with  a 
state  weighting  factor  of  10  and  control  weighting  factor  of  0.1.  Modeled  as  zero-mean 
white  Gaussian  noise  and  zero-mean  non-white  Gaussian  noise  respectively,  measure¬ 
ment  noise  and  plant  disturbance  noise  models  were  dependent  upon  their  simulation 
correlation  time,  which  was  required  to  be  an  integer  multiple  of  the  simulation  step 
time,  due  to  simulation  restrictions.  Since  various  switching  times,  chosen  as  mul¬ 
tiples  of  the  FSM  2%  settling  time  of  0.002s,  did  not  yield  any  noticeable  positive 
or  negative  effects,  a  time  of  0.02s  was  chosen  for  simulation  purposes.  Although, 
switching  control  revealed  slight  disturbance  rejection,  the  accuracy  of  the  frequency 
response  data  was  suspect  above  a  1000  Hz.  In  addition,  based  on  desired  regulation 
time  duration,  relatively  large  to  infinitely  many  switching  times  can  occur. 

Temporal  averaging  of  the  radial  tracking  error  time  response  data  revealed  dual 
fine  tracking  control  did  not  achieve  the  residual  radial  tracking  error  limit.  Although 
dual  control  obtained  a  better  result  than  no  control,  on  average,  single  control  of 
one  beam  achieved  a  temporal  mean  of  1.6  /Trad,  which  more  than  sufficiently  meets 
the  required  tracking  error.  However,  this  performance  is  achieved  at  the  cost  of 
approximately  10.11  /trad  mean  radial  tracking  error  on  the  other  beam. 

An  implementation  that  could  exploit  these  results  is  to  only  utilize  commu¬ 
nication  data  from  each  beam  during  periods  of  applied  control.  This  still  would 
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result  in  a  single  aperture  system  with  data  received  in  short  bursts.  In  particular, 
a  time  division  multiple  access  (TDMA)  scheme,  a  common  communication  resource 
allocation  method,  could  be  the  optimal  solution. 

5. 2  Contributions 

A  dual  fine  tracking  control  methodology  has  been  developed  for  propagation  of 
two  satellite  laser  communication  uplink  beams.  The  concept  incorporated  common 
uplink  beam  propagation  disturbance  sources  with  Kalman  filtering  and  LQR  control 
to  yield  a  dual  fine  tracking  control  system. 

LQR  weighting  matrix  identification  methodology  was  developed  based  on  track¬ 
ing  error  performance.  The  concept  entailed  calculating  the  ensemble  average  and 
temporal  average  tracking  error  statistics  and  evaluating  the  associated  mean  and 
dispersion  for  each  weighting. 

Based  on  the  defined  geometry  between  the  satellite  receiver  and  UAV  transmit¬ 
ters,  the  optical  turbulence  induced  random  beam  wander  affects  on  plant  disturbance 
were  investigated  and  revealed  limited  impacts  to  tracking  error  due  to  spatial  sepa¬ 
ration  between  transmitters  within  the  receiver  field  of  view. 

5.3  Recommendations  and  Future  Work 

Based  upon  the  software  used,  simulation  of  the  control  system  required  using 
a  fixed-step  solver  methodology,  to  analyze  response  statistics.  This  required  a  fixed 
time  step  of  0.00005s  and  required  large  amounts  of  memory  to  collect  and  process 
the  data.  However,  this  simulation  was  performed  in  the  continuous-time  simulation 
environment.  In  addition,  modeled  noise  sources  required  a  correlation  time  as  an 
integer  multiple  of  the  fixed-step  time.  Although,  an  arbitrarily  chosen  correlation 
time  may  simplify  time  response  analysis,  it  may  be  inadequate  for  accurate  frequency 
response  data.  A  possible  remedy  to  this  would  be  do  simulate  the  system  in  a 
discrete-time  environment . 
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Future  work  should  include  developing  and  implementing  a  TDMA  communi¬ 
cation  scheme  within  the  simulation  environment  to  determine  impacts  to  system 
performance.  In  addition,  the  simulation  model  should  incorporate  modulation  and 
demodulation  models  to  determine  corresponding  bit  error  rate  performance. 

A  higher  level  fidelity  model  of  aircraft  attitude  dynamics  is  required  to  in¬ 
vestigate  disturbance  affects  on  transmitter  pointing  error.  In  particular,  this  could 
be  applied  to  a  two-way  communication  link  between  the  satellite  and  transmitter 
to  incorporate  the  necessary  pointing  commands  in  the  transmitter  attitude  control 
system. 

Experimental  level  evaluation  is  recommended  in  order  to  validate  the  adequacy 
of  this  approach.  This  would  require  optical  hardware  and  software  components 
assembled  in  a  controlled  vibration  environment.  Propagation  disturbance  could  be 
induced  by  either  filtering,  adaptive  optics,  or  fast  steering  mirrors. 
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Appendix  A.  Optical  Communication  System 

A.l  Data  Transmission 

During  transmission,  all  user  information  (UI)  is  converted  into  digital  informa¬ 
tion  binary  digits  (bits),  which  are  formatted  into  message  bits,  6*.  Pnlse  modulation 
converts  these  message  bits  into  a  nonreturn-to-zero  level  (NRZ-L)  pulse  code  modu¬ 
lation  (PCM)  baseband  waveform,  g(t),  as  shown  in  Figure  A.l.  For  a  given  message 
bit  rate,  Rb,  the  baseband  waveform  data  rate  and  symbol  duration  are  determined 
respectively  by 


Rg  ~  Rb  (A-l) 

T,  =  f=Tb  (A. 2) 

Kb 

The  message  sequence  is  inherently  a  random  sequence  of  bits,  in  which  it  is 
assumed  the  occurrence  of  a  1  or  0  is  equally  likely  within  the  symbol  interval,  Tbl 
and  independent  of  all  other  intervals  [12,  page  247].  Channel  symbols  1  and  0 
are  denoted  by  bipolar  pulses  of  amplitudes  +V  and  —  V  respectively,  which  occur 
with  equal  probability,  P(+V)  =  P(— V)  =  1/2.  Hence,  the  NRZ-L  PCM  baseband 
waveform,  g(t),  is  a  sample  of  a  binary  random  process  denoted  by 

[  +V  :  bi  =  1 

Git)  =  (A. 3) 

{  -V  :  b2  =  0 

for  nTb  <  t  <  {n  +  1  )Tb,  n  =  0, 1,  2,  3, . . . 

A.  1.1  NRZ-L  PCM  Baseband  Waveform  Statistics.  The  mean  and  au¬ 
tocorrelation  functions  of  a  random  process  are  the  two  moments  of  interest  for  a 
communication  system  [29,  page  23].  The  mean  is  determined  by  [10,  pages  334] 

/OO 

gfc(g)dg  (A. 4) 

•OO 
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Figure  A.l:  PCM  Baseband  Waveform:  (a)  Channel  Symbol  Sequence,  (b)  NRZ-L 
PCM  Waveform. 

where  fa{g)  is  the  probability  density  function  (pdf)  of  G(t).  Since  G(t)  is  either 
+V  or  —V  for  all  t,  its  pdf  is  determined  by  the  sum  of  delta  functions  denoted 
by  [18,  pages  58,  195] 


fG(9)  =  P(+V)5(g-V)  +  P(-V)S(g  +  V) 

=  Q)  [%  -  V)  +  S(g  +  l/)] 

Thus,  the  mean  is  calculated  as 

/OO 

gfc(g)dg 

-oo 

=  (0/  g[d(g -V)  +  5(g  +  V)\dg 

=  Q)m+(-^)] 

=  0 


(A.5) 


(A.6) 


The  autocorrelation  of  G(t)  is  determined  by  [10,  pages  335] 


RG(t  1,^2)  —  E  [G(ti)G(t2)\ — 


gig-2  [fGi,G2(gi,  gz)]  dgxdg2  (A. 7) 


—00 


—00 


where  fcxxnifh ,  g2)  is  the  joint  pdf  of  G(t\)  and  G(t2),  and  ti  and  t2  are  arbitrary 
time  constants.  Similar  to  the  pdf  derivation,  the  joint  pdf  is  determined  by  the  sum 
of  delta  functions  denoted  respectively  by  [18,  pages  72,  79,  195] 


fGl,G2(gi,g2)  —  fGl,G2(gi,g2\tg)fTb(tg) 

=  [P{91  =  +V,  g2  =  +V)5(gi  -V,g2-  V) 

+  P(9l  =  +V,  g2  =  -V)5(gi  -V,g2  +  V ) 

+  P(9l  =  -V,  g2  =  +V)6(gi  +  V,g2-  V) 

+  P(9l  =  -V,g2  =  -V)5(gi  +  V,g2  +  V)]  fn(tg )  (A.8) 


where  tg  is  a  uniform  random  variable  in  the  interval  [0,77],  which  accounts  for  the 
unsynchronized  realizations  of  G(t),  and  fTb(tg )  is  the  associated  pdf  defined  as  [12, 
page  248] 


fn(tg) 


1/ 77  :  0  <  tg  <  Tb 
0  :  elsewhere 


(A.9) 


If  1 12  —  ti|  =  |t|  and  |r|  +  tg  >  Tb,  G(ti)  and  G{t2)  occur  in  different  pulse  intervals 
and  are  independent.  Thus,  the  joint  conditional  pdf  is  determined  by  [18,  page  195] 


/Gl,G2(<7l,  92\tg) 


P(gi  =  +V)P(g2  =  +V)5(gi  -V,g2-  V ) 

P(9i  =  +V)P(g2  =  -V)S(g  1  -V,g2  +  V ) 

P{9i  =  ~V)P(g2  =  +V)S{91  +  V,g2-  V) 

P(gi  =  ~V)P(g2  =  -V)S{g  1  +  V,g2  +  V ) 

Q  S(gi  -V,g2-V)+  ^  8(gi  ~V,g2  +  V ) 

(^4)  92  —  V)  +  5(gi  +  V,  g2  +  V)  (A. 10) 
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which  implies  an  autocorrelation  function  defined  as 


1^2)  — 


+ 


poo  poo  poo 

/  /  .9 1 .92  fc  1  ,G2  (.9 1 ,  .92 1  tg)fn  (tg )  dg  i  dg2dtg 

—  OO  J  —OO  J —OO 

Tb-\r\  !■  00  /-oo 

T  J  J  I  I  gig2fGi,G2(gi,g2\tg)dgidg2dtg 

1\  ( Th  -  Irl 


— 00  */  — OO 

/»oo  /»oo 


9i92<%i  -  V,  g2  -  V)dg1dg2 


'  — OO  «/  — OO 


r»oo  poo 


'  — OO  «/  — OO 
poo  poo 


'  — OO  «/  — OO 

poo  poo 


S'132%1  -  v,  g2  +  V)dg1dg2 
gmSigi  +  v,g2-  V)dg1dg2 
gig2d{gi  +  V,g2  +  V)dgidg2 


'  — OO  «/  — OO 


j  -  hL  j  [(u)(v)  +  (v)(-v)  +  (-c)(c)  +  (-v)(-v)] 

n  p-M)  [2V2-2V2]=0  (A. 11) 

TbJ 


However,  if  |r|  +tg  <  Tf,,  G(ti)  and  G{t2)  occur  within  the  same  pulse  interval.  Thus, 
the  joint  conditional  pdf  is  defined  by  [18,  page  195] 


/c?l,G2(<7l,  92\tg) 


P(g2  =  +v\9l  =  +V)P{91  =  +V)6(g1  -V,g2-  V ) 

P(g2  =  ~V\9l  =  V)P{91  =  V)5(g i  -  Vj92  +  V ) 

P(g2  =  +v\9l  =  -V)P{g  1  =  -V)S(9l  +  V,g2-  V ) 

P(g2  =  -v\9l  =  -V)P(9l  =  -V)%,  +  V,g2  +  V ) 

(1)  Q)  <%i  -V,g2-V)  +  (0)  Q)  <%i  -  V,g2  +  V ) 

^2^]  —  ^ ^2  —  V)  +  +  V, g2  +  H)]  (A. 12) 


<%i  +  V,g2  -  V)  +  (1)  (  -  )  +  V,g2  +  V ) 
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which  implies  an  autocorrelation  function  defined  as 


RG(t  i,t2)  — 


poo  poo  poo 

/  /  9i92fGi,G2(9u92\tg)fD,{tg)dgidg2dtg 

—  OO  J  —  OO  J  —  OO 

T6-|r|  p  oo  po o 

T  i  i  II  gi92fGi,G2(9i,g2\tg)dg1dg2dtg 


i  o 


r 

1_  n 


—oo  J  —oo 

roo  poo 


'  —oo  J  —oo 


9i92$(gi  ~V,g2-  V)dg1dg2 
9i92d(gi  +  V,  g2  +  V)dg]dg2 
[(C)(V)  +  (-C)(-V)] 

3 11  -  !0  y ^ 


poo  poo 

—00  J  —00 

1 


=  VZ,1~Tk 


(A.13) 


Therefore,  the  autocorrelation  function  is 


RG{t  i,t2)  —  Rg(t)  — 


\r\/Tb)  : 

r 

T 

^  J-b 


(A.14) 


In  particular,  based  on  the  mean,  mG(t)  =  0,  and  the  autocorrelation  function, 
Reid  1,^2)  =  Rg(t),  the  autocovariance  of  G(t)  is  dehned  as  [10,  page  335] 


Cc(ti,t2)  —  E  [{G(fi)  —  mc(ti)}{G(t2)  — 

=  E  [{Gih)  -  0}{G(f2)  -  0}]  =  E  [G(h )G(t2)} 

=  RG(ti,t2)  =  Rg(t)  (A. 15) 


Since  the  mean  is  constant  and  the  autocovariance  is  a  function  of  r  =  t2  —  ti,  G(t) 
is  a  wide-sense  stationary  process  (WSS)  [10,  page  359]. 

A.  1.2  NRZ-L  PCM  Baseband  Waveform  Power  Spectral  Density.  The  base¬ 
band  waveform  power  spectral  density  (PSD),  determined  from  the  Fourier  transform 
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of  Rg(t),  is  derived  as  [10,  page  404] 


Saif) 


J  —  oo 

sin(7 xfTb) 


(v2n 


KfTb 


(V2Tb)smc2(fTb) 


Rc(j)  exp(— j27r/r)dr 
1  2 


(A. 16) 


Based  on  the  Chapter  III  defined  data  rate,  Rb  =  100  Kbps,  NRZ-L  PCM  wave¬ 
form  PSD  is  shown  in  Figure  A. 2,  scaled  for  unit  amplitude.  The  average  power  is 
determined  by  [29,  page  25] 

flc(o)  =  V"2  fl  -  4)  =  V2  (A.17) 


Figure  A. 2:  NRZ-L  PCM  Baseband  Waveform  PSD 
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A. 2  Optical  Receiver 


A  conceptual  00K  optical  receiver  is  illustrated  in  Figure  A. 3.  After  optical 
signal  processing  as  defined  in  Appendix  D,  an  avalanche  photodiode  (APD)  detects 
and  converts  the  the  received  optical  power  into  an  amplified  current  signal,  ir(t).  A 
matched  filter  maximizes  the  signal-to-noise  ratio  (SNR)  prior  to  sampling  a  current 
test  statistic,  z(Tb),  for  the  decision  circuit.  Due  to  noise,  N0,  distortion  effects 
inherent  within  the  received  signal,  the  decision  circuit  observes  the  received  signal 
and  determines  whether  si  is  sent,  denoted  by  hypothesis  Hi,  or  s 2  is  sent,  denoted 
by  hypothesis  H2  [6,  page  534],  After  symbol  decision,  based  on  the  threshold  7 D,  the 
estimated  received  symbol  is  assigned  to  its  corresponding  estimated  message  bits,  6j, 
which  in  turn  are  formatted  back  into  user  information. 

A. 2.1  Received  Current.  The  received  current  signals,  based  on  the  unipolar 
transmitted  signals  within  the  interval  t  E  [0,  Tb\  defined  in  Appendix  D,  are  defined 
respectively  as 


h(i) 

Oc  ~t~  in(t) 

=  RphPr  +  in{t)  [A] 

(A. 18) 

*2  (t) 

—  0  +  in(t) 

=  in(t)  [A] 

(A. 19) 

Figure  A. 3:  OOK  Optical  Receiver 
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where  idc  is  the  dc  current,  R.ph  is  the  photodetector  responsivity,  Pr  is  the  received 
optical  power,  and  i„  is  the  current  noise. 


A. 2.2  Matched  Filter.  In  order  to  maximize  the  SNR  for  sampling  and 
detection,  the  received  current  signal  is  filtered  by  a  matched  filter.  The  maximum 
output  SNR  is  defined  by  the  ratio  of  the  input  signal  energy  and  noise  power  spectral 
density,  defined  as  [29,  page  127] 


SNRt 


{Ejt^t)}  -  E[i2(t)}}2  _  2 Ed 

N0 


(A. 20) 


where  Ed  is  the  input  energy  difference  signal,  a20  is  the  output  noise  variance,  and 
N0/2  is  the  two-sided  power  spectral  density  of  the  input  noise.  N0  is  associated 
with  thermal  noise  and  is  the  dominant  photodetector  noise  source,  as  determined  in 
Appendix  D.  Thermal  noise  is  characterized  as  a  zero-mean  additive  white  Gaussian 
noise  (AWGN)  and  is  used  to  model  communication  system  noise  [29,  page  33].  Thus, 
based  on  the  noise  statistics,  the  matched  filter  output  SNR  is  determined  as 


SNRt 


i dc  _  -Ed 
O2o  ~  N0 


(A.21) 


The  ideal  energy  difference,  based  on  noiseless  unipolar  input  current  signals  within 
the  time  interval  [0,  T&],  assuming  a  unit  Ohm  resistance,  is  determined  as  [29,  page  128] 


Ed 


A  EfaW-Efoit)]  }2dt 
Jo 

(  Rph  Pr  )  2  Tb  [J] 


(- RphPrfdt 


(A. 22) 


A. 2. 3  Symbol  Detection.  Assuming  symbol  synchronization  is  maintained 
for  optimum  demodulation,  the  symbol  detection  goal  is  to  determine  whether  the 
sampled  current  test  statistic  from  the  matched  filter  output,  z(Tb),  satisfies  either 
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hypothesis  Hi  or  II2-  determined  respectively  as 


Hi  :  z(Tb)  —  idc  +  in(Tb)  (A. 23) 

H2  :  z{Tb)  =  in{Tb )  (A. 24) 

Since  z(Tb)  is  a  function  of  in(Tb),  a  zero-mean  Gaussian  random  variable,  it  is  also 
subject  to  Gaussian  statistics.  Thus,  the  corresponding  conditional  pdfs  for  either 
signal  being  sent  are  denoted  respectively  by  [10,  page  113] 

^  (A-25) 

p(*'i!,2)  =  (^)exp{-^}  (a-26) 

Assuming  the  transmitted  unipolar  signal  waveforms  have  an  equally  likely  a  priori 
probability,  P(s i)  =  P(s2 )  =  1/2,  the  conditional  pdfs  are  symmetric,  as  shown  in 
Figure  A. 4.  Thus,  the  decision  criterion,  to  estimate  which  symbol  was  sent,  is  based 
on  the  maximum  likelihood  criterion  as  [29,  page  110] 

z(Tb)  > 
z(Tb)  < 


7 o  =>■  symbol  1  (A. 27) 

7 o  symbol  0  (A. 28) 


where  70  is  the  optimum  decision  threshold.  The  optimum  threshold  is  determined 
from  the  minimum  error  criterion,  which  is  a  function  of  the  means  of  each  conditional 
pdf  denoted  by  [29,  page  182] 


7  o  = 


^ dc 

T 


(A. 29) 


A. 3  Probability  of  Bit  Error 

Based  on  the  maximum  likelihood  criterion,  there  are  two  ways  a  symbol  detec¬ 
tion  and  estimation  error  can  occur  for  a  binary  system:  false  detection  and  missed 


75 


detection.  The  probability  of  error  is  calculated  by  the  summation  of  the  probabilities 
of  these  two  occurrences  [29,  page  1044],  False  detection,  also  referred  to  as  false 
alarm,  occurs  when  the  sampled  voltage  exceeds  the  voltage  threshold  when  s2  is  sent. 
Conversely,  missed  detection,  referred  to  as  probability  of  fade,  occurs  when  si  is  sent 
but  the  sampled  voltage  does  not  exceed  the  voltage  threshold.  Thus,  the  probability 
of  symbol  error,  PE,  is  defined  as  [29,  page  1044] 

Pe  =  P(Si)PEade  +  P(S'2)PFalse  (A. 30) 

Since  OOK  modulation  is  an  orthogonal  binary  signaling  set,  the  probability  of  bit 
error  is  equivalent  to  the  probability  of  symbol  error,  denoted  by  [29,  page  551] 

Pb  =  Pe  =  P(Si)Ppade  +  P{s2)PFalse  (A. 31) 

A. 3.1  Gaussian  Communication  Channel.  Assuming  the  communication 
channel,  defined  in  Appendix  E,  behaves  as  an  AWGN  channel  the  corresponding 
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false  alarm  and  fade  probabilities  are  determined  respectively  by  [29,  page  1044] 


P , 


False 


p(z\s2)  di 


J  To 

/To 

-OO 


(A. 32) 
(A. 33) 


Thus,  the  probability  of  bit  error  is  determined  as 


Pb  —  P(Si)Ppade  +  P(yS2)PFalse 


=  P(s  l) 


f-To 


p  (z|si)  dz 


+  P{s2) 


p  (z|s2)  dz 


LJ7o 


(A. 34) 


Since  both  symbols  are  equally  likely  to  have  been  sent  and  symmetry  of  the  condi¬ 
tional  pdfs,  the  probability  of  bit  error  simplifies  to  [29,  page  121] 


Let 


Pb  = 


no 


no 

/•OO 


'7° 


p  (^|si)  dz 
p  (z|s2)  dz 

exp 


o0\/2n 


+  l^> 


p  {z\s2)  dz 


Lo'To 


2a„2 


dz 


(A. 35) 


Thus 


z 


u  = 


O0 

7o  <  z  <  oo 
7  o 


X  = 


On 


Pb  = 


oQdu  =  dz 


7o  ,  , 

—  <  u  <  oo 

Oo 

^dc 


2  On 


exp 


u 


du  =  Q  (x) 


(A. 36) 


where  Q(x)  is  the  complementary  error  function  [29,  page  1044],  Figure  A. 5  illustrates 
the  relationship  between  Q(x)  and  an  arbitrary  x  G  [0,6],  commonly  referred  to  as 
the  error  curve.  Substituting  in  the  predetermined  matched  filter  output  SNR,  the 
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Figure  A. 5:  Complementary  Error  Function  Curve 


arbitrary  Q-function  parameter  becomes 


(A. 37) 


where  Eb  is  the  average  energy  per  bit,  defined  as  Eb  =  Ed/ 2  for  orthogonal  sig¬ 
nals  [29,  page  132],  Thus,  the  probability  of  bit  error  for  OOK  modulation,  shown  in 
Figure  A. 6,  is  determined  as 


Pb 


(A. 38) 


A. 3.2  Fading  Communication  Channel:  Non- Ideal  Tracking.  Assuming  the 
incident  beam  irradiance,  defined  in  Appendix  D  as  I(pz,zp),  fluctuates  at  a  slower 
rate  relative  to  the  sample  time  interval,  Tb,  the  average  received  optical  power,  Pr,  is 
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Figure  A. 6:  OOK  Modulation  Theoretical  BER  Curve 

constant  [3,  page  268].  However,  tracking  errors,  quantified  by  the  tracking  loss  Lr, 
attenuates  the  received  power.  Thus,  the  random  tracking  error  fluctuates  incident 
beam  irradiance  about  the  propagation  axis,  zp,  creating  a  fading  effect  within  the 
communication  channel. 

As  defined  previously  in  Chapter  If,  Pr  is  determined  from  the  Friis  transmission 
equation  by  [15,  pages  94-105] 

Pr  Pt{iltGt.LtLao)  {LsLa)  (LrGrr]r) 

where  Pt  is  the  transmitted  power,  r]t  is  transmitter  efficiency,  Gt  is  the  transmitter 
antenna  gain,  Lt  is  the  pointing  loss,  Lao  is  aerodynamic-optical  loss,  Ls  is  the  path 
loss,  La  is  atmospheric  loss,  Lr  is  the  tracking  loss,  Gr  is  the  receiver  antenna  gain, 
and  7]r  is  the  receiver  efficiency  loss.  Let 

Oir  —  Pti,VtGtLtPao)(LsLa)(r)r)  (A. 39) 
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Assuming  ar  and  Gr  are  known  a  priori,  Pr  can  be  defined  as  a  function  of  the  radial 
receiver  tracking  error,  pri  denoted  by 

Pr(pr )  =  ctrGrLr  =  arGr  exp  (— Grp^)  (A. 40) 

As  defined  in  Chapter  II,  the  receiver  tracking  error  is  subject  to  Rayleigh  probability 
distribution,  denoted  by 

h(pr)  =  (A)  exp  (-^)  ■  Pr>0  (A.41) 

Since  Pr  is  an  exponential  function  of  a  Rayleigh  random  variable,  the  probability  of  an 
event  involving  Pr,  defined  on  the  interval  (0,  arGr],  is  equivalent  to  the  corresponding 
event  involving  pr,  defined  on  the  interval  [0,oo)  [10,  page  120].  Thus 

poo 

P[0  <  Pr  <  arGr]  =  P@=  fi{(pr)dpr  (A. 42) 

Jo 

Hence,  the  probability  of  fade  and  false  detection,  assuming  pr  and  in(t)  are  indepen¬ 
dent  random  variables,  are  determined  respectively  by  [10,  page  209] 


P 


Fade 


p{pr>  on  [z(Tb)  <  7o|si]} 

=  P\pr  >  0 \P[z(Tb)  <  7o|si] 

fR(Pr)dpr>  \  /  p(z\si)di 


[p(^|si)  fR{pr)}  dz  dpr 


(A. 43) 


P> 


False 


P{pr>  on  [z(Tb)  >  7o|s2]} 

=  P[pr  >  0 }P[z(Tb)  >  7o|s2] 

/*oo  ^  r  poo 

fR(pr)dpr>  \  /  p(z\s2)dz 


0 

/*oo  /*oo 

>0  J~fo 


'lo 


[p{z\s2)  fR{pr)\dz  dpr 


(A. 44) 
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The  probability  of  bit  error  is  determined  as 


Pb  ~  P(sl)PFade  +  P{S‘2)P False 

-y  (  roo  r'jo  roc  roo 

=  /  [p(z\s1)fR(pr)]dzdpr+  /  \p  (z\s2)  fR(pr)\  dzdpr 

Z  IJo  J- oo  JO  J  7o 

r°°  (  rio  roo  \ 

=  2  J  |  J  P(z\si)dz  +  J  p(z\s2)dzj  fR{pr)dpr 

^  roc  r  roo  'j 

=  2  J  [2j  P(Z\S2)dzj  fR(pr)dpr 

roo  (  roc 


p(z\s2)dz  \  fR(pr)dpr 

(  Period  )  f  R  ( Pr  )  dpr 


(A. 45) 


where  the  modulation  probability  of  bit  error  is  determined  as 


Pmod  —  Q  (  \l  I  —  Q 


Ed  \  „  I  li2dcTb 


2Nn 


=  Q 


2Nn 


=  Q 


=  Q 


[RphPr(Pr)]' 


( Rph ) 


2NaR,h 
21  1/2 


=  Q 


( Rph? 


1/2 


L(/V0)(2i?fe)J 
Pr(Pr )  }  =  Q{KrPr{pr)} 


Pr  ( pr ) 


=  Q  [7iravG>exp  (— Grp^)] 


(A. 46) 


Thus,  the  probability  of  bit  error  for  a  fading  communication  channel,  due  to  tracking 
loss,  is  defined  by 

[KrarGrexp  (-Grp2r)\  fR(pr)dpr  (A.47) 
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A. 4  Link  Margin 

Comparison  of  the  detected  received  signal  power  with  receiver  sensitivity,  Rxs 
determines  the  system  link  margin,  denoted  in  dB  by 


Link  Margin  =  Pr  [dB]  —  Rxs  [dB] 


(A. 48) 


A. 4-1  Received  Optical  Power.  The  parameters  used  to  determine  the  de¬ 
tected  received  signal  power  are  determined  from  the  Friis  transmission  equation,  and 
the  avalanche  photodiode  (APD)  gain  as 


Pr  —  ( PtVtGtLtLao)(LsLa)(LrGrr]r)Gd 


(A. 49) 


where  Gd  is  the  APD  gain  defined  in  Appendix  D.  Based  on  the  Chapter  III  defined 
operational  parameters,  the  theoretical  detected  received  optical  power  for  each  beam 
is  shown  in  Figure  A. 7.  Although,  a  larger  wavelength  is  less  susceptible  to  atmo¬ 
spheric  loss,  at  higher  elevation  angles  the  shorter  wavelength  achieves  slightly  better 
power  performance. 

A. 4-2  Receiver  Sensitivity.  Since  the  noise  statistics  are  Gaussian,  in  addi¬ 
tion  to  the  optimum  decision  threshold,  the  complementary  error  function  provides 
an  estimate  of  the  the  receiver  performance  [1,  page  42],  Thus,  the  receiver  sensitiv¬ 
ity  is  determined  as  the  minimum  received  power  required  to  maintain  the  desired 
system  BER,  for  a  given  message  data  rate  [1,  page  43].  For  the  purposes  of  this 
study,  the  minimum  received  power  which  produces  a  system  BER  of  10-6,  common 
for  free-space  communications,  will  be  the  receiver  sensitivity. 
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Figure  A. 7:  Theoretical  Detected  Received  Optical  Signal  Power 

Based  on  the  complementary  error  function  curve,  as  shown  in  Figure  A. 5,  the 
receiver  sensitivity  is  determined  as  [1,  page  246] 

Q(x)  =  1(T6  =>•  x  «  4.753 

Rx,  =  P1)  (A)(4) 1/2  TO  (A.50) 

where  i2n  is  the  total  equivalent  noise  power  at  the  APD.  Thus,  for  each  transmitted 
beam,  the  receiver  sensitivity  is  shown  in  Figure  A. 8. 

A. 4-3  Link  Margin.  Based  on  the  detected  received  optical  power  and 
receiver  sensitivity,  for  a  transmitted  symbol  1,  the  system  link  margin  is  shown  in 
Figure  A. 9.  Thus,  within  the  receiver  held  of  view  (FOV),  the  longer  wavelength 
achieves  approximately  0.1  dB  greater  link  margin  performance.  In  either  case,  both 
beams  achieve  a  positive  margin  above  1.5  dB,  assuming  no  tracking  loss  effects,  which 
will  satisfy  the  required  system  BER  performance  [29,  page  265]. 
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Figure  A. 8:  Receiver  Sensitivity 
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Figure  A. 9:  Link  Margin 
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A .  5  Summary 

This  appendix  presented  the  fundamental  communication  theory  aspects,  prob¬ 
ability  of  bit  error  and  link  margin,  of  the  laser  communication  architecture.  Proba¬ 
bility  of  bit  error,  commonly  referred  to  as  bit  error  rate  (BER),  provides  a  measure  of 
communication  system  performance,  which  is  dependent  upon  received  signal  power 
and  radial  tracking  error  statistics.  In  addition,  based  on  the  modulation  scheme, 
the  system  BER  is  used  to  determine  the  receiver  sensitivity  and  link  margin,  which 
determines  if  adequate  signal  power  is  available  for  communications. 
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Appendix  B.  Satellite  Dynamics 


B.l  Satellite  Position 

The  relative  satellite  orbital  motion  about  Earth  can  be  described  by  the  two 
body  problem  equation  of  motion  [34,  page  50] 


d2r , 
dt2 


(B.l) 


This  relation  satisfies  both  conservation  of  energy  and  angular  momentum  per  unit 
mass  denoted  by  [34,  pages  50-52] 


H  =  rs  x  vs 


(B-2) 

(B.3) 


where  rs  is  the  satellite  position  vector,  vs  is  the  satellite  velocity  vector,  E  is  the 
energy  per  unit  mass,  H  is  the  angular  momentum  vector,  and  y=  3.986  x  105  km3/s2 
is  the  Earth’s  gravitational  constant. 

B.1.1  Orbit  Location.  Orbit  orientation  relative  to  a  geocentric  inertial 
(GCI)  reference  frame,  defined  by  unit  vectors  i,  j,  and  k  as  shown  in  Figure  B.l, 
and  satellite  motion  within  the  orbital  plane  are  determined  by  the  following  classical 
orbital  elements  [34,  page  60]: 


e  =  eccentricity 
a  =  semimajor  axis 
i  =  inclination 
T0  =  time  since  perigee  passage 

=  right  ascension  of  ascending  node 
u  =  argument  of  perigee 


satellite 


k 


Figure  B.l:  GCI  Frame  Satellite  Orbit 
Orbit  energy  and  angular  momentum  magnitude  are  determined  by  [34,  page  56] 


E 

H 


2  a 

'//Ja  i  I  -e2) 


(B.4) 

(B.5) 


A  circular  orbit  assumption  implies  e  =  0  and  cu  can  not  be  determined.  How¬ 
ever,  an  alternate  orbital  element,  known  as  the  argument  of  latitude,  ua  ,  is  used 
to  determine  relative  orbital  position  from  the  line  of  nodes,  n,  also  referred  to  as 
the  nodal  unit  vector  [16,  page  134],  Within  the  equatorial  plane,  O  determines  the 
angular  position  of  the  nodal  vector  from  the  vernal  equinox,  i,  as  shown  in  Fig¬ 
ure  B.l.  Using  conic  section  geometry,  based  on  the  GCI  frame,  the  nodal  unit  vector 
is  [34,  page  62] 


h  = 


kx  H 
kxH 


cos(fl)i  +  sin(fl)j 


(B.6) 
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Figure  B.2:  ECF  Frame  Satellite  Orbit 

Accounting  for  Greenwich  sidereal  time,  0g(^o),  angle  between  vernal  equinox  and 
Greenwich  meridian  at  epoch,  satellite  longitude  location,  Xaz,  and  the  orbital  angular 
velocity,  ua,  the  argument  of  latitude  over  time  is  determined  by  [28,  page  43] 

uQ(t)  =  ua(t  - 10)  +  0G(to)  +  xaz-n  (B.7) 

B.1.2  Satellite  Location.  In  terms  of  the  Earth  centered  fixed  (ECF)  ref¬ 
erence  frame,  defined  by  unit  vectors  x,  y,  and  z  illustrated  in  Figure  B.2,  relative 
position  from  Greenwich  meridian,  x  =  9G(t0),  is  determined  by  [16,  page  100] 

rxyz  =  a  [cos(Aa2)  cos ((pei)x  +  sin(Aa2)  cos ((f)ei)y  +  sin {<j>d)z\  (B.8) 

where  Aa2  is  the  azimuth,  also  referred  to  as  longitude  location,  and  (j)ei  is  the  eleva¬ 
tion,  also  referred  to  as  latitude  location.  The  satellite  GCI  reference  frame  position 


vector,  fs,  is  related  to  the  ECF  reference  frame  position  vector,  fxyz,  by  the  following 
rotational  transformation: 

rxyz  =  REGrs  (B.9) 


where 

Reg  =  RsMt^R^R^n)  =  ( Rge)T  (B.10) 

The  corresponding  direction  cosine  matrices,  for  positive  counter  clockwise  rotations, 
are  determined  by 


Ri(a) 


R2{pi) 


R3(a) 


1  1  0  0  ^ 
0  cos(a)  sin(a) 

^  0  —  sin(a)  cos(ct)  y 

^  cos(a)  0  sin(a)  ^ 
0  1  0 
l  —  sin(a)  0  cos(a)  y 

^  cos(a)  sin(a)  0  ^ 
—  sin(a)  cos(a)  0 

0  0  V 


Thus,  GCI  position  and  velocity  vectors  are  determined  as 


rs  =  RGErxyz  (B.ll) 

dv 

vs  =  tl/'  -d,  x  rs  (B.12) 


where  CJS  is  the  satellite  angular  velocity  vector. 


B.2  Satellite  Attitude  Dynamics 

As  illustrated  in  Figure  B.3,  satellite  attitude  within  orbit  is  defined  about  either 
the  orbit  reference  frame  or  satellite  body  reference  frame.  Vectors  xQ ,  ya,  zQ  represent 
the  orbit  reference  axis  frame:  za  is  the  negative  unit  vector  of  fs ,  xQ  is  the  unit  vector 


Figure  B.3:  Satellite  Attitude 

of  va,  and  ya  is  normal  to  the  orbit  plane.  Unit  vectors  xs,  ys,  zs  define  the  satellite 
body  axis  frame.  Attitude  motion  of  the  satellite  platform  about  its  center  of  mass 
(cm)  is  defined  by  Euler’s  rotational  equation  of  motion  as  [28,  page  95] 


hs  +  usxhs 


T, 


(B.13) 


where  hs  is  the  satellite  body  angular  momentum  vector  and  Ts  is  the  total  torque 
vector  acting  on  the  satellite  body.  Assuming  an  axisymmetric  rigid  body  about  the 
zs  axis  and  alignment  with  the  principal  axes  of  inertia,  the  satellite  body  angular 
momentum  vector  is  determined  by 


h°  I  c  ecu 
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UJy 
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IxLOy 
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Uz  / 

IzUz 

(B.14) 
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where  I*  is  the  satellite  moment  of  inertia  matrix.  Tims,  attitude  equations  of  motion 
can  be  represented  by  the  following  nonlinear  equations: 


UJy 

Uz 


(A  ~  h) 

lx 

( lx  -  Iz) 

lx 


Tx 

+  — 

T 
.  y 

^X^Z  “1“  —j— 
*x 


Tz 

Iz 


(B.15) 


B.2.1  Attitude  Kinematics.  Applying  attitude  kinematics,  ujs  is  dehned 
as  [28,  page  102], 

ujs  a) so  T  u)0is  (B.16) 

where  ujso  is  the  satellite  body  relative  to  orbit  frame  angular  velocity  vector  and  a )ois 
is  the  angular  velocity  vector  of  the  orbit  reference  frame  relative  to  the  GCI  frame  in 
satellite  body  coordinates.  Satellite  attitude  with  respect  to  the  body  axis  reference 
frame,  as  shown  in  Figure  B.3,  is  dehned  by  Euler  angles:  roll  angle,  0S,  about  the 
xs  axis,  pitch  angle,  9S,  about  the  ys  axis,  and  yaw  angle,  if;s,  about  the  zs  axis. 
Assuming  a  yaw,  pitch,  and  roll  order  of  rotational  transformation,  ujso  is  determined 
as  [28,  page  103] 


0 1 

(°) 

UJso  =  Ri{(j)s)R2(0s)R3{i>s) 

0 

+  RMR2(9s) 

+  Ri(cj)s) 

0 

(B-17) 

U*  J 

0  J 

0 ) 
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For  small  Euler  angles,  where  cos(a)  ~  1  and  sin(a)  ss  a,  the  rotational  transforma¬ 
tion  matrices  are  approximated  respectively  by 


R1{<j>s)R2{6s)R^s) 


RMR2(os) 


Ri{(ps) 


1  i  Vb  —o8  ^ 

0 s  1  0s 

Y  0S  0S  1  J 

1  1  0  9S^ 

0  1  0S 

y  os  0s  i  y 

(  1  0  0  ^ 

0  1  0s 

y  0  0s  1  J 


(B.18) 


(B.19) 


(B.20) 


Thus, 


\ 


W*/ 


(B.21) 


To  determine  cn0j,  angular  rate  of  change  of  each  orbit  reference  frame  axis,  as  shown 
in  Figure  B.4,  respectively  is  determined  by, 


dx0 

dt 


dy0 

dt 


dzQ 

dt 


=  VoU  3 


=  -Z0UJ  2 


dxQ  dy 
dy  dt 
dxQ  dj3 
d/3  dt 
dyQ  da 

da  dt 
dya  dy 

dy  dt 
dz0  df3 

d0  dt  =  XoUJ2 
dzQ  da 


da  dt  °  1 


XqUJs 


da  dt 


=  -y0v  i 
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Figure  B.4:  Orbit  Reference  Frame  Rate  of  Change 


After  performing  a  dot  product  between  each  rate  of  change  and  associated  orbit 
reference  frame  axis,  the  respective  angular  velocities  are  [28,  page  106] 


dx0 

dt 


■Vo 


dy0  A 


dt  'Z‘ 
dza  , 

~dt'Xl 


dza  „ 


dt  'Vo 
dXn  , 


dt 


\y0U2,\  \y0\  cos(0)  =  u;3 
\z0UJ\  \  \z0\  cos(0)  =  U>1 
\x0UJ2\  \xa\  cos(o)  =  u2 


(B.22) 


Substituting  in  the  predetermined  GCI  coordinates  for  each  respective  orbit  reference 
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frame  axis 


w  i  = 


dza  „  d  f  —rs 
dt  dt  1  IrJ 


(vs)  ■  (rs  xvs)  =  0 


dfx 


rsxvs\  J  rs  [r svs  sin(90)]  dt 


r2vs 


UJ2  — 


w  3  — 


rsv2s 


dzQ  „  (  1  \  dry  /  ry  \  -1  _s  -ry 

—  -a;0= - — .  —  = - (ry  •  v8)  = - cos(O)  =  -u0 

dt  V  rsJ  dt  \|ry|/  r'svs  rs 

dx0  „  _  /  f s  x  vs  \  _ _ 1 _ d2rs  _ 

dt  \VSJ  dt  \|ryxfy|/  vs  [r3vs sin(90)]  dt 2  s  s 

y#)  ^  x  =  (^)  ^  x  v*)  =  0 


Thus, 


UJl 

r  0  ' 

^ oi 

U)2 

= 

\  W3  / 

l  0  J 

(B.23) 


In  terms  of  the  satellite  body  axis  reference  frame,  for  small  Euler  angles, 


^ ois  Rl{4>s)R2iPs)Rz{}lJs)^oi 


^  ~UJ0ljjs  ^ 

tl^o 

y  J 

Thus,  the  satellite  angular  velocity  vector  is 


^  4>s  -  ^ 

LUS  = 

Uy 

= 

0  s  Ld 0 

\U*  7 

^  ^Ps  H“  ^O^ps  J 

(B.24) 


(B.25) 
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Substituting  into  the  nonlinear  attitude  dynamic  equations,  accounting  for  small  Euler 
angles  where  aa  ~  0  and  a/3  ~  0,  satellite  attitude  motion  is  defined  by  linearized 
rotational  equations  of  motion  as, 


fs 

Os 

4>s 


UJ0 


i  + 


(4  -  4) 


Ty 

4 

Tz 

—UJq^s  +  ~j~ 


Ips  +  ul 


(4  -  Q 

4 


(B.26) 


B.2.2  Satellite  Field  of  View.  The  satellite  FOV  geometry,  as  shown  in 
Figure  B.5,  is  determined  by  [16,  page  108] 


Pe 

=  sin  1  ( RE/rs ) 

(B.27) 

Po 

=  90°  —  pe 

(B.28) 

De 

=  Re  tan(p0) 

(B.29) 

where  pe  is  the  radial  distance  of  Earth  visible  from  orbit,  Re  is  the  Earth’s  radius,  pQ 
is  the  radial  distance  of  Earth  from  the  Earth’s  cm,  and  De  is  the  Earth  horizon  dis¬ 
tance,  also  referred  to  as  slant  range.  Based  on  the  satellite’s  operational  constraints, 


Figure  B.5:  Satellite  Field  of  View  Geometry 
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defined  in  Chapter  III,  these  parameters  are  determined  as, 


Pe 

=  sin"1  (6, 378/42, 164)  =  8.7°  [0.152  rad] 

(B.30) 

Po 

=  90°  -  8.7°  =  81.3°  [1.419  rad] 

(B.31) 

De 

=  (6,  378)  tan(81.3)  =  41,  679  km 

(B.32) 

Due  to  the  aircraft’s  operational  constraint  on  minimum  elevation  angle,  defined  in 
Chapter  III  as  20°,  the  actual  satellite  FOV  is  determined  by, 

FOVs  =  2rja  (B.33) 


where  rja  is  the  nadir  angle,  radial  distance  from  SSP  to  aircraft  visible  from  orbit.  The 
angular  relationships  between  aircraft  and  satellite  longitude  and  latitude  locations, 
(A az,(pei )  and  (A a,  (pa)  respectively,  are  defined  as  [16,  page  111] 


Pas 

^ as 


=  cos  1  [sin 


(0Q)  -  cos(pas)  sin (<f)el)/  (sin(pas) 
nfmlA 


=  cos 


"sin  (rfg) 
,sin(pe) 


COb^e^Jj 


(B.34) 


where  pas  is  the  radial  distance  from  SSP  to  aircraft  along  the  Earth’s  surface,  defined 
as  the  Earth  central  angle,  Aas  is  the  azimuth  direction  from  SSP  to  aircraft,  and  esa 
elevation  angle  from  aircraft  to  satellite.  The  satellite  azimuth  direction  from  the 
aircraft,  measured  from  true  north,  is  determined  by 


0°  <  Aas  <  180°  =>  Asa  =  180°  +  \as  (B.35) 

180°  <  Aas  <  360°  ^  Asa  =  Aas  -  180°  (B.36) 


96 


For  a  given  aircraft  position,  ra,  the  nadir  angle  and  maximum  range  are  determined 
by  [16,  page  111] 


Pa 


D 


a 


90  Pas  &sa 

sin"1  [sin(pe)  cos(esa)] 
'sin  (pas)' 

ra  .  ,  , 

Sin  (77a) 


(B.37) 

(B.38) 

(B.39) 


where  Da  is  the  maximum  range  between  the  satellite  and  aircraft.  Thus,  for  the 
predefined  operational  constraints,  the  maximum  nadir  angle,  satellite  FOV,  Earth 
coverage,  and  maximum  range  are  determined  as 


?/a  =  sin  1  [sin(8. 7)  cos(20)]  =  8.172°  [0.143  rad]  (B.40) 

FOVs  =  2(8.172)  =  16.343°  [0.285  rad]  (B.41) 

Pas  90  ?7a  &sa 

=  90°  -  8.172°  -  20°  =  61.828°  [1.079  rad]  (B.42) 

Da  =  (6396)  [Sm(r6R118729?l  -  39,668  km  (B.43) 

sm(8.172) 


B.3  Disturbance  Environment 

Orbit  perturbing  forces  are  considered  negligible;  hence,  attitude  disturbances 
are  the  primary  concern.  The  satellite  structure  is  the  primary  load  path  and  de¬ 
signed  to  compensate  for  expected  minimum  stiffness  requirements  for  rigid  body 
purposes  [35,  page  2],  Due  to  alignment  with  the  positive  yaw  axis,  +2S,  platform 
disturbances  creating  roll  and  pitch  attitude  errors  affect  receiver  pointing.  Assum¬ 
ing  the  attitude  determination  and  control  system  (ADCS)  adequately  compensates 
against  expected  external  environmental  platform  disturbance  torques,  maintaining 
three-axis  stabilization,  satellite  internal  disturbances  affect  pointing  error.  In  par¬ 
ticular,  common  internal  disturbances  are  cm  uncertainty,  control  hardware  misalign¬ 
ment  or  output  error,  liquid  fuel  sloshing,  rotating  machinery,  and  thermal  shock  [16, 
page  355]. 
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Spacecraft  structural  loads  are  comprised  of  structurally  transmitted  vibration, 
shock,  and  thermal  stress  [11,  page  417].  Vibration  and  shock  induced  structural 
loads  result  from  translational  and  angular  accelerations  produced  by  control  hard¬ 
ware.  Thermal  stress,  due  to  differential  heating  or  cooling,  results  in  structural 
deformations. 

Control  hardware  misalignment  or  output  error  and  cm  uncertainty,  produce 
disturbance  torques  during  both  translational  and  rotational  motion  [16,  page  355]. 
These  disturbances  are  modeled  as  constant  disturbance  torques  during  thruster  br¬ 
ings.  Liquid  fuel  sloshing  and  rotating  machinery  effects  are  modeled  as  sinusoidal 
disturbance  torques.  Thermal  shocks  arise  from  thermal  stresses  on  solar  arrays  when 
the  satellite  enters  or  leaves  eclipse  [16,  page  355].  Assuming  each  solar  array  behaves 
as  a  cantilever  beam,  thermal  shocks  are  modeled  as  impulse  disturbance  torques. 

B.4  Summary 

This  appendix  presented  the  fundamental  satellite  orbital  and  attitude  dynam¬ 
ics  principles  which  impact  the  dual  FTCS  design.  Location  of  the  satellite  geosyn¬ 
chronous  orbit  maximizes  the  held  of  view  (FOV)  coverage  area  over  the  Earth’s 
surface  and  determines  the  relative  uplink  beam  propagation  distances.  It  is  assumed 
the  ADCS  provides  three-axis  stabilization  against  expected  environmental  distur¬ 
bances,  which  limits  the  platform  disturbances  to  be  from  internal  sources.  These 
internal  disturbances  are  considered  to  be  random  torque  inputs  defined  in  the  lin¬ 
earized  rotational  equations  of  motion.  Although  not  explicitly  implemented  in  the 
simulation,  these  equations  of  motion  are  the  basis  for  the  plant  disturbance  noise  af¬ 
fecting  motion  about  the  x-axis  and  y-axis,  as  shown  in  Chapter  III  and  Appendix  G. 


Appendix  C.  Aircraft  Dynamics 


C.l  Aircraft  Position 

As  illustrated  in  Figure  C.l,  aircraft  attitude  is  defined  about  either  a  true  North 
reference  frame  or  aircraft  body  axis  reference  frame.  Based  on  aircraft  location,  true 
North  reference  frame  is  defined  by  unit  vectors  xn,  yw,  zv:  xn  aligned  with  true 
North,  yw  faces  due  West,  zv  local  unit  vector  of  fa.  Unit  vectors  xa,  ya,  za,  roll, 
pitch,  and  yaw  axis  respectively,  define  the  aircraft  body  axis  reference  frame. 

Since  the  UAV  is  operating  within  the  satellite  receiver  field  of  view,  defined  in 
Appendix  D  as  FOVr,  and  assuming  a  spherical  Earth,  the  geocentric  longitude  and 
latitude  over  time  are 


Aa(^)  —  x>0{t  —  t0)  +  dc(t0)  +  5Xa(t)  (C.l) 

(paif)  =  (pel  +  S(pa(t)  (C.2) 


where  5Xa{t)  is  the  longitude  displacement  and  5(pa(t )  is  the  latitude  displacement. 
In  terms  of  the  GCI  frame,  defined  by  unit  vectors  i,  j,  and  k,  the  relative  position 
and  velocity  vectors  are  denoted  by 


ra 

Va 


r, 
dfa 
dt 


cos(A a)  cos (<pa)i  +  sin(Aa)  cos (4>a)j  +  sin {(pa)k 
=  Aaxra 


(C.3) 

(C.4) 


where  ra  is  the  position  vector  magnitude  and  ixa  is  the  aircraft  angular  velocity 
vector. 


C.2  Aircraft  Attitude  Dynamics 

C.2.1  Attitude  Dynamics.  Attitude  dynamics  about  the  aircraft  cm  is  de¬ 
fined  from  Euler’s  rotational  equation  of  motion  as 

K  +  u}axha  =  Ta  (C.5) 
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Figure  C.l:  Aircraft  Attitude 
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where  ha  is  the  aircraft  angular  momentum  vector  and  Ta  is  the  total  aircraft  torque 
vector.  Only  one  aircraft  body  axis,  ya,  is  aligned  to  a  principal  axis;  thus,  the  aircraft 
body  angular  momentum  vector  is  defined  by  [39,  page  374] 


ha  I  a^a 


hi  0  Ji3 


0  I22  0 

7;  1  0  /■ 


UJn 


(C.6) 


33  / 


where  Ia  is  the  aircraft  moment  of  inertia  matrix.  The  aircraft  angular  velocity  vector, 
uja.  is  equivalent  to  the  sum  of  the  body  axis  and  true  North  reference  frame  angular 
velocity  vectors  defined  respectively  by 


^xa 


^ ya 

\  Uza  ) 


^ ab  “t-  ^ nb 


(C.7) 


where  0ab  is  body  axis  reference  frame  angular  velocity  vector  and  a hb  is  the  true 
North  reference  frame  angular  velocity  vector  in  body  axis  coordinates. 


C.2.2  Attitude  Kinematics.  ffeading  angle,  (3a ,  and  angle  of  attack,  aa , 
determine  body  axis  reference  frame  orientation  relative  to  true  North  reference  frame. 
Attitude  with  respect  to  the  body  axis  reference  frame,  as  shown  in  Figure  C.l,  is 
defined  by  Euler  angles:  roll  angle  (ftb ,  pitch  angle  9b,  and  yaw  angle  hb-  Both  reference 
frames  are  related  to  the  GCI  frame  by  the  following  rotational  transformations 


Rna  =  Ri(180)  R2(—aa)  R3(— Pa) 

RNG  =  RMRs{-[  90-AJ) 

ftAG  _  rAN^NG  _  j^NG 

where  RNA  denotes  the  rotational  transformation  from  the  body  axis  reference  frame 
to  the  true  North  reference  frame,  RNG  denotes  the  rotational  transformation  from 
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the  GCI  frame  to  the  true  North  reference  frame,  and  RAG  denotes  the  rotational 
transformation  from  the  GCI  frame  to  the  body  axis  reference  frame. 

Incorporating  the  corresponding  direction  cosine  matrices,  RNA  and  RNG  are 
determined  as 


R 


NA 


^l(180)i?2(— 1 (y.a)R^,{— Pa) 
f  1  0  0  \  f  cos(aa)  0  —  sin(ci!0)  ^ 


/ 


V 


o-io 

0  0-1 

1  0  0 
0-10 
0  0-1 


0  1  0 

y  sin(aa)  0  cos(aa)  J 


(  cos(/3a)  -  sin (/3a)  0  ^ 
sin  {(5 a)  cos(/3a)  0 

0  0  1 


\  (  cos(aa)  cos(^a)  -  cos(aa)  sin(/3a)  -sin(aa)  ^ 


/ 


sin  (f3a)  COS  (/3a) 

\  sin(aa)  cos ((3a)  -  sin(aa)  sin(/3a 


^  cos(aa)  cos (13a)  -  cos (aa)  sin (/3a)  -  sin(aa)  ^ 

-  sin(/3a)  -  cos (f3a)  0 

y  -  sin(aa)  cos(/?a)  sin(o;0)  sin(/3a)  -  cos(aa) 


0 


cos 


iaa)  J 


(C.8) 
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Ri{^>a)Rz  (—  [90  —  Aa]) 

(  i  n  n  \  ( 
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0  0 

cos(0a)  sin(0a) 

\  0  —  sin(0a)  cos(0a)  J 

1  0  0  \  ( 

0  cos(0a)  sin  (0Q) 

\  0  —  sin(0a)  cos(0a)  J 

-  sin(Aa) 


V 


7 


V 


-  cos (0a)  cos(Aa)  -  cos(0a)  sin(Aa)  sin {(j)a) 
sin(0a)  cos(Aa)  sin(0a)  sin(Aa)  cos(0a) 


cos(Xa  —  90)  sin(Aa  —  90)  0  ^ 

—  sin(Aa  —  90)  cos(Aa  —  90)  0 
0  0  1 

—  sin(Aa)  cos(Aa)  0  ^ 

—  cos(Aa)  -  sin(Aa)  0 

0  0  1 

cos(Aa)  0  ^ 


(C.9) 
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Assuming  a  roll,  pitch,  and  yaw,  order  of  rotational  transformation,  the  body  axis 
reference  frame  angular  velocity  vector  is  defined  by 


tiab  =  R3(lpb)R2(0b)Rl((pb) 


0 

V  0  / 


■R3(A)R2(eb) 


/o\ 

ob 

V  0  / 


R3{i>b 


0  \ 

0 

V*/ 


The  true  North  reference  frame  angular  velocity  vector  is  determined  by, 


c on  —  R^^uje 

'  -  sin(Aa) 

-  COS (0a)  COs(Aa) 
y  sin  (0a)  cos(Aa) 

/  0  \ 

UJE  sin(0a) 
y  UE  COS (0a)  y 


COs(A a)  0  ^ 

-  cos (0a)  sin(Aa)  sin(0a) 
sin (0a)  sin(Aa)  cos(0a)  / 


/  0  \ 

0 

V  UE  ) 


(C.10) 


(C.ll) 


where  uje  is  the  Earth’s  angular  velocity  vector  and  u>e  =  7.289  x  10~5  rad/s.  In 
terms  of  the  aircraft  body  axis  reference  frame,  the  true  North  angular  velocity  vector 
becomes 


^ nb 


=  RagSe  =  (Rna)T  RngSe  =  (Rna) 


■)NA\ T  . 


UJr] 


(  cos(aa)cos (/3a)  —  sin(/3a)  -  sin(aa)  cos(0a)  ^ 


-  cos(aa)  sin(0a)  -  cos(0a)  sin(aa)  sin(0a' 

—  sin(aa)  0  —  cos(aa) 

'■  ~  sin (0a)  sin (0O)  -  sin(aa)  cos (0a)  cos(0a)  ^ 
=  UE  -  cos(0a)  sin (0a)  +  sin(aa)  sin(0a)  cos (0a) 

-  cos(aa)  COS (0a) 
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\ 


UJE  sin (0a) 

y  U)E  COS  (0a)  y 


(C.12) 
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Thus,  aircraft  attitude  equations  of  motion  can  be  represented  as 


tia  ~  ~  Ia  1^ala^a  +  I  a  1^a 


(C.13) 


where  f2a  is  the  angular  velocity  vector  skew  symmetric  matrix  defined  as, 


na  = 


o 

~Uza 


^ za 

0 


^ xa 


\ 


^x.a  0 


(C.14) 


C.3  Transmitter  Module 

It  is  assumed  that  all  transmitter  components  are  housed  within  a  telescope 
mounted  on  a  gimbaled  platform,  as  shown  in  Figure  C.2.  The  telescope  body  behaves 
as  a  cylindrical  axisymmetric  rigid  body,  defined  by  xt,  yt,  zt  principal  axes  and  body 
rates  9t,  'ipt-  The  outer  gimbal  is  aligned  with  the  za  axis  for  azimuth  control,  while 
the  inner  gimbal  is  aligned  with  the  yt  axis  for  elevation  control.  The  transmitted 
beam  is  aligned  with  the  xt  axis. 

C.3.1  Transmitter  Position.  In  terms  of  the  GCI  frame,  the  relative  trans¬ 
mitter  position  and  velocity  vectors  are  denoted  by 

r  i  T 


ft  =  rgt 

1  0  0 

(C.15) 

^  dft 

Vt  =  ~dt  = 

tit  x  ft 

(C.16) 

where  RGT  is  the  rotational  transformation  from  the  transmitter  principal  axis  frame 
to  the  GCI  frame  and  tit  is  the  transmitter  angular  velocity  vector. 

C.3. 2  Transmitter  Attitude  Dynamics.  Attitude  dynamics  about  the  trans¬ 
mitter  cm  is  defined  from  Euler’s  rotational  equation  of  motion  as 

ht  +  u>t  x  ht  =  Tt  (C-17) 
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Figure  C.2:  Transmitter  Platform 
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where  ht  is  the  transmitter  angular  momentum  vector  and  Tt  is  the  total  transmitter 
torque  vector. 

C.3.3  Transmitter  Attitude  Kinematics.  Transmitter  attitude  relative  to  the 
true  North,  aircraft  body  axis,  and  GCI  reference  frames  are  determined  respectively 
by, 


Rnt  =  R3(Xsa)R2(£sa )  (C.18) 

RAT  =  RAN RNT  =  (RNAjT  rnt  (C.19) 

Rtg  =  RtnRng  =  (Rnt)T  Rng  (C.20) 

where  RNT  denotes  the  rotational  transformation  from  the  transmitter  principal  axis 
reference  frame  to  the  true  North  reference  frame,  RAT  denotes  the  rotational  trans¬ 
formation  from  the  transmitter  principal  axis  reference  frame  to  the  aircraft  body  axis 
reference  frame,  and  RTG  denotes  the  rotational  transformation  from  the  GCI  frame 
to  the  transmitter  principal  axis  reference  frame. 

Incorporating  the  corresponding  direction  cosine  matrices,  RNT  is  determined 
as, 


R 


.NT 
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/  cos(Asa)  sin(Aso) 

-  sin(Aso)  cos(Asa)  0 
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COb^csay  U  blll^c 

0  1  0 
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-  sin(Asa)  cos(esa)  cos(Asa)  -sin(Asa)si 

'l  Q  nnafc 


\  -  sin(eso) 


a)  sin(esa) 

)  sin(eSQ) 
cos(ssa) 


(C.21) 
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C.3.4  Transmitter  Attitude  Equations  of  Motion.  Transmitter  angular  ve¬ 


locity  and  angular  momentum  vectors  are  determined  respectively  by 
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(C.23) 


where  \t  is  tlie  transmitter  moment  of  inertia  matrix. 

Thus,  transmitter  attitude  equations  of  motion  can  be  represented  as 


ujt  —  —If  1EltItuJt  +  If  lTt 


(C.24) 


where  Qt  is  the  transmitter  angular  velocity  vector  skew  symmetric  matrix  defined 


as, 


Elt  — 


\  ^ yt 


0  UJ  zt  U) yt 
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~^xt  0 


\ 


(C.25) 


C. 4  Disturbance  Environment 

Aircraft  perturbing  forces  are  defined  by  external  disturbance  forces  from  two 
flight  reference  regimes:  steady  and  unsteady  flight.  Steady  flight  is  classified  as 
nonaccelerated  and  nonrotating  flight  dynamics,  while  unsteady  flight  is  classified  as 
rotational  flight  dynamics  [39,  pages  222-223].  For  the  purposes  of  this  study,  emphasis 
is  placed  on  the  receiver  fine  tracking  performance  and  it  is  assumed  each  UAV  is  a 
rigid  body  in  the  steady  flight  regime  and  loitering  within  FOVr.  Hence,  aircraft 
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attitude  stabilization  and  constant  flight  parameters  are  controlled  by  an  autopilot 
[39,  page  307]. 

C-4-1  Pointing  Loss.  Assuming  the  autopilot  compensates  for  low  frequency 
external  disturbances  while  the  transmitter  gimbal  system  maintains  coarse  pointing, 
high  frequency  disturbances,  also  referred  to  as  jitter,  are  the  primary  concern  for  this 
study.  Previous  work  has  shown,  for  a  direct  detection  system,  pointing  jitter  induces 
optical  power  losses  of  -1.8  dB  [0.66]  and  -0.1  dB  [0.98]  respectively  for  1  yurad  and 
0.2  /irad  radial  pointing  error  [8,  page  256,259].  Thus,  based  on  the  tracking  error 
requirements  defined  in  Chapter  III,  it  is  assumed  the  pointing  loss  is 


Lt  «  0.8  [-0.969  dB] 


(C.26) 


At  the  satellite  receiver,  this  pointing  loss  is  assumed  to  create  random  beam  wander, 
assumed  to  have  zero-mean  white  Gaussian  statistics,  within  the  receiver  field  of  view 
(FOV).  As  such,  the  variance  of  the  radial  pointing  error  can  be  determined,  from 
the  pointing  loss  and  antenna  gain  equations  in  Chapter  II,  as  [8,  page  253] 


hi  (Lt)  In  (0.8) 

Gt  Gt 

-  In  (0.8)  ( ttAt ) 2  [rad2]  (C.27) 


Assuming  the  azimuth  and  elevation  errors  have  equal  variances,  the  corresponding 
azimuth  and  elevation  variance  respectively  is  determined  as 


In  (0.8)  /  XA2 
8vr2  \W0) 


[rad2] 


(C.28) 


Thus,  the  corresponding  azimuth  and  elevation  pointing  error  variance,  based  on  the 
Chapter  III  defined  transmitter  parameters,  is  calculated  respectively  for  each  beam 
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as 

(C.29) 
(C.30) 

C-4-2  Aerodynamic- Optical  Loss.  Aerodynamic-optical  effects  are  due  to 
turbulent  airflow  about  an  aircraft  inducing  air  density  fluctuations  resulting  in  phase 
distortion  of  the  transmitted  beam.  Phase  distortion  affects  on-axis  gain,  denoted 
by  the  Strehl  ratio,  which  relates  the  beam  irradiance  to  an  ideal  Gaussian  beam 
irradiance.  Thus,  aerodynamic-optical  loss  can  be  approximated  by  the  Strehl  ratio 
as  [15,  page  115] 

Lao  pc  Strehl  Ratio  =  exp(— kcae)2  (C.31) 

where  ae  is  the  rms  wavefront  error. 

For  wavelengths  from  visible  to  infrared  (IR),  with  an  aircraft  at  35,000  ft  and 
speed  Mach  0.85,  simulation  revealed  a  severe  strehl  ratio,  approximately  0.3  to  0.6, 
for  the  visible  spectrum,  while  IR  suffered  a  minor  strehl  ratio  near  unity  [32,  pages  6, 
11].  Based  on  the  chosen  laser  source  wavelengths  and  UAV  operational  parameters 
defined  in  Chapter  III,  it  is  assumed  the  strehl  ratio  is  0.7  [32,  page  11].  Thus,  the 
aero  optical  loss  is  determined  as 

Lao  ~  0.7  [-1.549  dB]  (C.32) 


G),  = 

0)B  = 


In  (0.8)  / 1510  x  10 


—9  \  2 


0.04 


87T2  y 

In  (0.8)  / 1550  x  10~ 


87T2  \ 


0.04 


=  4.0274  x  10~12  [rad2] 
=  4.2436  x  10~12  [rad2] 


C.  5  Summary 

This  appendix  presented  the  fundamental  flight  and  attitude  dynamics  for  the 
UAV  platform  and  transmitter  telescope.  Assuming  the  autopilot  maintains  steady- 
flight  conditions  and  coarse  tracking  is  maintained  by  the  transmitter  gimbal  mech¬ 
anism,  attitude  disturbances  are  considered  to  be  random  torque  inputs  defined  in 
the  transmitter  attitude  equations  of  motion.  Although  not  explicitly  implemented 
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in  the  simulation,  these  equations  of  motion  are  the  basis  for  the  pointing  jitter  af¬ 
fecting  motion  about  the  x-axis  and  y-axis,  as  shown  in  Chapter  Iff  and  Appendix  G. 
As  a  result,  pointing  jitter  induces  a  pointing  loss  of  approximately  0.8,  based  on 
zero-mean  white  Gaussian  statistics.  In  addition,  aerodynamic-optical  effects  impart 
a  power  loss  of  approximately  0.7. 
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Appendix  D.  Optical  Components 


D.  1  Transmitter 

D.1.1  Laser  Source.  Propagation  of  electromagnetic  (EM)  energy  per  unit 
area  per  unit  time,  within  a  medium  based  on  direction  zp,  is  described  by  the  Poynt- 
ing  vector  as  [13,  page  48] 


where 


S  =  c2e 


E  (^p,  t)  x  B  ( Zp,  t) 


c  =  speed  of  light 
3  x  108  m/s 

e  =  electric  permittivity 
E  =  electric  field  vector 
B  =  magnetic  field  vector 
t  =  time 


(D.l) 


The  electric  and  magnetic  field  vector  functions  are  related  by  satisfaction  of  partial 
differential  equations,  known  as  Maxwell’s  equations  [26,  page  159].  In  particular, 
E,  commonly  referred  to  as  the  optical  field,  is  used  to  define  propagation,  denoted 
by  [13,  page  50] 

E  (zp,t)  =  E0cos(iuct  +  4>c)  [V/m]  (D.2) 

where  E0  is  the  amplitude,  a ;c  is  angular  temporal  frequency,  and  </>c  is  the  phase. 
Angular  temporal  frequency  is  defined  by 

27TC 

ujc  =  - —  =  27 rfc  [rad/s]  (D.3) 

^ C 

where  Ac  is  the  wavelength  and  fc  is  the  temporal  frequency. 
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Irradiance,  equivalent  to  optical  power  per  unit  area,  is  proportional  to  the  time 
average  of  the  Poynting  vector  and  commonly  denoted  as  [13,  page  50] 

I=(S)t  =  ^  [W/m2]  (D.4) 


As  the  beam  propagates,  irradiance  is  subject  to  Fraunhofer  diffraction  as  a  function 
of  the  propagation  distance,  zp,  and  radial  distance,  pz,  from  the  propagation  axis, 
zp.  Based  on  the  defined  operational  parameters  in  Chapter  III,  the  laser  source  is 
assumed  to  be  a  lowest  order  transverse  electromagnetic  (TEM0o)  collimated  Gaussian 
beam.  Thus,  irradiance  is  defined  respectively  as  [26,  pages  85-86] 
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(D.7) 


where  IQ  is  the  peak  irradiance,  WQ  is  the  initial  beam  waist  radius,  W(zp)  is  the 
diffractive  beam  waist  radius  at  zp,  and  z0  is  the  half  the  depth  of  focus.  Irradiance  is 
modeled  as  a  Gaussian  function  of  pz  and  has  a  peak  value  at  pz  =  0  [26,  pages  83-85]. 
Hence,  the  irradiance  mean,  and  variance  are  determined  respectfully  as 
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VAR  [. I(p2 ,  zp)\  =  E  [ I{pz ,  zp)2]  -  E  [I(Pz,  zp)f 
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Transmitted  optical  power  is  determined  by  [26,  page  85] 


(D.9) 


P  =  5^  pW|)  [W] 


Thus,  in  terms  of  transmitted  optical  power,  irradiance  is  defined  as 
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As  a  function  of  propagation  distance,  the  peak  irradiance  is  denoted  by 
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(D.12) 


Based  on  the  Chapter  111  defined  operational  parameters  for  the  communi¬ 
cation  uplink,  where  Pt  =  100  mW,  Wa  =  40  mm,  esa  G  [77.3°,  90°],  and  h  G 
[18.3  km,  35786  km],  Figures  D.l  thru  D.3  illustrate  the  theoretical  free-space  prop¬ 
agation  characteristics  for  a  collimated  Gaussian  beam  at  arbitrary  Near  IR  carrier 
wavelengths,  Ac  G  [800  nm,  1100  nm,  1300  nm,  1500  nm].  Based  on  the  satellite  re¬ 
ceiver  field  of  view,  defined  later  in  this  section,  as  the  elevation  angle  increases  the 
uplink  beam  propagates  through  less  of  the  atmosphere.  However,  shorter  wavelengths 
suffer  less  diffraction  and  achieve  a  higher  free  space  peak  irradiance. 
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Figure  D.l:  Propagation  Distance  vs  Elevation  Angle 
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Figure  D.2:  Diffractive  Beam  Waist  Radius  vs  Elevation  Angle 
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Figure  D.3:  Transmitted  Peak  Irradiance  vs  Elevation  Angle 
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D.1.2  Optical  Modulation.  For  the  purposes  of  this  study,  binary  ampli¬ 
tude  shift  keying  (ASK),  also  known  as  on-off  keying  (OOK),  is  the  bandpass  digital 
modulation  format.  For  RF  systems,  ASK  is  a  continuous  wave  (CW)  amplitude 
modulation  technique:  varies  carrier  wave  amplitude  as  a  linear  function  of  the  base¬ 
band  waveform  [12,  page  114].  However,  in  optical  communications  the  modulation 
variable  is  optical  power,  which  is  assumed  to  be  linearly  proportional  to  an  input 
voltage  [22,  pages  149,  255].  Thus,  the  uncoded  NRZ-L  PCM  baseband  waveform, 
defined  in  Appendix  A  as  a  g(t )  and  shown  in  Figure  D.4,  pulses  the  laser  either  on 
or  off  transmitting  an  OOK  modulated  waveform,  s(t),  derived  from  the  optical  held 
equation  as  [12,  page  114] 


s(t) 


En 


1  + 


9(t) 

V 


cos  (a )ct  +  (pc) 


(D.13) 


Based  on  this  waveform,  we  can  assign  a  unique  transmission  symbol,  which 
has  a  corresponding  channel  symbol  meaning,  to  each  modulation  type.  Hence,  two 
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Figure  D.4:  OOK  Digital  Bandpass  Modulation 
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possible  unipolar  signal  waveforms  are  transmitted  within  the  interval  t  G  [0 ,Tg\, 
denoted  by 


Si(t)  =  2E0cos(u>ct  +  (j)c)  ==>-  symbol  1 
s2(t)  =  0  ==>-  symbol  0 


(D.14) 


Correlation  of  these  two  signals  reveals 


si(t)s2(t)dt 


[2E0cos(ujct  +  (f>c)\  [0 }dt 


0 


Thus,  the  transmitted  signal  waveforms  are  orthogonal  [29,  page  129]. 


(D.15) 


D.  1.2.1  Transmitted  Optical  Signal  Statistics.  Since  s(t)  is  a  linear 
function  of  g(t),  assuming  (j)c  is  an  independent  uniform  random  variable  distributed 
in  the  interval  [0,  27t],  the  mean  is  calculated  as  [10,  page  232] 
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(D.16) 
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The  autocorrelation  of  S(t )  is  determined  as 
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(D.17) 


The  autocovariance  of  S'(t)  is  defined  as  [10,  page  335] 


Cs(t ut2)  =  E  [{S'(ti)  -  ms(ti)}{S(t2)  -  ms(t2)}] 

=  E  [{5(40  -  0 }{S(t2)  -  0 }]=E  [S'(ti)S'(t2)] 

=  Rs(ti,t2)  —  Rs{t)  (D.18) 


Since  the  mean  is  constant  and  the  autocovariance  is  a  function  of  t  =  t-2  —  t\,  S(t ) 
is  a  wide-sense  stationary  process  (WSS)  [10,  page  359]. 

D.  1.2.2  Transmitted  Optical  Signal  Power  Spectral  Density.  The 
power  spectral  density  (PSD)  is  determined  from  the  Fourier  transform  of  Rs{t) 
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as  [10,  page  404] 


Ss(f)  =  {ccttWc)  F  {Rs(t)}  =  (e0CTrWa)  Rs{r)  exp(-j27r/r)dr 
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(D.19) 


Based  on  the  Chapter  III  defined  operational  parameters,  an  OOK  modulated 
PSD  for  a  1510nm  beam  is  shown  in  Figure  D.5,  accounting  for  positive  frequencies 
and  scaled  to  unit  amplitude.  The  null-to-null  bandwidth  of  s(t)  is  determined  as  [29, 
page  49] 
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The  average  power  is  determined  by  [1,  page  63] 
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Figure  D.5:  Transmitted  Signal  PSD 


Thus,  for  each  transmitted  binary  symbol,  the  average  power  is  either  R$( 0)  =  2 Pt 
for  symbol  1  or  R$( 0)  =  0  for  symbol  2.  Assuming  transmission  of  either  symbol  is 
equally  likely,  the  average  transmitted  power  is  Rs( 0)  =  Pt  [24,  page  472], 


D.1.3  Transmitter  Antenna  Gain.  As  defined  in  Chapter  2,  the  transmitter 
antenna  power  gain  for  a  Gaussian  beam  is  denoted  by  [8,  page  253] 


Gt 


2t tW0 


(D.22) 


Figure  D.6  illustrates  the  theoretical  relationship  between  antenna  gain  and  near 
infrared  (IR)  wavelengths.  Larger  Wa  yields  a  higher  gain,  however,  shorter  wave¬ 
lengths  achieve  higher  gain  than  longer  wavelengths.  On  the  other  hand,  as  shown 
in  Appendix  E,  a  larger  waist  radius  suffers  increased  atmospheric  loss  from  optical 
turbulence. 
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Near-IR  Wavelength  [nm] 

Figure  D.6:  Transmitter  Antenna  Gain  vs  Wavelength 

D.2  Receiver 

D.2.1  Optical  Signal  Processing. 

D.2. 1.1  Cassegrain  Telescope.  Based  on  telescope  dimensions,  the  re¬ 
ceiver  FOV  is  determined  from  the  effective  focal  length  and  primary  mirror  diameter, 
ft,  and  Dpm  respectively,  defined  as  [13,  pages  171,  217] 

FOVr  =  2  tan-1  (D.23) 

Focal  length  is  determined  from  the  telescope  focal  ratio,  //#,  which  is  determined 
as  [13,  page  174] 

//#  =  5  (D.24) 


121 


Effective  focal  length  and  primary  mirror  focal  length,  fpm,  are  related  by  the  mag¬ 
nification  ratio,  Mf,  which  is  denoted  by  [13,  page  221] 


fL  —  Mffpm  (D.25) 

Thus,  given  the  Chapter  III  defined  telescope  specifications,  300  mm  aperture  diam¬ 
eter,  15  focal  ratio,  and  10  magnification  ratio,  the  remaining  telescope  parameters 
are  calculated  as  [15,  pages  118-122] 

Dpm  =  300  mm 
/l  =  -Dpm//#  =  4500  mm 
fPm  =  h/Mf  =  450  mm 
DSm  0.25Dpm  7 o  mm 

L  =  ~/pm  =  150  mm 

where,  Dsm  is  the  secondary  mirror  diameter  and  L  is  the  length  of  the  telescope. 
Thus,  the  receiver  FOV  is  determined  as 

FOVr  =  3.818  deg  [0.066  rad]  (D.26) 

The  corresponding  Earth  coverage  about  the  subsatellite  point  (SSP),  defined 
as  the  Earth  central  angle,  is  determined  as  [16,  page  111] 

Ar  =  90°  -  rjar  -  er  (D.27) 
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Figure  D.7:  Receiver  Field  of  View  Geometry 


where  par  is  the  nadir  angle  visible  within  FOVr ,  and  er  is  corresponding  minimum 
transmitter  elevation  angle,  determined  respectively  as  [16,  page  111] 
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where  Pe  =  8.7°,  which  was  determined  in  Appendix  B,  is  the  radial  distance  of  Earth 
visible  from  the  satellite  orbit.  Thus,  as  illustrated  in  Figure  D.7,  the  Earth  central 
angle  is  determined  as 


Ar  =  90°  -  1.909°  -  77.278°  =  10.813  deg  [0.189  rad]  (D.30) 


The  corresponding  maximum  aircraft-to-satellite  range,  also  referred  to  as  propaga¬ 
tion  distance,  is  determined  as  [16,  page  111] 


2p(max) 
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ra  —( - r 
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(D.31) 


where  ra,  determined  in  Appendix  C,  is  the  magnitude  of  the  aircraft  position  vector. 
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D.2.2  Receiver  Antenna  Gain.  As  defined  in  Chapter  2,  the  receiver  antenna 
power  gain,  calculated  for  each  beam  respectively,  is  determined  as  [8,  page  253] 


(Gr)A 


W(0-3) 

1510  x  10-9 


3.8957  x  1011  [W] 


(D.32) 


(Gr)B 


W(0-3) 

1550  x  10-9 


3.6973  x  1011  [W] 


(D.33) 


D.2.3  Optical  Signal  Detection. 


D.2.3.1  Photodetector.  In  order  to  determine  the  received  power  for 
demodulation  into  the  communication  signal,  each  incident  beam  is  directed  to  a  cor¬ 
responding  avalanche  photodiode  detector  (APD).  A  photodiode  detector,  commonly 
referred  to  as  a  photodetector,  detects  and  converts  the  received  optical  power,  Pr, 
into  an  output  dc  current  determined  by  [17,  page  936] 


i  dc  RphPr  [A] 


(D.34) 


where  Rph  is  the  photodetector  responsivity.  Thus,  in  terms  of  the  two  transmitted 
unipolar  signal  waveforms,  within  the  interval  t  G  [0,77],  the  received  current  is 
determined  respectively  as 

h  =  RphPr  +  in  [A]  =>  symbol  1  ^ 

i2  =  in  [A]  =>■  symbol  0 

where  in  is  the  equivalent  noise  current.  The  photodetector  responsivity,  Rph,  is 
determined  by  [17,  page  937] 


Rph  =  Gd H  [A/W] 


(D.36) 
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where 


Gd  =  detector  gain 

r]e  =  quantum  efficiency 

qe  =  electronic  charge 

=  1.602  x  10~19  coulombs 

h  =  Plants  constant 

=  6.626  x  10~34  joules  •  s 

fd  =  detection  frequency 
c 

a; 


Responsivity  is  constant  while  operating  within  the  photodetector’s  linear  dy¬ 
namic  range  (DR):  region  between  noise  equivalent  power  (NEP)  and  saturation 
power,  Psat,  denoted  by  [17,  page  942] 

Dfl=101og(2|P)  [dB]  (D.37) 


Photodetector  noise  equivalent  power  (NEP),  defined  as  the  minimum  received  power 
for  a  SNR  of  unity,  is  determined  by  [17,  page  938] 


NEP  =  rms (in)  ( [W] 
\VeQeJ 


(D.38) 


Quantum  shot  noise  determines  the  maximum  capable  received  SNR  for  an  optical 
receiver,  which  defines  the  photodetector  saturation  limit  as  [1,  page  64] 

Psat  =  (D.39) 

^ qsh 

where  i2qsh  is  the  quantum  shot  current  noise  power.  Thus,  based  on  the  Chapter  III 
defined  operational  parameters  and  assuming  rje  =  0.8,  the  photodetector  responsivity, 
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operating  within  the  linear  dynamic  range,  is  calculated  for  each  beam  respectively 


as 


(0.8)(1.602  x  10“19) 


(6.626  x  10_34)(1.9868  x  1014) 
20  [0.9735]  =  19.4709  [A/W] 


(D.40) 


(0.8)(1.602  x  10“19) 


(6.626  x  10_34)(1.9355  x  1014) 
=  20  [0.9993]  =  19.9867  [A/W] 


(D.41) 


where  Gd(0)  is  the  dc  multiplicative  gain. 

Additional  photodetector  response  characteristics  are  rise  time,  fr,  fall  time,  tf , 
and  3  dB  cutoff  frequency,  fph .  Rise  time  is  the  time  response  from  10  to  90%  output 
peak  value,  while  fall  time  is  the  decay  interval  from  90  to  10%  peak  value.  Cutoff 
frequency,  which  defines  the  detection  bandwidth  BWdl  is  related  to  rise  time  and 
signal  bandwidth,  BWnn,  respectively  by  [17,  pages  943-944] 


0.35 


fph  =  BWd 


(0.886 )RhRvAr  [Hz] 


(D.42) 


Based  on  the  Chapter  111  defined  operational  parameters,  the  photodetector  gain 
frequency  response  is  calculated  for  each  beam  respectively  as  [1,  page  108] 


1/2 


20 


(D.43) 


(D.44) 
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Thus,  based  on  the  operational  parameters  defined  in  Chapter  III,  assuming 
symbol  1  is  sent  within  the  detector  bandwidth,  the  theoretical  received  optical  power 
for  each  beam  is  shown  in  Figure  D.8.  Although,  a  larger  wavelength  is  less  susceptible 
to  atmospheric  loss,  at  higher  elevation  angles  the  shorter  wavelength  achieves  slightly 
better  power  performance.  After  photodetection  and  conversion  to  a  current  signal, 
the  theoretical  received  dc  current  for  each  beam  is  illustrated  in  Figure  D.9.  In 
particular,  a  longer  wavelength  achieves  a  higher  dc  current  due  to  greater  detector 
responsivity. 


D.2.3.2  Quadrant  Cell  Detector.  In  order  to  determine  the  tracking 
error,  each  incident  beam  is  directed  to  the  corresponding  quadrant  cell  detector 
(QD)  field  of  view  (FOVd).  Similar  to  a  photodetector,  the  QD  detector  absorbs 
incident  power  from  an  optical  signal  and  produces  a  corresponding  current,  denoted 
by  [17,  page  936] 

V  =  RqdPr  [A]  (D.45) 

where  Rqd  is  the  quadrant  cell  detector  responsivity.  QD  responsivity  is  determined 
by  [17,  page  937] 

Rqd  =  Gd H  [A/W]  (D.46) 

Divided  into  four  equivalent  quadrant  cells,  as  shown  in  Figure  D.10,  current  signals 
from  each  cell  are  used  to  determine  azimuth  and  elevation  signal  errors.  In  the 
presence  of  random  quadrant  signal  noise,  rij,  assumed  to  have  zero-mean  Gaussian 
statistics,  error  signals  are  determined  by  [15,  pages  180-181] 


xe 

Ve 


Azimuth  Error  = 
Elevation  Error  = 


1 a  T  Fla  'lb  T  Fib  1c  Flc  Id  Fid 
ia  +  na  +  ib  +  nb  +  ic  +  nc  +  id  +  nd 

la  T  Fla  +  ?c  +  Tlc  lb  Fib  Id  Fld 

ia  +  na  +  ib  +  nb  +  ic  +  nc  +  id  +  nd 


(D.47) 

(D.48) 


Assuming  quadrant  signal  noises  are  independent  of  each  other,  the  mean  errors  are 


127 


1.48 


x  Id'9 


Figure  D.8:  Theoretical  Received  Optical  Signal  Power 


x  id'8 


Figure  D.9:  Theoretical  Received  DC  Current 
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Figure  D.10:  Quadrant  Cell  Detector 
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determined  by, 


E[xe\ 


E[ye } 


E  [ig  +  na  +  ib  +  nb  -  ic  -  nc  -  id  -  nd\ 

E  [ia  +  na  +  ib  +  nb  +  ic  +  nc  +  id  +  nd\ 
ia  +  ib  ~  ic~  id  +  E[na]  +  E[nb]  -  E[nc\  -  E[nd\ 
ia  +  ib  +  ic  +  id  +  E[na ]  +  E[nb]  +  E[nc\  +  E[nd] 
ia  R  ib  ic  id 

ia  +  ib  +  A  d~  id 

ia  d"  ic  ib  id 

d“  'I'b  d~  d~ 


(D.49) 

(D.50) 


Assuming  a  large  SNR  and  independent  quadrant  noises,  azimuth  and  elevation  vari¬ 
ances  are  calculated  as  [15,  pages  181-182] 


a 


2 

X 


a 


2 

y 


E[x2e }  -  E[xef 


E 

E 


{ig  +  na  +  ib  +  nb-ic-nc-id-  ndf 
{ia  +  na  +  ib  +  nb  +  ic  +  nc  +  id  +  ndf 

f  {ia  +  ib-ic~  id)  +  {na  +  nb-nc-  nd)  \ 
V  {ig  +  ib  +  ic  +  id)  +  {na  +  nb  +  nc  +  nd) ) 


f  ia  d ~  ib  ic  id  \ 

\ia  +  ib  +  ic  +  id) 

{ia  d-  ib  ic  id) 
{ia  +  ib  +  ic  +  id )2 


{ig  +  ib-ic-id)2  +  E[n2a]  +  E[n2b]  -  E[n2c]  -  E[n2d] 
{ia  d~  ib  d~  ic  d~  id)" 

E[nl}  +  E[nl}  +  E[nl}  +  E[nl] 


{ia  +  ib-ic-  id)2 

{ia  +  ib  +  ic  +  id)2 


{ia  d-  ib  d-  ic  d-  id) 

+  °d 

{ia  d-  ib  +  iC  d“  id)2 


i2 

°nq 

1 

(D.51) 

i2 

lqd 

SNRqd 

i2 

1 

i2 

Lqd 

SNRqd 

(D.52) 

where  iqd  is  the  current  signal  power  and  i2nq  is  the  current  noise  power,  also  referred  to 
as  current  noise  variance.  Since  azimuth  and  elevation  error  variances  are  equivalent, 
the  root  mean  square  (rms)  error  is  defined  by 


rms(<rT)  =  rms^)  =  arms 


1 

y/SNR^d 


(D.53) 
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In  addition,  quadrant  detector  noise  imparts  residual  angular  jitter,  known  as  noise 
equivalent  angle  (NEA),  defined  as  [15,  page  71] 


NEA  =  ^  = - -i - 

SF  SF^SNR^d 


(D.54) 


where  SF  is  the  slope  factor.  Thus,  the  accuracy  of  the  QD  sensor  is  inversely  pro¬ 
portional  to  SNRqd ,  which  affects  the  ability  of  the  tracking  control  system  to  reject 
platform  and  signal  propagation  disturbances  [5,  page  676]. 


D.2.3.3  Photodetector  Noise.  Three  types  of  noise  sources  occur 
within  a  photodetector:  thermal  noise,  dark  current  noise,  and  quantum  shot  noise. 
Thermal  noise,  also  known  as  Johnson  noise,  arises  from  electron  thermal  motion 
within  dissipative  elements  and  is  modeled  as  a  zero-mean  white  Gaussian  noise  pro¬ 
cess  [29,  pages  30-31].  Thermal  current  noise  power  is  denoted  as  [17,  page  932] 


VAR[ith] 


i2th. 

Rl 


[A2] 


(D.55) 


where  a2h  is  the  thermal  noise  power,  and  Rl  is  the  load  resistance.  The  thermal 
noise  power  is  determined  by  [17,  page  932] 


a2th  =  A{kBT°)BWd  [W] 


(D.56) 


where  kB  =  1.38  x  10”23  joules/°K  is  the  Boltzmann’s  constant,  and  T°  is  the  receiver 
noise  temperature.  Thus,  based  on  the  Chapter  III  defined  operational  parameters 
and  assuming  a  1000  Ohm  load  resistance,  thermal  current  noise  power  is  calculated 
as 


4(1.38  x  10"23) (400) (0.886) (200  x  103) 
1000 

=  3.9126  x  10“18  [A2] 


(D.57) 
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Due  to  the  inherent  nature  of  a  photodetector,  dark  current  is  an  additional 
current  source  present  whether  or  not  the  detector  is  illuminated.  Consisting  of  two 
components,  dark  current  is  modeled  as  [1,  pages  146] 

idk  =  idu  +  Gd(0)idm  [A]  (D.58) 

where  idu  is  the  unmultiplied  dark  current,  idm  is  the  multiplied  dark  current,  and 
Gd( 0)  is  the  DC  photodetector  gain.  For  the  purposes  of  this  study,  it  is  assumed 
both  the  unmultiplied  and  multiplied  dark  current  within  the  APD  are  approximately 
1  nA  [1,  page  116]. 

Quantum  shot  noise  occurs  within  a  photodetector  due  to  quantum  nature  of 
incident  photons  and  charge  carriers,  background  radiation,  and  dark  current.  Inci¬ 
dent  photons  arriving  upon  a  photodetector  exhibit  Poisson  statistics  [17,  page  929]. 
However,  assuming  a  high  photon  counting  capability,  photodetection  statistics  can 
be  approximated  by  a  zero-mean  white  Gaussian  distribution  [1,  page  240].  The  vari¬ 
ance,  referred  to  as  quantum  shot  noise  current  power,  is  determined  as  [17,  page  931] 

VAR[iqsh]  —  iqsh 

=  (^QeF„)  [idc  +  Gdiib  +  idm)\BWd  [A2]  (D.59) 

where  Fn  is  the  excess  noise  factor,  4  is  the  background  current,  and  idm  is  the 
multiplied  dark  current.  The  excess  noise  factor  can  be  approximated  by  [1,  page  143] 

Fn  «  Gd01  (D.60) 

Background  current  arises  from  optical  background  radiation,  which  causes 
background  current  flow,  4,  determined  as  [5,  page  679] 

n  =  (fg)  *2  [A]  (D.61) 
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where  c2  is  the  optical  background  noise  power,  defined  in  Appendix  E.  As  such, 
based  on  the  Chapter  III  defined  operational  parameters  and  assuming  r)e  =  0.8,  the 
optical  background  current  respectively  for  each  beam,  is  calculated  as 


(ib)A 


(0.8)(1.602  x  10“19) 
(6.626  x  10_34)(1.9868  x  1014) 
1.2816  x  10“14  [A] 


(1.3164  x  10“14) 


(D.62) 


{h)b 


(0.8)(1.602  x  10"19) 
(6.626  x  10~34)(1.93  55  x  1014) 
1.2816  x  10~14  [A] 


(1.2825  x  10“14) 


(D.63) 


Thus,  based  on  the  operational  parameters  defined  in  Chapter  III,  previously  defined 
dc  current,  dark  current  and  background  current,  assuming  symbol  1  is  sent  within 
the  detector  bandwidth,  the  maximum  theoretical  quantum  shot  noise  power  for  each 
beam  is  illustrated  in  Figure  D.ll.  In  particular,  as  the  elevation  angle  increases, 
which  decreases  the  propagation  distance,  the  noise  power  increases,  due  to  the  in¬ 
creased  dc  current.  In  addition,  the  quantum  shot  noise  power  is  at  least  three  orders 
of  magnitude  less  than  thermal  noise  power;  hence,  the  detector  is  thermal  noise 
limited  [17,  page  933]. 

During  normal  operations,  the  photodetector  equivalent  current  noise  power  is 
determined  as  [1,  page  154] 

=  ilsh  +  4  +  ZqelduBWd  [A2]  (D.64) 

Based  on  the  operational  parameters  defined  in  Chapter  III  and  previously  defined 
current  noise  sources,  the  theoretical  maximum  noise  equivalent  power,  respectively 
for  each  beam,  is  shown  in  Figure  D.12,  which  inherits  the  elevation  angle  versus 
power  trend  as  quantum  shot  current  noise  power. 
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Figure  D.ll:  Maximum  Quantum  Shot  Noise  Power 


x  10'18 


Figure  D.12:  Maximum  Photodetector  Noise  Equivalent  Current  Power 
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D.  3  Summary 

This  appendix  presented  the  fundamental  electromagnetic  theory  aspects,  prop¬ 
agation  characteristics,  and  optical  components  of  the  laser  communication  architec¬ 
ture.  A  near  infrared  (IR)  on-off  keying  (OOK)  modulated  Gaussian  beam  is  selected 
as  the  transmitter  carrier  wave,  which  is  subject  to  channel  propagation  effects  de¬ 
fined  in  Appendix  E.  Receiver  optical  components  collect  and  detect  each  incident 
beam  for  tracking  error  and  communication  signal  processing. 
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Appendix  E.  Communication  Channel 


The  communication  channel,  propagation  medium  between  transmitter  and  receiver, 
consists  of  Earth’s  atmosphere  and  free  space.  Although  a  majority  of  optical  dis¬ 
turbances  occur  within  the  atmosphere,  most  of  the  propagation  path  occurs  in  free 
space.  Within  free  space,  an  optical  signal  is  subject  to  free  space  diffraction,  which  is 
characterized  by  path  loss.  In  addition,  optical  background  radiation  from  extraneous 
sources  attenuates  the  received  signal  power. 

E.l  Atmospheric  Loss 

The  atmosphere  is  a  nonhomogeneous  medium,  which  affects  laser  beam  prop¬ 
agation  by  absorption,  scattering,  and  optical  turbulence.  These  effects  are  found  in 
four  types  of  atmospheric  conditions:  clear-air,  clouds,  rain,  and  fog  [4,  page  480]. 
Based  on  the  operational  parameters  defined  in  Chapter  III,  clear-air  atmospheric 
condition  is  the  primary  concern  for  this  study.  For  a  direct  detection  system,  clear- 
air  atmosphere  produces  the  following  optical  effects:  attenuation,  beam  spreading, 
beam  wander,  and  beam  scintillation  [2,  page  22].  Attenuation,  beam  spreading  and 
beam  wander  contribute  to  an  atmospheric  loss  to  the  received  optical  signal  power, 
while  beam  scintillation  affects  signal  probability  of  fade. 

E.  1.1  Attenuation.  Attenuation  is  caused  by  the  interaction  of  a  medium 
with  an  incident  electromagnetic  wave.  In  particular,  atoms  within  the  medium  in¬ 
teract  with  an  incoming  wave,  depending  upon  its  wavelength,  by  either  absorbing 
or  scattering  incident  photons  [13,  page  67].  A  majority  of  this  interaction  occurs 
up  to  20  km  above  the  Earth’s  surface,  however,  some  optical  turbulence  exists  up 
to  40  km  [4,  page  480].  Assuming  absorption  and  scattering  occur  independently, 
attenuation  loss  is  defined  by  the  Beer  Lambert  Bouguer  Law  as  [14,  page  25] 


(E.l) 
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where  /(/?-,  zp)  is  the  observed  irradiance,  I0(zp)  is  the  peak  irradiance,  pz  is  the  radial 
distance  from  the  propagation  axis  zp,  zp  is  the  propagation  distance,  and  a(z)  is  the 
extinction  coefficient.  If  the  extinction  coefficient  is  assumed  constant,  for  a  given 
wavelength,  Ac,  attenuation  loss  becomes 

Lat  =  exp  [-a(Ac)zp]  (E.2) 

Based  on  clear-air  atmospheric  conditions,  which  is  a  low  attenuation  environment, 
optical  and  infrared  (IR)  laser  beam  propagation  is  predominantly  affected  by  optical 
turbulence  [4,  page  480].  Thus,  for  the  purposes  of  this  study,  the  attenuation  loss  is 
approximated  as 

Lat  «  0.9  [-0.46  dB]  (E.3) 

E.1.2  Optical  Turbulence.  Optical  turbulence  is  defined  as  atmospheric  re¬ 
fractive  index  fluctuations  caused  by  temperature,  which  induces  beam  spreading, 
beam  wander,  and  beam  scintillation  upon  a  propagating  optical  wave.  For  the  pur¬ 
poses  of  this  study,  optical  turbulence  effects  are  based  on  the  Kolmogorov  power 
law  spectrum  [4,  pages  490].  The  strength  of  optical  turbulence  is  determined  by  the 
refractive  index  structure  constant,  C2,  which  is  assumed  to  be  constant  along  a  hor¬ 
izontal  path.  However,  for  propagation  along  a  vertical  or  slant  path,  the  structure 
constant  varies  as  a  function  of  altitude,  commonly  defined  by  the  Hufnagel- Valley 
(H-V5/7)  profile  model  as  [4,  pages  480-481] 

C2n{h)  =  0.00594(u/27)2(10“5h)10exp(— h/1000) 

+  2.7  x  10~16  exp(— h/1500)  +  Ah  exp(-h/100)  (E.4) 

where  h  is  the  altitude,  v  —  21  m/s  is  the  rms  windspeed,  and  A} l  =  1.7  x  10-14  m~2'3 
is  the  nominal  value  of  C'/(0)  at  the  ground.  Figure  E.l  illustrates  the  structure 
constant  profile  based  on  the  H-V5/7  model.  Due  to  operational  parameters  defined 
in  Chapter  III,  the  altitude  range  being  considered,  h  G  [18.3  km, 40  km],  results 
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Figure  E.l:  Refractive  Index  Structure  Constant  Profile 


in  an  extremely  small  structure  constant,  as  illustrated  in  Figure  E.2  For  a  com¬ 
munication  uplink,  C^(h')  is  used  to  determine  the  following  optical  turbulence  con¬ 
stants  [4,  pages  502,523] 


A*o 

Ah 

/J'2 

lh 

Z 


C2n(h)dh 


c2n(h) 


0-0 


h  -  ha 
hs  ha 


5/3 


dh 


Re 


1  - 


Cl(h)  (if3dh 

f  '  Cl(h)  {y/s  [A?  +  i  (1  -  0£)]5/a  -  A5/y5/sJ5/3(iA 


h  -  ha 
K  -  hn 


(E.5) 

(E.6) 

(E.7) 

(E.8) 


where  ha  is  the  aircraft  altitude,  and  hs  is  the  satellite  altitude.  Parameters  0,  0, 
and  A  are  the  output  plane  beam  parameters,  which  where  derived  in  Appendix  D 
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Figure  E.2:  Relative  Refractive  Index  Structure  Constant  Profile 
as  [4,  page  489] 


0  =  00 

(E.9) 

0  =  1-0 

(E.10) 

A  =  2% 

kcW2(zp 

y  (E.ll) 

where  zp  is  the  propagation  distance,  W(zp)  is  the  diffractive  beam  radius,  and  kc 

is  the  wave  number.  Parameters  0O  and  A0  are 

the  input  plane  beam  parameters, 

defined  in  Appendix  D  as  [4,  page  488] 

©o  =  1 

(E.12) 

A  -  2Zp 
kcW2 

(E.13) 

where  Wa  is  the  initial  beam  radius  at  transmission. 
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E.  1.2.1  Beam  Spreading.  Optical  turbulence  produces  additional 
beam  spreading,  beyond  normal  diffraction  occurring  within  free-space  propagation. 
Thus,  for  a  Gaussian  beam  the  effective  beam  radius  is  defined  by  [4,  page  500] 


We{zp) 


W(zp)  [1  +  (D0/rat)^] 1/2  ,  0  <  Da/rat  <  1 

W(zp)  [1  +  (D0/ratf/3]3/5  ,  0  <  D0/rat  <  oo 


where  D2  =  8 W2  and  rat  is  the  atmospheric  coherence  width,  defined  as  [4,  page  492] 


Tat  =  [0.42  SeC  (CasJ^c/d  3/5 


(E.15) 


where  (as  is  the  zenith  angle,  defined  as  the  complement  to  the  elevation  angle.  The 
first  condition  satisfies  the  condition  for  weak  irradiance  fluctuations  theory,  while  the 
second  condition  satisfies  strong  irradiance  fluctuation  theory.  As  similarly  derived 
in  Appendix  D,  the  expected  Gaussian  beam  irradiance  profile,  subject  to  We(zp),  is 
determined  as  [4,  page  494] 


E[I(pz,zp)\ 


W2 


I  K2(zP) 


exp 


2  Pi 


[  We2(zp) 


(E.16) 


where  pz  is  the  radial  distance  from  the  propagation  axis.  Comparison  with  free  space 
peak  irradiance,  defined  in  Appendix  D,  the  beam  spreading  power  loss  is  defined  as 


Lhc 


E[I(pz,zp )] 

^o(^p) 


w2(zp) 

W? 


exp 


2  pi 


W2(zp)  W2J_ 


(E.17) 


where  the  maximum  power  loss  occurs  on-axis,  pz  =  0,  denoted  by 


Lf,s  (max ) 


e\H^zp)\ 

lo(Zp) 

~w2(zpy 

.  ^e2  . 


W2(zp) 

w2 


exp 


(0) 


W2(zp) 


(E.18) 


Figure  E.3  illustrates  the  theoretical  irradiance  fluctuation  condition  of  beam  spread¬ 
ing  at  the  satellite  receiver,  for  a  collimated  Gaussian  beam  based  on  Chapter  III  input 
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Figure  E.3:  Irradiance  Fluctuation  Condition 


parameters  Pt  =  100  mW,  Wa  =  40  mm,  esa  e  [77.3°,  90°],  h  e  [18.3  km,  35786  km], 
and  arbitrary  Near  IR  carrier  wavelengths,  Ac  G  [800  nm,  1100  nm,  1300  nm,  1500  nm]. 
For  a  given  initial  beam  radius,  shorter  wavelengths  approach  the  strong  irradiance 
fluctuation  regime,  while  long  wavelengths  are  less  susceptible  to  optical  turbulence. 
On  the  other  hand,  the  corresponding  beam  spreading  effects,  shown  in  Figures  E.4 
thru  E.6,  share  a  similar  trend  with  free-space  diffraction  effects  revealed  in  Ap¬ 
pendix  D.  Shorter  wavelengths  are  less  susceptible  to  beam  spreading,  which  results 
in  a  higher  mean  peak  irradiance. 


E.  1.2.2  Beam  Wander.  Optical  turbulence  induces  random  radial 
displacements,  based  on  zero-mean  Gaussian  statistics,  of  the  incident  beam  centroid 
within  the  receiver  boresight,  commonly  referred  to  as  beam  wander.  For  zenith  angles 
less  than  60  degrees,  optical  turbulence  is  characterized  by  weak  irradiance  fluctuation 
theory  based  on  Rytov  approximation  [4,  page  478].  As  such,  based  on  the  previously 
defined  operational  parameters,  the  beam  wander  displacement  variance  for  an  uplink 
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Figure  E.4:  Effective  Beam  Radius 
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Figure  E.5:  Mean  Peak  Irradiance 
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Figure  E.6:  Beam  Spreading  Loss 


propagation  path  is  defined  as  [4,  page  502] 


VAR[Pbw]  =  [7.25Po(hs-ha)2sec3(Cas)W;1/3] 


h  fm2l 


(E.19) 


=  cr, 


bw 


After  conversion  to  radians,  pbw  =  crbw/zp ,  beam  wander  power  loss  can  be  interpreted 
as  a  pointing  loss,  determined  by  [8,  page  253] 

Lbw  =  exp  (~Gtp2bw)  (E.20) 

Figures  E.7  and  E.8  illustrate  the  observed  theoretical  beam  wander  effects.  Although 
beam  wander  is  independent  of  wavelength,  a  similar  trend  results  for  increasing 
elevation  angle  as  previously  shown  for  beam  spreading,  ffowever,  the  power  loss 
is  affected  by  the  signal  gain,  which  is  a  function  of  wavelength  as  illustrated  in 
Appendix  D.  Thus,  shorter  wavelengths  suffer  more  power  loss  due  to  beam  wander. 
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Figure  E.7:  Root  Mean  Square  Beam  Wander 


Figure  E.8:  Beam  Wander  Loss 
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E.  1.2.3  Beam  Scintillation.  Random  index  of  refraction  fluctuations, 


caused  by  optical  turbulence,  produce  corresponding  irradiance  fluctuations,  known 
as  beam  scintillation.  Classified  as  weak,  medium,  or  strong,  irradiance  fluctuations 
are  analyzed  by  the  normalized  variance  of  irradiance,  commonly  known  as  scin¬ 


tillation  index,  denoted  by  erf  [4,  page  261].  Based  on  the  previously  defined  uplink 


communication  operational  parameters  and  ideal  tracking,  within  the  weak  irradiance 
fluctuation  regime,  aj  is  defined  by  [4,  page  504] 


(E.21) 


where  U  ( pr  —  pbw )  is  the  unit  step  function. 

Ideally,  a]  <  1  satisfies  the  conditions  for  weak  fluctuation  theory  [4,  page  264], 
Based  on  the  previously  defined  parameters  for  the  uplink  beam,  the  tracked  on-axis, 
pr  =  0,  scintillation  index  is  shown  in  Figure  E.9.  Sharing  a  similar  trend  for  elevation 
angle  as  other  optical  turbulence  effects,  transmitting  from  an  altitude  of  18.3  km, 
results  in  a  scintillation  index  roughly  between  an  order  of  magnitude  of  10~4'38  and 
10-4  28.  Thus,  optical  turbulence  induced  beam  scintillation  effects  can  be  considered 
negligible  beyond  an  altitude  of  20  km  [4,  page  504], 

E.1.3  Atmospheric  Loss.  Atmospheric  loss  is  comprised  of  power  losses  from 
attenuation,  beam  wander,  and  beam  spreading,  denoted  by 


La  Lat  T  Lbs  T  L iru-  [dB] 


(E.22) 


Based  on  the  precalculated  values,  the  maximum  theoretical  atmospheric  loss  is  shown 
Figure  E.10.  Retaining  the  trend  in  elevation  angle  and  wavelength  dependency, 
longer  wavelengths  are  less  susceptible  to  power  loss  induced  by  atmospheric  effects. 
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Figure  E.9:  Tracked  On-Axis  Scintillation  Index 


Figure  E.10:  Atmospheric  Loss 
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E.2  Path  Loss 


As  shown  in  Appendix  D,  the  decrease  in  irradiance  of  an  electromagnetic  wave, 
subject  to  free-space  diffraction,  is  a  function  of  propagation  distance.  As  such,  the 
effective  transmitted  power  loss,  commonly  referred  to  as  path  loss,  is  defined  by  [29, 
page  254] 


Based  on  the  previously  defined  parameters  for  the  uplink  beam,  the  path 
shown  in  Figure  E.ll.  Although  wavelength  dependent,  clearly  path  loss  is  the 
contributor  to  attenuation  of  transmitted  signal  power. 


(E.23) 

loss  is 
largest 


CD 

T3 


CO 

CO 

o 


0) 


:§; 

o 


Q_ 


-289  | - , - . - . - . - . - . - 

-290  - 

-291  -  . _______ 

- 800 nm 

- lOOOnm 

'292  '  1300nm  ' 

-  —  1500nm 

-293  -  _ _ - 

-294  - 

-295  -  - - 

-296  - 1 - 1 - 1 - 1 - 1 - 1 - 

76  78  80  82  84  86  88  90 


- 800nm 

- lOOOnm 

1300nm 
-  1500nm 


Elevation  Angle  [deg] 
Figure  E.ll:  Path  Loss 


E.3  Optical  Background  Radiation 

Extraneous  optical  sources,  such  as  the  sun,  planets,  stars,  and  blackbody  ra¬ 
diation,  are  the  main  contributors  of  optical  background  radiation.  Depending  upon 
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wavelength  and  the  satellite  receiver  field-of-view  (FOV),  these  sources  are  modeled 
as  additive  white  Gaussian  noise  (AWGN)  for  a  direct  detection  optical  receiver  [4, 
page  444],  However,  optical  filtering  can  restrict  background  radiation  to  a  manage¬ 
able  level,  at  the  cost  of  signal  attenuation.  For  an  Earth  observing  satellite  receiver 
in  geosynchronous  orbit  (GEO),  optical  background  radiation  is  approximately  104 
photons  per  second  per  Angstrom  of  optical  bandwidth  [1,  page  151].  Depending  upon 
optical  filter  bandwidth,  denoted  as  AXp,  the  observed  optical  background  radiation 
power  is  defined  by  [5,  page  679] 


.2 


(Nb)(AXF)  [W] 


(E.24) 


where 


[W/m] 


(E.25) 


h  =  Plank's  constant 


=  6.626  x  10  31  joules  •  s 
c  =  speed  of  light 


«  3  x  10s  m/s 

Ac  =  carrier  signal  wavelength 
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Based  on  the  operational  parameters,  defined  in  Chapter  III,  the  optical  background 
noise  power,  respectively  for  each  beam,  is  calculated  as 


(< ?l)A  =  (IxlO14) 


(6.626  x  1CT34)(3  x  108) 
1510  x  10-9 
1.3164  x  10"14  [W] 


(1  x  10 


-9> 


(E.26) 


Mb  =  (1  x  10  ) 


(6.626  x  10_34)(3  x  108) 


1550  x  10 


-9 


(1  x  10“9) 


=  1.2825  x  10 


-14 


[W] 


(E.27) 


E.4  Summary 

This  appendix  presented  the  channel  propagation  effects  encountered  within  the 
laser  communication  architecture.  Assuming  clear-air  atmospheric  conditions,  path 
loss  is  the  greatest  contributor  to  signal  attenuation.  However,  optical  turbulence 
induces  additional  signal  attenuation  due  to  beam  spreading  and  beam  wander  effects. 
Optical  background  radiation  creates  an  additional  noise  source  for  signal  detection 
processing. 
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Appendix  F.  Controller  Characterization 

As  determined  in  Chapter  III,  the  control  system  is  required  to  regulate  radial  tracking 
error  within  the  error  performance  limit,  1.68  /rrad  and  1.74  ^nrad  respectively  for  each 
beam,  to  ensure  the  communication  system  meets  a  bit  error  rate  (BER)  of  10~6. 
On  the  other  hand,  due  to  a  single  actuator  limitation  of  one  fast  steering  mirror 
(FSM),  the  controller  must  switch  between  state  estimates  for  each  beam.  Thus,  the 
controller  characterized  by  the  weighting  matrices  and  switching  time  which  yields 
the  best  performance. 

Analysis  of  controller  performance  is  based  on  applying  Monte  Carlo  analysis, 
taking  samples  generated  from  multiple  simulation  runs,  and  computing  an  ensemble 
average  of  the  outputs.  In  order  perform  ensemble  averaging,  the  system  model  is 
converted  to  a  fixed-step  system  using  the  SIMULINK  Runge-Kutta  solver  solution. 
After  trial  and  error,  a  fixed  time  step  of  Ts  =  5  x  ICC5  sec  was  used  to  ensure 
against  SIMULINK  processing  errors.  Only  ten  simulation  runs  were  performed  to 
ensure  against  creating  memory  overloads,  due  to  the  small  fixed  time  step  and  final 
time  selection,  in  processing  the  data  in  MATLAB. 

In  addition,  a  rate  transition  block,  defaulted  to  act  as  a  unit  delay  for  con¬ 
tinuous  sample  time,  is  used  on  all  plant  disturbance  and  measurement  noise  inputs 
to  ensure  proper  timing  of  these  inputs  into  the  system.  Correlation  time  for  for  the 
bandlimited  white  noise  blocks,  due  to  SIMULINK  requirements,  is  changed  to  a  inte¬ 
ger  factor  of  the  fixed  time  step  as  tc  =  20 Ts  =  0.001  sec.  The  response  signal  means 
and  standard  deviation  statistics  are  computed  to  produce  outputs  for  evaluating  the 
modeled  fine  tracking  control  system. 

F.l  Weighting  Matrices 

Since  each  quadrant  cell  detector  determines  relative  azimuth  and  elevation 
tracking  error,  an  initial  state  weighting  matrix  based  on  this  tracking  error  is  given 


150 


by 
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As  derived  in  Chapter  III,  the  corresponding  system  state  weighting  matrix  is  denoted 
by 
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(F.2) 


The  fast  steering  mirror  (FSM)  is  modeled  as  an  uncoupled  second  order  SISO  system 
for  each  axis,  thus  the  control  weighting  matrix  is  defined  by 


R  = 


1  0 
0  1 


(F.3) 


In  order  to  determine  the  relative  LQR  weighting  matrix  affects  on  controller  perfor¬ 
mance,  the  control  system  was  implemented  by  varying  Q  and  R  independently  and 
evaluated  against  the  tracking  error  performance  limit  with  the  maximum  plant  dis¬ 
turbance  due  to  maximum  rms  beam  wander  of  5.406  x  10-6  /irad.  The  implemented 
weighting  matrices  are  shown  in  Table  F.l,  and  the  scalar  coefficients  of  Q  and  R  will 
be  referred  to  as  the  weighting  scale  factors. 
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Table  F.l:  Weighting  Matrix  Variables 


State  Matrix 

Ql 

Q2 

Q3 

Q4 

Scale  Factor 

1000 

100 

10 

1 

Control  Matrix 

R1 

R2 

R3 

R4 

Scale  Factor 

1 

0.1 

0.01 

0.001 

F.1.1  State  Weighting  Matrix  Response.  Based  on  simulation  final  time, 
Tf  =  1  second,  and  an  arbitrary  switching  time  based  on  the  FSM  2%  settling  time, 
Tsw  =  0.4  seconds,  tracking  error  and  control  input  data  was  collected  for  each  state 
weighting  matrix,  associated  with  the  chosen  tracking  error  weighting  matrix,  paired 
with  a  unit  control  weighting  matrix.  Figures  F.l  thru  F.3  show  a  magnified  view 
of  the  time  response  ensemble  average  mean  and  standard  deviation  data.  Based 
on  the  time  interval  illustrated,  Beam  B  tracking  error  is  regulated  while  the  other 
beam  is  unregulated.  In  particular,  the  standard  deviation  transient  peaks,  during 
regulation  of  Beam  B,  are  due  to  the  noise  correlation  time.  An  initial  comparison  of 
the  response  for  each  state  weighting  scale  factor  response  reveals  relatively  similar 
mean  values.  However,  a  larger  state  weighting  matrix  encounters  more  dispersion 
as  seen  in  the  standard  deviation  plots.  In  terms  of  control  input  response,  which 
reveals  similar  behavior  for  each  axis,  a  larger  state  weighting  results  in  greater  mean 
control  authority,  but  at  the  cost  of  increased  dispersion. 

Taking  a  temporal  average  of  the  ensemble  averaged  mean  and  standard  devi¬ 
ation  data,  illustrated  in  Figures  F.4  thru  F.6,  reveals  each  state  weighting  matrix 
yields  relatively  similar  temporal  mean  tracking  error  values,  but  higher  order  weight¬ 
ing  has  slightly  greater  dispersion.  Analysis  of  the  temporal  standard  deviation  data 
shows  an  increasing  temporal  mean  trend  for  lower  order  weighting,  but  the  disper¬ 
sion  decreases.  A  larger  control  weighting  results  similar  in  temporal  mean  control 
and  larger  dispersion  due  to  more  regulation  being  applied.  In  addition,  a  decreasing 
mean  and  dispersion  trend  in  control  standard  deviation  is  seen  for  decreased  weight- 


152 


5 
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Figure  F.l:  State  Weighting  Matrix  Ensemble  Average  Mean  Tracking  Error  Re¬ 

sponse 

ing.  Based  upon  these  results,  the  optimum  state  weighting  matrix  is  Q3,  due  to 
better  overall  standard  deviation  performance. 


F.l. 2  Control  Weighting  Matrix  Response.  Applying  the  same  method  used 
for  determining  the  state  weighting  matrix,  tracking  error  and  control  input  data 
was  collected  for  for  each  control  weighting  matrix  paired  with  a  unit  tracking  error 
matrix.  Figures  F.7  thru  F.9  show  a  magnified  view  of  the  time  response  ensemble  av¬ 
erage  mean  and  standard  deviation  data.  In  particular,  a  greater  tracking  error  trend 
is  evident  for  decreased  control  weighting.  Thus,  lower  control  weighting  increases 
control  authority,  at  the  cost  of  greater  deviation.  Temporal  average  of  the  ensemble 
averaged  mean  and  standard  deviation  data,  illustrated  in  Figures  F.10  thru  F.12, 
reveals  each  control  weighting  matrix  yields  a  relatively  similar  temporal  mean  track¬ 
ing  error  response.  However,  comparison  of  the  standard  deviation  response  shows  a 
lower  order  control  weighting  has  slightly  less  mean  deviation  at  the  cost  of  greater 
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Figure  F.2:  State  Weighting  Matrix  Ensemble  Average  Tracking  Error  Standard 

Deviation  Response 


X-axis  Y-axis 


Figure  F.3:  State  Weighting  Matrix  Ensemble  Average  Control  Response 
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Figure  F.4:  State  Weighting  Matrix  Temporal  Mean  Tracking  Error  Response 
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Figure  F.5:  State  Weighting  Matrix  Temporal  Standard  Deviation  Response 
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Figure  F.6:  State  Weighting  Matrix  Temporal  Control  Response 


dispersion.  As  seen  with  the  state  weighting  matrix  data,  each  control  weighting  case 
achieves  a  similar  mean  control  input  response,  but  smaller  weighting  results  in  larger 
deviation  due  to  less  regulation  being  applied.  Based  upon  these  results,  the  optimum 
control  weighting  matrix  is  R2,  due  to  better  overall  standard  deviation  performance. 


F.2  Switching  Time 

Applying  a  similar  simulation  approach  to  characterize  switching  time  effects, 
the  control  system  was  implemented,  with  the  predetermined  optimum  weighting 
matrices  and  Tsw  was  varied  and  evaluated  against  the  tracking  error  performance 
limit.  The  implemented  switching  times,  based  on  multiples  of  the  FSM  2%  settling 
time  of  approximately  0.002  seconds,  are  shown  in  Table  F.2.  Applying  the  same 
method  used  for  determining  the  weighting  matrices,  tracking  error  and  control  input 
data  was  collected  for  each  switching  time  over  a  simulation  time  of  5  seconds. 
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Beam  A  Beam  B 


Figure  F.7:  Control  Weighting  Matrix  Ensemble  Average  Mean  Tracking  Error 

Response 


Beam  A  Beam  B 


Figure  F.8:  Control  Weighting  Matrix  Ensemble  Average  Tracking  Error  Standard 
Deviation  Response 
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Figure  F.9:  Control  Weighting  Matrix  Ensemble  Average  Control  Response 
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Figure  F.10:  Control  Weighting  Matrix  Temporal  Mean  Tracking  Error  Response 
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Control  Weighting  Matrix  Temporal  Standard  Deviation  Response 
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Figure  F.12:  Control  Weighting  Matrix  Temporal  Control  Response 
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Table  F.2:  Candidate  Switching  Time  Variables 


Variable 

T1 

T2 

T3 

T4 

T 

1  sw 

0.002  sec 

0.02  sec 

0.2  sec 

1.0  sec 

F.2.1  Ensemble  Average  Tracking  Error  Response.  Since  each  FSM  axis 
is  modeled  by  the  same  SISO  system,  controlled  independently,  and  yielded  similar 
weighting  performance,  the  tracking  error  time  response  results  are  only  given  for 
the  X-axis,  as  illustrated  in  Figures  F.13  thru  F.16.  Each  figure  shows  the  ensemble 
mean  and  standard  deviation  results  for  the  time  domain  propagated  through  to  each 
switching  time.  As  time  progresses,  the  controller  switches  tracking  error  regulation 
between  Beam  A  and  Beam  B.  In  either  case,  at  time  of  regulation  for  each  switching 
time,  an  initial  transient  spike  is  seen  in  the  mean  response.  In  particular,  each 
standard  deviation  time  response  exhibits  initial  transients  occurring  at  intervals  of 
the  noise  block  correlation  time.  Based  on  these  results,  it  is  apparent  that  all  of  the 
switching  times  have  similar  performance. 

F.2. 2  Ensemble  Average  Control  Input  Response.  Comparison  of  the  en¬ 
semble  mean  and  standard  deviation  of  the  control  input,  shown  in  Figures  F.17  thru 
F.20,  reveals  that  all  switching  times  generate  a  similar  time  response.  As  the  time 
domain  propagates,  regulation  standard  deviation  is  subject  to  the  noise  block  cor¬ 
relation  time  interval,  as  seen  in  the  tracking  error  response.  However,  none  of  the 
candidate  switching  times  show  noticeable  benefits  over  other  switching  times. 

F.2. 3  Temporal  Average  Response.  Temporal  average  of  the  ensemble  aver¬ 
aged  mean  and  standard  deviation  data,  illustrated  in  Figures  F.21  thru  F.23,  reveals 
each  switching  time  case  yields  relatively  similar  temporal  mean  tracking  error  re¬ 
sponse.  In  addition  the  standard  deviation  share  similar  range  of  dispersion  about 
the  mean.  The  slight  increase  at  T4  in  Beam  A  mean  deviation,  while  Beam  B  de¬ 
creases,  is  due  to  more  Beam  B  tracking  error  regulation  occurring  over  the  5  second 
simulation.  However,  the  range  of  dispersion  about  the  mean  deviation  does  not 
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Beam  A  Beam  B 


Figure  F.13:  Tracking  Error  Time  Response,  te  (0,0.01s) 

Beam  A  Beam  B 


Time  [sec]  Time  [sec] 

Figure  F.14:  Tracking  Error  Time  Response,  te(0. 02, 0.03s) 
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Beam  A  Beam  B 


Time  [sec]  Time  [sec] 

Figure  F.15:  Tracking  Error  Time  Response,  t£(0. 20, 0.21s) 


Beam  A 


Beam  B 


Time  [sec]  Time  [sec] 

Figure  F.16:  Tracking  Error  Time  Response,  te(l. 00, 1.01s) 
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Beam  A 


Beam  B 


Figure  F.17:  Control  Input  Time  Response,  t€ (0,0.01s) 

Beam  A  Beam  B 


Figure  F.18:  Control  Input  Time  Response,  te (0.02,0.03s) 
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Beam  A 


Beam  B 


Figure  F.19:  Control  Input  Time  Response,  te (0.20,0.21s) 

Beam  A  Beam  B 


Figure  F.20:  Control  Input  Time  Response,  te(l. 00, 1.01s) 


164 


change.  In  terms  of  the  associated  control  input,  all  cases  once  again  exhibit  similar 
mean  response.  However,  dispersion  increases  at  shorter  time  cases  for  both  the  mean 
and  standard  deviation  control  input  responses,  which  is  due  to  more  instants  in  time 
of  control  authority  is  being  applied  during  the  simulation  duration. 
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Figure  F.22:  Switching  Time  Temporal  Average  Tracking  Error  Standard  Deviation 
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Figure  F.23:  Switching  Time  Temporal  Average  Control  Input 
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F.  3  Summary 

Controller  characterization,  based  on  weighting  matrix  and  switching  time,  has 
been  examined  using  SIMULINK  simulations.  Simulation  results  show  that  an  op¬ 
timum  regulator  performance  is  achieved  with  a  state  weighting  scale  factor  of  10 
and  control  weighting  scale  factor  of  0.1.  In  general,  initial  transients  occurred  ap¬ 
proximately  every  0.001s  due  to  the  correlation  time  for  the  bandlimited  white  noise 
blocks.  On  the  other  hand,  neither  switching  time  case  positively  or  negatively  af¬ 
fected  controller  performance.  For  each  case,  the  controller  adequately  compensated 
for  the  disturbance.  Hence,  for  the  purposes  of  this  study,  the  defined  optimum  regu¬ 
lator  weighting  parameters  will  be  implemented  with  a  switching  time  of  0.02s,  which 
allows  the  noise  correlation  time  to  remain  a  significant  order  of  magnitude  faster. 
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Appendix  G.  Dual  Fine  Tracking  Control  System  Simulation  Model 
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Figure  G.l: 


Dual  Fine  Tracking  Control  System  Simulation  Model 
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Figure  G.2:  Optical  Signal  Processing  Block 


Figure  G.3:  Plant  Disturbance  Noise  Block 
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Figure  G.4:  Measurement  Noise  Block 
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Figure  G.5:  Kalman  Filter  Block 


Figure  G.6:  Controller  Block 
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